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PREFACE AND ACKNOWLEDGMENTS 


The first experiment to combine airborne Doppler lidar and ground based 
dual Doppler-radar measurements of wind to detail the lower tropospheric flows 
in quiescent and stormy weather was conducted in Central Oklahoma during four 
days of June-July 1981. Data from these unique remote sensing instruments, 
coupled with those from conventional and novel in -situ facilities such as a 
500 m tall meteorologically instrumented tower, rawinson.de, and a dense 
network of surface based sensors, have been analyzed to enhance the under- 
standing of wind, waves, and turbulence. The data collected for this study 
and the analyses presented in this report had the multifaceted purposes of (1) 
comparing winds mapped by ground-based dual Doppler radars, by anemometers on 
a tall tower, and by NASA's innovative airborne Doppler lidar, (2) comparing 
measured atmospheric boundary layer flow with flows predicted by theoretical 
models; (3) investigating the kinematic structure of air mass boundaries that 
precede the development of severe storms and; (4) studying the kinematic 
structure of thunderstorm phenomena (e.g. downdrafts, gust fronts, etc.) that 
produce wind shear and turbulence hazardous to aircraft but are difficult to 
observe with conventional instrumentation. 

The results of experiments reported herein show that it is feasible to 
have detailed but remotely acquired measurements of the cloudless atmospheric 
motion in many provocative weather situations, whereas such observations were 
previously limited to scaled-down laboratory experiments. We have carried out 
these detailed measurements in nature's own laboratory so we can improve fore- 
casts using data that may become available to operational meteorologists. 
Meeting the challenge of improving weather forecasts and storm hazard warnings 
demands an integrated and well -conceived experimental design that builds upon 
a base of theoretical and experimental results and uses the technological 
expertise, leadership, and cooperation of several research groups, bound 
together with strong and enduring commitments. NASA's Marshall Space Flight 
Center (MSFC), NOAA's National Severe Storms Laboratory (NSSL), and the Uni- 
versity of Oklahoma's Department of Meteorology and Cooperative Institute for 
Mesoscale Meteorological Studies (CIMMS) have this commitment. This report is 
the first one showing the results of the combined efforts from this group 
dedicated to improve our understanding of weather. 
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The report is divided into 3 parts. Part I, "Intercomparison of Wind 
Data From Airborne Lidar, Ground-Based Radars and Instrumented 444 m Tower," 
gives results of the first intercomparison of atmospheric boundary layer wind 
measured by dual Doppler radars and airborne Doppler lidar. The second part, 
"The Structure of the Convective Atmospheric Boundary Layer as Revealed by 
Lidar and Doppler Radars" relates remotely sensed properties of homogeneous 
turbulent flow in the atmospheric boundary layer to those properties predicted 
by theories. Part III, "Doppler Lidar Observations in Thunderstorm Environ- 
ments," gives an analysis of the first observations with an airborne Doppler 
lidar of thunderstorm outflows and flow of environmental air into a complex of 
cumulus clouds. Appendix A describes the Doppler lidar instrumentation, its 
operation, and methods used to deduce the Doppler velocity component of hori- 
zontal wind. 

These first results are due to the dedicated efforts of many and the 
authors are indebted to NASA for its support, under the management of 
Drs. James Dodge and William Vaughan, and to both the NSSL and NASA staffs for 
their fine engineering support for this large data collection effort. Jim 
Bilbro and Ed Weaver of NASA, Chuck DiMarzio of Raytheon, and Bob Lee of 
Lassen Research, helped us to better understand the airborne lidar system, its 
unique features and the errors that it is vulnerable to. Dr. Dan Fitzjarrald, 
NASA, is to be especially thanked for presenting us with the lidar data in an 
easily usable form, and for constructive discussions. Drs. Douglas Lilly and 
Claude E, Duchon of the University of Oklahoma gave beneficial comments and 
suggestions during the course of this study. Encouragement and support from 
Dr. E. Kessler and staff at NSSL were critical to the quality of the work. 

Joan Kimpel of CIMMS and Robert Goldsmith of NSSL provided graphic services, 
and Michelle Foster of NSSL typed the manuscript. 

The work reported herein has already led to the publication of nine 
papers and to two others in preparation. The article "The structure of the 
convective atmospheric boundary layer as revealed by lidar and Doppler radars" 
has appeared in a recent issue of the international journal: Boundary-Layer 

Meteorology and another article "Estimation of the average surface heat flux 
over an inhomogeneous terrain" has been accepted for publication in that jour- 
nal. The paper "Comparison of winds, waves, turbulence as observed by air- 
borne lidar, ground based radars, and instrumented tower" has appeared in the 



American Geophysical Union's journal: Radio Science . The research on "An 

observational study of a mesoscale area of convection under weak synoptic- 
scale forcing" has been published by the American Meteorological Society's 
Monthly Weather Review . Papers to Applied Optics and The Journal of 
Atmospheric and Oceanic Technology are in preparation. This work has also 
supported two Master's Theses. 


Editor. Richard J. Doviak 
National Severe Storms Laboratory 
1313 Halley Circle 
Norman, OK 73069 
April 1985 
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PART I: INTERCOMPARISON OF WIND DATA FROM 

AIRBORNE LIDAR, GROUND-BASED RADARS AND 


INSTRUMENTED 444 m TOWER 



ABSTRACT 


On June 29, 1981 two ground-based Doppler radars, an airborne Doppler 
lidar, a very tall (444 m) instrumented tower, and a rawinsonde collected wind 
data in the Planetary Boundary Layer (PBL) in Central Oklahoma. This allowed, 
for the first time, intercomparison of wind fields synthesized from airborne 
Doppler lidar data with those from dual Doppler radar data. The vertical 
profile of wind in the PBL measured by the radars compared favorably with the 
profiles measured by the tower and rawinsonde while the one obtained from 
lidar data differed from the other three by as much as 3 m*s -1 in wind speed 
and 38° in direction. The time dependence of differences in wind estimates 
from radar and lidar suggested that these discrepancies could be attributed to 
a Schuler resonance in the aircraft's inertial navigation system which caused 
an erroneous component of the aircraft velocity vector to be subtracted from 
the lidar radial velocities, thus creating errors in the synthesized wind 
speed and direction. The vertical profile of turbulent fluctuations of the 
horizontal wind detected by the different sensing systems compared well. 

Also, spectra from the different sensing systems compared well in both 
magnitude and shape, suggesting that the lidar and radar detected similar 
turbulent structure. On July 2, 1981 an experiment was conducted to compare 
the point -by-point wind measurements made with the lidar, radar, and INS 
systems. The results of this experiment also suggest that Schuler resonance 
may account for discrepancies between the lidar and radar measurements. 



1.1. INTRODUCTION 


During June 29 - July 2, 1981, four research flights totaling 11 hours 
were made in Oklahoma by the NASA Convair 990, Galileo II aircraft, instru- 
mented with a coherent Doppler lidar. These flights were made within the 
range (i.e. < 100 km) of the NSSL Doppler radars and most of the time within 
the area of optimum dual Doppler analysis for the radars (one at Norman and 
the other Cimarron radar located at Page Airfield 40 km northwest of 
Norman). The data collected during these flights besides giving information 
on the structure of the cloud-free boundary layer, the prestorm environment, 
and the airflow in and around the clouds, afforded an opportunity for the 
first time intercomparison of wind data from totally different systems -- 
airborne Doppler lidar, ground-based Doppler radars, rawinsonde, and instru- 
mented tower. For this study, the principal instruments used in the data 
collection process were the Airborne Doppler Lidar and Ground-Based Doppler 
Radar whose operation and performance are described in Appendix A. In this 
part intercomparison data collected on June 29 and July 2 will be discussed. 
1.2 THE EXPERIMENT 

On June 29, 1981 data were collected between 1150-1410 CST. At this time 
skies were mostly clear and winds were moderate (5-8 m*s"^) out of the south 
and tended to back (turn counter-clockwise) and increase during the hours of 
data collection. Real time displays of Doppler velocities measured by radar 
indicated that the wind field was relatively uniform to ranges of at least 
80 km in all directions which was consistent with hourly surface observations 
during the time of data collection. There was significant daytime heating 
(maximum temperature of 35°C) and abundant gulf moisture, making the atmo- 
sphere convectively unstable; the height of the capping inversion was approxi- 
mately 1140m. Doppler lidar radial velocity estimates are calculated in 
resolution volumes spaced 320 m in range from the aircraft out to approxi- 
mately 15 km in range (see Appendix A for complete description). However, on 
this day attenuation of the lidar beam by atmospheric water vapor only allowed 
accurate velocity estimates to 5 km whereas on drier days a range of 10 km was 
achieved. Velocity vectors are calculated at intersections of fore and aft 
directed beams (Fig. 1.1) in a similar way in which velocity vectors are 
synthesized from the two Doppler radars (Doviak and Zrnic', 1984). 
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Figure 1.1 Schematic of Airborne Doppler Lidar system showing intersections 
of forward and aft lidar beams where wind vectors relative to the ground are 
calculated. Also shown are the other sensing systems that were used for data 
collections on June 29, 1981. 

The lidar instrumented aircraft was flown in oblong "race tracks" 
approximately 65 km x 15 km, southwest of Norman (Fig. 1.2) where wind 
observation by the two NSSL Doppler radars are relatively free of ground 
clutter interference (Doviak and Berger, 1980). Complete circuits were made 
at four different heights: 900, 750, 600, and 400 meters above ground 

level. The aircraft was flown in a counter-clockwise direction so that the 
Doppler lidar, which is located on the left side of the aircraft, would scan 
the same approximate area when going in opposite directions. However, on this 
day targets were only detected out to 5 km; so no intercomparison of lidar 
data collected just minutes apart was possible. All in all, eight data runs 
were made (each of approximately 5 minutes duration) during the period 1300- 
1410 CST. 
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Figure 1.2 Location of two Doppler 
radars, the instrumented tower T, 
the rawinsonde release point R, the 
approximate lidar instrumented 
aircraft track +, the grid area 
(square box) to where 1300-1302 CST 
data were interpolated. 


During this period, NSSL's two Doppler radars had their range resolution 
set at 450 m and were scanning the same volume in which the aircraft was being 
flown. Their Nyquist velocities were about 11 m*s _1 and their unambiguous 
ranges were 345 km (Appendix Table A.l). 

Every 5 minutes each radar scanned a volume subtended by 60° azimuth and 
3.4° elevation. This was done by scanning the radar 60° azimuthal ly, starting 
at 0.4° elevation angle and incrementing the elevation angle 0.5° at the end 
of each azimuthal sector. The Doppler radars are able to detect returns from 
refractive index fluctuations (cause by turbulent mixing of gradients of 
temperature and moisture) in the clear air boundary layer. Doviak and Jobson 
(1979) were able to observe clear air wind perturbations using these Doppler 
radars. Doviak and Berger (1980) were able to follow perturbations in the 
wind field as they were advected downwind by the mean flow. When data from 
consecutive radar scans (approximately 3.5 minutes apart) were compared, it 
was noticed that perturbations had moved downwind by a distance that the mean 
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wind would have advected them during this time. They concluded that velocity 
estimates synthesized from NSSL's dual Doppler radars are accurate enough (a 
few tenths of a ms -1 ) to map the wind field and its kinematic properties in 
the clear air out to ranges of 60 km. Sixteen volume scans were recorded from 
each radar (1250-1400 CST). 

Data were also collected throughout June 29th by instruments on a tower, 
operated and maintained by NSSL. This 444 m tall tower is located 39 km north 
of Norman. On this day wind speed, wind direction, temperature and vertical 
velocity data were collected at seven different heights between the surface 
and 444 m. 

A rawinsonde was released at 1150 CST from Tinker Air Force Base which is 
located approximately 30 km north-northeast of Norman. Data from the 
rawinsonde allowed us to determine the height of the inversion as well as the 
wind profile throughout the boundary layer. 

1.3. SYNTHESIS OF WIND FIELDS 

In order to compare data from the airborne Doppler lidar and the Doppler 
radars, data from both were interpolated to common horizontal Cartesian grids 
(20 km x 20 km, grid spacing 500 m) using a Cressman weight (Cressman, 1959) 
with horizontal radius of influence of 1.5 km. A vertical radius of influence 
of 300 m was used for the radar data but no vertical interpolation was needed 
for the lidar data because data were collected along beams on nearly 
horizontal surfaces. 

Before interpolation, data from both sensing systems were edited using 
reflectivity, spectral width, and signal -to-noise ratio thresholds. This 
removed data that were erroneous due to weak signal or other measurement 
uncertainties. The Doppler radar data were also edited along each radial to 
remove erroneous data points caused by point targets (aircraft, birds, 
etc.). This was done using an algorithm that compared the velocity value in 
question with the eight nearest points along the radial. If the velocity 
value was more than 3.5 nrs -1 from this mean, the velocity value was not used. 

Data from each radar (interpolated to common grids) were vectorially 
combined to arrive at a velocity vector at each grid point. In a similar way, 
data from the airborne lidar forward and aft directed beams were vectorially 
combined. 
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Grid volume locations were determined by plotting the Inertial Navigation 
System's (INS) derived track of the NASA aircraft. The gridded volumes were 
placed such that data collection by both Doppler radars and the lidar were 
acquired together. All in all, nine grid volumes were chosen to compare lidar 
and radar data. Location of the grid to where the 1300-1302 data were 
interpolated is shown in Fig. 1.2. 

1.4. COMPARISON OF WIND DATA 
1.4.1. Mean Wind 

Fig. 1.3 is a sample of comparable vector wind fields from both the lidar 
and radars at about 1400 CST. Mean wind given at the top of the figures is 
the vector average over the entire grid. Differences in mean wind speed of 
2.4 ms -1 and direction of 16° can easily be noticed. Comparison of mean wind 
for all of the nine common grid volumes is shown in Table 1.1. A systematic 
difference between the lidar and radar mean winds can be noticed. Lidar 
detected winds at earlier times and higher heights were lighter and had a more 
westerly component while at later times and lower heights the lidar detected 
winds were stronger and had a more easterly component than the radar detected 
winds. 
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Figure 1.3 Comparison of 1400 CST radar and lidar vector wind fields (the 
radar wind field is on the left). 





TABLE 1.1 


Comparison of Lidar and Radar Estimated Mean Wind 
(north and west distances from NRO are positive) 


Lidar Radar 


TIME 

(CST) 

Grid Origin 
(from NRO, km) 

Height 

(km) 

Vel 

(m/s) 

Di r 
(deg) 

Vel 

(m/s) 

Di r 
(deg) 

1300 

(-40, +30) 

.90 

6.8 

197.9 

6.6 

185.7 

1310 

(-30, +40) 

.88 

6.8 

218.6 

7.8 

180.6 

1320 

(-30, +25) 

.79 

7.4 

199.4 

8.0 

177.7 

1325 

(-15, +35) 

.78 

7.0 

204.9 

8.3 

181.1 

1330 

(-40, +40) 

.76 

7.1 

192.8 

8.4 

177.0 

1335 

(-50, +35) 

.58 

8.7 

169.4 

8.5 

173.6 

1340 

(-25, +25) 

.61 

9.5 

165.8 

8.9 

176.1 

1355 

(-45, +25) 

.40 

11.4 

156.5 

9.0 

171.3 

1400 

(-20, +20) 

.44 

11.6 

158.4 

9.2 

174.5 


1.4,2. Wind Profiles 

Another way of comparing data from the two sensing systems is to compare 
the vertical profile of wind that the two systems measure. Wind profiles can 
be constructed from lidar and radar data using the mean vector winds listed on 
Table 1.1 because these were computed at different heights throughout the 
boundary layer. The wind profile in the boundary layer as measured by lidar, 
radar, rawinsonde, and tower is shown in Fig. 1.4. The radar wind profile is 
nearly linear and also nearly constant with height, as is the rawinsonde 
profile. The lidar wind profile shows more vertical wind shear and non- 
1 inearity. 

Because the profile data were collected over an hour time period (the 
data near the top nearly an hour before the data at the bottom), changes in 
the wind speed and direction during this time would cause the wind profiles to 
be skewed. From the tower data continuously collected at the 444 m height, we 
detected a trend in both wind speed and direction between 1300-1400 CST. The 
wind speed increased 1.3 ms" 1 during this time and the wind direction had a 
counter-clockwise shift (i.e., backing) of 13°. Using the observation that 
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Figure 1.4 Comparison of wind 
profiles measured by lidar, radar, 
tower, and rawinsonde. 


the wind field was horizontally homogeneous, we assumed that this time trend 
observed by the tower occurred throughout the lowest kilometer of the atmos- 
phere in the air mass probed by the lidar and radars. We then adjusted the 

wind profiles relative to 1330 CST to remove this time trend (Fig. 1.5). The 
tower profile was constructed from a 12-minute average of data from 
1321-1333 CST. The rawinsonde profile was adjusted using the two-hour trend 
from the tower data. These wind profiles are now estimates of the actual 

vertical profiles as detected by the different sensing systems at the common 

reference time. The trend removal does not change any of the differences 
between the lidar and radar profiles, but it shows that the radar profile 
compares better in wind speed and direction with both rawinsonde and tower 
profiles than it does with the lidar profile. 

1.4.3. Intensity of Turbulence 

2 2 

The variances o u and o v of the orthogonal wind components (u is in the 
direction of the mean surface layer wind) were computed for each lidar and 
radar-estimated vector wind field, and then combined, <jy = (o^ + o^) to 
get the standard deviation of the horizontal wind velocity fluctuations. The 
total variance is 
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( 1 , 1 ) 


where a £ is the standard deviation due to errors in the velocity estimates and 
c? t is the standard deviation caused solely by turbulence and small scale 
organized flows such as waves and convective cells (Doviak and Zrnic', 

1984). The variances in (1.1) relate to the velocity estimated at grid points 
and includes errors introduced by interpolation and variance reduction due to 
the Cressman interpolation scheme (Doviak, et al., 1976). 

We estimated the radar error standard deviation, <j £ , by assuming that 
errors were due only to the radar's electronic noise. From signal -to-noise 
ratio measurements we calculated radial velocity estimate standard deviation 
using a technique described by Doviak and Oobson (1979). Following the 
methods of Doviak et al., (1976), we computed about a 1.1 nvs"* standard error 
for the single radial velocity estimate from each radar. However, about 25 



Figure 1.5 Comparison of wind 
profiles measured by lidar, radar, 
tower, and rawinsonde. Profiles 
have been adjusted to remove a time 
trend relative to 1330 CST. 
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data from each radar were interpolated to each grid point using a Cressman 
weighted filter function. Doviak et al., (1976) have shown that a spherically 
symmetric Cressman influence region reduces the variance of interpolated data 
by an amount. 


R 


2 


1.69 


( 1 . 2 ) 


where 1 - R 2 is the ratio of the grid point radial velocity variance a ^ to 

the variance of each radial velocity estimate used in the Cressman filter. N 

is the number of velocity estimates interpolated to a grid point and although 

the influence region used here is ellipsoidal we assume that (1.2) approxi- 

2 

mates the variance reduction. Thus, the error variance of the interpolated 

radial velocity is approximately 7.4 x 10“^ m^s"^. Now the variance o' of 

£ 

horizontal wind speed is a function of wind direction and for the case under 
consideration here, in which the wind intersects the line connecting the 
radars by about 45°, o £ is nearly equal to = 0.3 m s" 1 (Doviak et al., 
1976). All standard deviations larger than this value are attributed to 
actual variations in the windfield. 

An estimation of the lidar error standard deviation was made assuming 
that errors were caused by: (1) electronic noise in the lidar system which 

was estimated from the signal -to-noise ratio (Lee, 1980); (2) errors in ground 
speed and ground track angle as given in specifications for the Inertial 
Navigation System (INS); and (3) errors in the line-of-sight angle of the 
laser beam. Specified accuracies of the aircraft ground speed and the ground 
track angle are approximately ±0.5 m*s -1 and ±0.4° respectively. Pointing of 
the laser beam is accurate to within ±0.3° (Raytheon, 1983). Error estimates 
calculated from the signal -to-noise ratio averaged 0.35 m*s"l within 5 km of 
the aircraft. If we assume independent errors, the SD of the lidar radial 
velocity estimates due to errors is equal to 1.2 m s _1 sec. After the wind 
vectors are synthesized with approximately 35 radial velocity estimates inter- 
polated to each grid point, the standard deviation caused by errors in the 
velocity estimates is equal to 0.6 m s"*. 

Figure 1.6 shows the comparison of the standard deviation of velocity 
fluctuations measured by the three sensing systems. Standard deviations for 
wind measured by the tower instruments were averaged for approximately 
10 min. The total standard deviation o^. detected by the radar is nearly con- 
stant with height, although noticeable time trends are evident. From 1300 
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Figure 1.6 Standard deviation of 
the horizontal velocity fluctuations 
from lidar, radar, and tower. 


to 1320 CST there is a general decrease in standard deviation with time, while 
from 1320 to 1335 there is a general increase. From 1335 to 1400 CST no 
apparent trends are noticed. 

The total standard deviation measured by the lidar compares well in 
magnitude with that measured by the radars. Neglecting variations with 
height, a time trend similar to that of the radar Cy is noticed from 1300 to 
1320 CST (decrease in standard deviation with time), while from 1320 to 1400, 
little variation is noticed. 

The total standard deviation calculated from the tower data shows no 
apparent trends and tends to be larger than that estimated by either the lidar 
or radars. Two reasons may be responsible for these larger values. (1) the 
tower is close to the surface, where velocity fluctuations may be expected to 
be larger because of surface effects (Mason and Sykes, 1980), and (2) the 
tower instruments make point measurements while the lidar and radar both 
average through a resolution volume, thus, in effect, smoothing the data and 
decreasing the estimate of Oy. 
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1.4.4. Spectra of Horizontal Velocity Fluctuations 

Spectra of horizontal velocity fluctuations of the u and v wind compo- 
nents were computed from the lidar, radar, and tower data. For comparison of 
radar and tower data, spectra were calculated in the direction of the mean 
wind while for comparison of lidar and radar data, spectra were calculated in 
the direction of the flight path (Fig. 1.7). 

Spectra were calculated from radar data by interpolating data to a line 
of 32 grid points spaced 500 m apart. After the mean and trend were removed, 
Fourier analysis was done to resolve wavelengths between 1 and 16 km. Twenty 
spectra were averaged to give the spectra shown in Figs. 1.8 and 1.9 that were 
used for comparison. 

Spectra were calculated in a similar way from the lidar data. Several 
lines of 32 grid points spaced 500 m apart and parallel to the flight path 
were selected. Data were interpolated to these grid points and after the mean 
and trend were removed, Fourier analysis was done to resolve wavelengths 
between 1 and 16 km. Six spectra were averaged to give the spectra shown in 
Fig. 1.8. Spectra from lidar and radar data compare well in both magnitude 
and shape, especially in the v-component where peaks in both spectra occur at 
about the same wavelength. 



Figure 1.7 Schematic showing the 
direction along which spectra were 
calculated to compare data. 
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SPECTRA CALCULATED IN DIRECTION OF FLIGHT PATH SPECTRA CALCULATED IN DIRECTION OF FLIGHT F 



Figure 1.8 Comparison of spectra of the u and v velocity component 
fluctuations from lidar and radar data. Spectra are for scales 
(wavelengths A) along directions parallel to the flight path. 


Figure 1.9 is the comparison of spectra in the direction of the mean wind 
from tower and radar data. Spectra were calculated from the tower data 
assuming that all perturbations in the wind field were advected with the mean 
wind and did not change in the time frame needed to get enough data for 
spectral analysis (Taylor's hypothesis). Thus, a line of grid points was set 
up in the direction of the mean wind and radar data were interpolated to 
them. Likewise, radar data interpolated to lines parallel to the mean wind 
direction were used in the comparison shown in Fig. 1.9. Radar and tower 
spectra compare well in both magnitude and shape. 

1.5. EXPLANATION OF DIFFERENCES BETWEEN LIDAR AND RADAR MEASURED WINDS 

Two factors may be responsible for the differences between the lidar and 
radar measured mean winds (Table 1.1) and the difference between the lidar 
wind profile and the other wind profiles (Fig. 1.5). These are: (1) the 
lidar has a very small vertical resolution (20 cm) which may have allowed it 
to detect small vertical scale variations in the wind field that were not 
detected by the radar because of its larger vertical resolution (-700 m) or by 
the tower or rawinsonde because of time and location differences, or (2) an 
inherent error in the Doppler lidar system or in the INS of the aircraft 
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Figure 1.9 Comparison of wind spectra measured by tower instruments and 
radar. Spectra are for scales A along the direction of the mean wind. 


(which is used to subtract out the motion of the aircraft from the lidar 
measured radial velocities) may have caused the lidar measured winds to 
oscillate about the actual mean winds. 

In order to test whether the different vertical resolutions of the 
sensing systems may have been the cause of differences in the estimated wind 
profiles, we constructed comparable wind profiles from the lidar and radar 
data. This was done using radar data as close as possible to the radars to 
construct a wind profile with the highest vertical resolution possible. Data 
about a 20 km range were used to construct a wind profile, because ground 
clutter overwhelmed air velocity data from closer ranges. The vertical reso- 
lution of the radars at this range is approximately 280 m. Wind profiles at 
20 km range were constructed by finding a representative velocity value for 
each elevation angle at range 20 km. The representative velocity was calcu- 
lated by finding the modal velocity (the average was found to be biased 
towards zero due to ground clutter effects) for the data points in a 4 km x 
10° sector. Hence, two radial velocity profiles were constructed (one from 
each radar) using 1330 tilt sequence data. Assuming horizontal homogeneity, 
we synthesized these two profiles arriving at mean wind speed and direction 
profiles. 
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The lidar profile was altered for better comparison by taking a vertical 
average of the trend-removed profile (valid at 1330 CST) over a 280 m vertical 
distance centered at the heights where radar data were available. Even with 
the same vertical resolution there were still differences as large as 

1.8 m*s"l in wind speed and 18° in wind direction (Fig. 1.10). Thus, we 

concluded that the different vertical resolution of the systems cannot explain 
the large differences between the 1 idar-detected and the radar-detected wind 
profiles. 

We then turn to the possibility of an error in either the lidar system or 

in the aircraft's INS. Looking at Fig. 1.4 it is evident that winds measured 

by the lidar between 1300-1330 CST had a more westerly component and were 




l 1 — 1 I 
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WIND DIRECTION (deg) 


Fi oon 6 1,10 I rend removed lidar wind Profile valid at 1330 CST (averaged over 
a 280 m depth) and radar wind profile (vertical resolution 280 m) at 1330 CST. 



lighter than winds measured by the radars. Between 1330-1400 CST the lidar 
measured winds had a more easterly component and were stronger than winds 
measured by the radar. 

To investigate reasons for this apparent bias of the lidar winds, we 
calculated the correction that had to be added to the forward and aft radial 
velocities to make the lidar and radar-measured winds consistent. The 
corrections to the forward and aft radial velocities that were needed at 
different times throughout the hour of data collection are given in 
Table 1.2. Unexpectedly, it was found that in all cases the correction needed 
was nearly the same for both forward and aft radial velocities. The only wind 
vector that could cause equal addition to both radial velocities is a vector 
perpendicular to the aircraft path. Hence, we plotted the corrections perpen- 
dicular to the aircraft track vs. time, setting corrections to the west nega- 
tive and to the east positive (Fig. 1.11). Also plotted in Fig. 1.11 is a 
sinusoidal wave with amplitude 4 nrs"* and period 84 minutes, the period asso- 
ciated with the Schuler resonance, which is an inherent source of error in an 
INS (Frye, 1958; Britting, 1971). The calculated corrections follow this 
curve remarkably well. 

TABLE 1.2 

Corrections that needed to be added to forward and aft radial 
velocities to make the lidar measured mean wind consistent with 
the radar measured mean wind. 


Time Correction 

Forward 

1300 CST -1.45 m/sec 


1310 

+4.75 

1320 

-2.75 

1330 

+2.49 

1340 

+1.50 

1350 

-2.67 

1400 

+3.52 


Ai rcraft 


Aft 

Heading 

-1.87 m/sec 

20° 

+4.27 

205° 

-2.81 

0 

O 

CVJ 

+2.09 

ro 

o 

o 

o 

+1.72 

25° 

-2.63 

200° 

+3.48 

25° 


To support the possibility of a Schuler resonance, we examined 
photographs that were taken from the aircraft every 10 seconds looking 



straight down. We were able to locate certain landmarks found on the ground 
and compare these locations with that given by the INS. We found that the 
aircraft location given by the INS was consistently offset by -15 km to the 
northeast of its actual location (Fig. 1.12). When this offset was removed, 
we found that the actual positions of the aircraft oscillated about the INS 
locations in a way consistent with the position curve (integral of the velo- 
city curve) associated with the Schuler resonance shown in Fig. 1.11. Table 
1.3 is a comparison of position errors of the INS (after the offset was 
removed) vs. the estimated position error due to the Schuler resonance. We 
zeroed the difference of data in the first of Table 1.3. 

The agreement seen in Table 1.3 corroborates the evidence presented 
earlier, leading us to believe that a Schuler resonance of the aircraft's INS 
caused the subtraction of erroneous components of the aircraft's ground- 
relative velocity vector from the lidar measured forward and aft radial 
velocities resulting in errors in the synthesized mean wind speed and 
direction. Furthermore, the maximum position error of the Schuler resonance 
shown in Fig. 1.11 is 3.2 km (disregarding the offset position error which may 
be due to other sources) after 4 hours of flight from California and is well 
within the stated accuracy of the INS of 3.2 km per hour of flight (NASA, 
1979). 


TABLE 1.3 

Difference between actual and INS positions compared with 
estimated difference due to the Schuler resonance from 
Fig. 1.11. 


Position 


Time 


Actual Position Error 
Estimated From Photos 

Error Estimated 
From Schuler Resonance 

1358 

CST 

+0.7 (km) 

+0.7 (km) 

1338 


-3.1 

-3.0 

1320 


-0.6 

-1.6 

1.6. 

COMPARISON 

OF DATA COLLECTED ON JULY 2, 1981 



On July 2, 1981, the CV-990 aircraft circled the Norman Doppler radar at 
a range of about 50 km, at 1.1 km height (Fig. 1.13) in a counter-clockwise 
direction. This was done so that the Airborne Doppler Lidar System (ADLS), 


VELOCITY (m/s) 



TIME (CST) 


Figure 1.11 Plot of aircraft 
movement (undetected by the INS) 
perpendicular to the aircraft 
heading vs. time. Velocities and 
distances to the northwest are 
negative and to the southeast 
positive. 



Figure 1.12 Map showing difference 
between actual aircraft location and 
INS indicated location. 
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which was located on the left side of the plane, could collect data on the 
side of the plane nearest to the radar where Doppler radar velocity estimates 
were not contaminated by echoes of the plane. The 1 idar-instrumented aircraft 
collected data from 1353-1438 CST, while the plane was flown in a one and one 
half circle around the Norman Doppler radar. 

During this time, data were also collected in two different modes by the 
Doppler radar as it tracked the aircraft. Data were collected in the time 
series mode from 1350-1439 CST (i.e., in-phase and quadrature-phase echo 
samples at 16 range gates are recorded) and in the pulse pair mode (i.e., real 
time estimates of radial velocity are made by the pulse-pair processor and 
recorded in addition to the echo samples of I and Q) from 1423-1439 CST. How- 
ever, the pulse pair velocity estimates were found to be in error because of 
contamination by overlaid echoes and thus were not used for this comparison. 

A rawinsonde was released at 1150 CST on this date from Tinker Air Force 
Base, located 30 km NNE of Norman. The wind profile in the lowest 2 km of the 
atmosphere as detected by the rawinsonde is shown in Fig. 1.14. Wind measured 
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Figure 1.13 Flight path of the ADLS 
on July 2, 1981. Tinker AFB is the 
rawinsonde release site. 
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RAWINSONDE WIND PROFILE 
1152 CST, 2 JULY 81 



(deg) (m/sec) 


Figure 1.14 Wind profile as 
detected by rawinsonde at 1152 CST 
July 2, 1981. 


by the rawinsonde at the height at which data collection by the lidar and 
radar was made, was approximately 2 nrs"* from 55°. 

In order to compare data from the lidar and radar, lidar data from the 
forward and aft-pointed beams were interpolated to separate (but commonly 
located) Cartesian grids using a Cressman weight with radius of influence of 
500 m (Cressman, 1959). The radial velocities of the two grids of data were 
then vectorially combined to determine wind speed and direction, and then the 
component of this wind towards the radar was computed. All lidar data within 
a 2° azimuthal sector from the radar were then averaged, and assumed valid at 
the mid-point, to enable comparison with the radar data. 

The radar-measured radial velocity values were determined from averages 
of 100 Doppler spectra (Fig. 1.15). Spectra at the same range were averaged 
over an azimuth sector of approximately 10° (~20 seconds of data) and the peak 
of this average spectra was used as the radial velocity value representative 
for that azimuthal sector at that range. The width o (the square root of the 
second moment) of the Doppler spectra is 2 m s“*. However, this should not be 
taken to be indicative of the velocity fluctuations within the beam because 
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VELOCITY (m/s«c) 


Figure 1.15 Example of an average 
of 100 Doppler velocity spectra 
showing approximate 20 dB signal -to- 
noise ratio in the spectral power 
densiti es. 


window bias needs to be removed (Waldteufel, 1976). The velocities in sixteen 
sector volumes (spaced 600 m in range) were then averaged together to arrive 
at a mean radial velocity for a volume approximately 10° (azimuth) by 10 km 
(range) by 0.7 km (vertical). 

Winds derived from the aircraft's INS were also available every minute 
during the time the aircraft was flown. The component of the INS-derived 
winds towards the radar was used for intercomparison. 


Fig. 1.16 is a plot of the radial velocities from the three different 
sensors through a full 360° arc about the Norman Doppler radar. The points 
are connected for visual clarity and do not represent interpolation. The 
lidar and INS-derived winds have both been converted to radial velocities 
towards the radar so that these winds can be compared to the radial velocity 
estimates made with the Doppler radar. The mean wind detected by the three 
sensing systems is approximately out of the east at 5 m s"V. 

Although we can expect deviations about the sine wave because of non 
uniformities in the wind field. Fig. 1.16 shows that the deviations found in 
any one sensor's velocity estimates are not consistent with those found in 


s 


22 




Figure 1.16 Comparison of radial velocities detected by radar, lidar, and the 
aircraft's INS. 

either of the other sensors' estimates over more than a few tens of degrees. 
Possible reasons for this disparity between sensors are numerous. For one, 
the Doppler radar velocity estimates are inconsistent and biased between 180°- 
300° due to overlaid second trip echoes (Doviak and Zrnic', 1984) of storms 
that were 100-180 km west of Norman (Fig. 1.17). The apparent departure of 
the radial velocities from a sine curve expected for uniform wind (e.g., 
velocities in the azimuthal sector 180 - 0 deg, reach values about 3 m s“* 
whereas velocities in the sector 0 - 180 deg. reach a peak value of about 
4 m s" 1 for the radar and INS data and almost 6 m s" * for the lidar data), 
measured by all of the sensing systems could be due to perturbations in the 
wind field caused by storms which were 50 km west of the data collection 
circle. Furthermore, the resolution volume over which the radial velocities 
are estimated differs vastly for the lidar and radar measurements. The INS 
measurements are point values smoothed by the response of the aircraft. And 
finally, the lidar measures the velocity of air weighted by the reflectivity 
distribution of aerosols, whereas the radar measures air velocity weighted by 
the reflectivity distribution of refractive index irregularities or by 
insects, if they contribute to target reflectivity. 
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Figure 1.17 Photo of NSSL's WSR-57 radar display (1427 CST , July 2, 1981) 
depicting storms to the west which caused second trip echoes in the Doppler 
radar data. Range rings are spaced 40 km. 
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It can be seen from Fig. 1.16 that the INS-derived radial velocities tend 
to fluctuate more than either the lidar or radar radial velocities, possibly 
because the INS velocities are point measurements whereas the lidar and radar 
velocities are averages (over an azimuth sector of about 2° and 10° 
respectively). The obvious bias between lidar and radar radial velocities 
between 0° and 180° could be possibly related to a Schuler resonance of the 
aircraft's INS, as was found in the lidar data collected on June 29, 1981. 

1.7. CONCLUSIONS 

A major objective of this experiment was to compare data from NASA's air- 
borne Doppler lidar and NSSL's dual Doppler radar system. Discrepancies in 
the mean wind and wind profile detected by the different sensing systems were 
explained as being caused by a Schuler resonance of the aircraft's INS, which 
caused an erroneous component of the aircraft's ground-relative velocity vec- 
tor to be subtracted from the 1 idar-measured velocities. Furthermore, it is 
believed that these errors were as large as 4.5 nrs -1 due to the lengthy time 
that the aircraft was in flight before collecting data. 

Standard deviation of velocity fluctuations measured by the lidar and 
radar compared well, as did spectra of the velocity fluctuations of lidar, 
radar, and tower data. Thus, it is concluded that the lidar and radars are 
accurate enough to map velocity perturbations in the clear-air boundary layer. 

Finally, it is concluded that NASA's airborne Doppler lidar system is 
capable of measuring wind fields in the clear air environment on a smaller 
scale than was previously available. An improved version of the airborne 
Doppler lidar system will be used to conduct research flights in the future. 
This improved lidar system will incorporate changes designed to correct all of 
the major error sources encountered in the flights conducted in 1981. Accele- 
rations of the aircraft, which are not detected by the INS (such as a Schuler 
resonance), will be accounted for by alternately pointing the lidar beam below 
the horizon so that ground returns can be observed and their apparent Doppler 
velocities removed in later analysis. This system will be used to observe 
regions surrounding convective phenomena which at present is not possible 
using conventional Doppler radars. The airborne Doppler lidar system will be 
especially useful when used in conjunction with Doppler radars because radars 
are capable of mapping the wind field inside a convective storm whereas the 
lidar maps the wind field just outside the storm. This should give us a 
better insight into convective phenomena. The airborne Doppler lidar system 
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will also be useful to examine turbulence and waves in the clear air planetary 
boundary layer, especially near the capping inversion where few observations 
have been conducted. 
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PART II: THE STRUCTURE OF THE CONVECTIVE ATMOSPHERIC 

BOUNDARY LAYER AS REVEALED BY LIDAR 
AND DOPPLER RADARS 
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ABSTRACT 


The turbulent structure of the convective atmospheric boundary layer, 
based on the analyses of data from the instrumented NSSL-KTVY tower, airborne 
Doppler lidar, and ground-based Doppler radars, is presented. The vertically 
averaged winds over the baroclinic boundary layer agree well with those for a 
barotropic atmosphere supporting the hypothesis of Arya and Wyngaard (1975) 
that vertically averaged winds are insensitive to barocl ini city. The computed 
momentum flux profiles are affected by barocl ini city. The momentum flux pro- 
files deduced from wind measured with radar agree well with those numerically 
predicted. Application of the assumptions that measured wind equals the geo- 
strophic wind V g at levels above the inversion and that the V g profile is 
linear yields surface pressure gradients consistent with those estimated from 
a network of pressure sensors. Horizontal wind spectra from lidar, radar, and 
tower data compare well with each other in shape and magnitude. A consistent 
peak near 4 km in all the computed spectra might have been caused by hori- 
zontally symmetric cells with a horizontal wavelength 4 times the boundary 
layer height as shown in Kuettner (1971) for the case of weak wind shear. 
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2.1. INTRODUCTION 

Our present knowledge of the structure of the convective atmospheric 
boundary layer comes mainly from the results of numerical models (Deardoff, 
1972, Wyngaard, et al., 1974). Although data from the Wangara expedition 
(Clarke, et al., 1971) and the Minnesota boundary layer experiment (Izumi and 
Caughey, 1976; Kaimal, et al., 1976) confirm some of the numerical model pre- 
dictions, there is still an urgent need for observational studies on the con- 
vective atmospheric boundary layer (ABL). 

Results on the structure of the convective ABL based on the analyses of 
data from the instrumented 444 m NSSL-KTVY tower, airborne Doppler lidar, and 
ground-based Doppler radars will be presented here. A brief description of 
the NASA airborne Doppler lidar system and the NSSL Doppler radars is given in 
Appendix A. 

The data collected on June 29, 1981 were found to be especially suitable 
for the study of the convective ABL structure. The skies during this data 
collection period were mostly clear and void of precipitation targets. Also, 
there was a well-defined inversion lid for making a reliable estimate of the 
height of the convective ABL, which is an important scaling length for the 
generalization of observational results on ABL structure. 

2.2. DATA COLLECTION 

On June 29, 1981 the 1 idar-i nstrumented aircraft was flown in oblong 
"race tracks" (Fig. 1.2 of Part I) and data collected for the intercomparison 
experiment described in Sect. 1.2 is used here to deduce the kinematic struc- 
ture of the ABL. The radar data were collected simultaneously as described in 
Section 1.2 and data from the instrumented NSSL-KTVY tower, located about 39 
km north of Norman, were also collected throughout June 29, 1981, A rawin- 
sonde was also released at 1150 CST on this day from Tinker Air Force Base 
which is located approximately 30 km north-northeast of Norman. 

2.3. DATA ANALYSIS 

Data from the radars, lidar, and tower were edited to screen out data 
that were obviously erroneous. Editing was done using signal -to-noise ratio, 
spectral width, and reflectivity thresholds, if data at a point exceeded any 
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of these thresholds, they were discarded. The Doppler radar data were also 
edited along a radial to remove erroneous data due to anomalous targets such 
as aircraft and birds. This edit procedure is fully described by Eilts 
(1983). 

2.3.1. Construction of the Wind Profile 

An observed wind profile (Fig. 2.1) was constructed using tower data 
below 444 m and dual Doppler radar data from 500 m up to 2 km. Lidar data 
were not used in the construction of the wind profile because of the errors in 
the mean wind detected by the lidar, possibly due to a Schuler resonance of 
the aircraft's INS which caused erroneous components of the aircraft's ground 
relative motion to be subtracted from the measured radial velocities (see Part 
I). The Cartesian coordinate system was rotated such that the x-axis and the 
u-component of wind were aligned along the direction of surface wind and the 
y-axis and v-component of wind were 90° counterclockwise to the surface 
wind. The mean wind data from the tower represent a time average of the data 


OBSERVED (h/L=-15.7) 

— - — calculated geostrophic wind profile 
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Figure 2.1 Observed wind profile 
and estimated geostrophic wind 
profile. The asterisks at the 
surface denote the surface 
geostrophic wind components 
estimated from surface pressure 
gradients. 
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collected between 1320-1335 CST, while those from the dual Doppler radars 
represent a space average of 1330 CST tilt sequence data that had been inter- 
polated to 20 x 20 km grids placed every 250 m from 500 m up to 2 km. The 
geostrophic wind profile is calculated using the methods described in section 
2.4.2 under the assumption that it is linear with height and equals the 
observed wind above (z=2 km) the inversion. 

2.3.2. Calculation of Velocity Spectra 

Spectra of the longitudinal and transverse components of wind were calcu- 
lated from the tower, radar, and lidar data. Spectra were calculated from the 
radar data by interpolating the data to a series of 32 grid points spaced 500 
m apart. After the mean and trend were removed, Fourier analysis was done to 
resolve wavelengths between 1 and 16 km. Thirty-two spectra were averaged to 
give the spectra shown in Fig. 2.2. 

Spectra from the lidar data were calculated in a similar way by interpo- 
lating data to a series of 32 grid points aligned along the flight path. Six 
spectra were averaged to yield the spectra shown in Fig. 2.3. 

Spectra from the tower data were calculated assuming Taylor's hypo- 
thesis. A series of 32 grid points spaced at 500 m was aligned along the 
direction of the mean wind, and data were interpolated to the grid points by 
assuming that all of the data were advected with the mean wind at 444 m. Four 
spectra were averaged to give the spectra shown in Fig. 2.4. 

2.3.3. Determination of Micrometeorological Parameters 

The tower data collected enabled us to determine u* the friction velo- 
city, Q q = (w'«i ’) the surface kinematic heat flux, and hence L the Obukhov 
length. Estimation of these parameters is necessary for comparing our obser- 
vational results with numerical model results. 

The surface kinematic heat flux Q 0 was estimated using a technique 
described by Sundara-Rajan et al., (1984) for highly convective conditions. 
This technique (Appendix B) is based on the assumption that the normalized 
vertical velocity variance remains constant between the levels 0.2h and 0.6h 
where h is the height of the convective ABL — an assumption adequately sup- 
ported by observations of the convective ABL over horizontally homogeneous 
sites (Kaimal et al., 1976). Sundara-Rajan et al., (1984) found that this 
method gives reliable estimates of the average surface heat flux even over an 
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Figure 2.3 Spectra of horizontal wind velocity fluctuations from lidar 
data. LOS angle is the line-of-sight angle of the lidar beam. The two 
spectra represent data from the forward and aft looking beams. 
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Figure 2.4 Spectra of horizontal wind velocity fluctuations from tower data. 

inhomogeneous urban area like Oklahoma City. Th us, it is straightforward to 
compute Q 0 from the vertical velocity variance (w 1 ) and the height (h) of the 
convective ABL using the assumption that 



1/3 

where w* is the convective scaling velocity given by w* = ( 3hQ Q ) 
where g = g/0 Q is the buoyancy parameter; g is the acceleration of gravity 
and 9 Q is the average virtual potential temperature of the surface layer. We 
had to resort to this indirect way of estimating Q 0 because the surface-layer 
temperature profile from the tower observations was not sufficiently accurate 
for the "profile" method (Appendix B). 

The friction velocity u*, used as a scaling parameter for our horizontal 
wind, was computed using an iterative technique assuming the following 
Businger-Dyer wind profile (Businger et al., 1971; Businger, 1973) for the 
convective surface layer: 


_ u * 

u = Y~ Un (z/z Q ) - <KS)] 


( 2 . 2 ) 


where 5 = z/L and z 0 is the roughness scale. Eq. (2.2) was applied to the 
lowest levels (7 and 26 m) of the tower, yielding 


j 


34 





(2.3) 


u* 

U 2 " U 1 = iT“ twi ^ z 2 /l 0 + + - *» U 2 )3 

where u^ - u^ are the mean wind speeds at the heights of z 2 and z\ 
respectively. The function f is given by 

* = 2tn [(1 + x ) / 2 ] + 2n [(1 + x 2 )/2] - 2 tan' 1 (x) + -n/2 
1/4 

where x = (1 - 15c) andq and c 2 are equal to Zj/L and z 2 /l respectively. 
The von Karman constant k is assumed to be equal to 0.35 in accordance with 
Businger et al., (1971). Eq. (2.3) was used iteratively, starting with the 
assumption that tUi) and \p(c ? ) = 0 (neutral stability). Approximately 10-15 
iterations were needed for convergence. The Obukhov length L = (-u* /k6Q Q ) 
was determined using the computed values of u* and Q 0 . 

The height h of the convective ABL was determined from the virtual 
potential temperature profile computed from the rawi nsonde data (Fig. 2.5). 
For the two hours of data collection on June 29, 1981, the parameters were 
estimated to be: 

Q q a 0.13 m s u* a 0.42 m s' 1 , L a -70 m, and h a 1140 m. 

2.3.4. Estimation of Momentum Flux Profiles 

From the observed wind and geostrophic wind profiles we can infer the 
profiles of momentum fluxes using the velocity defect equations: 


f(V-V)=du'w' 

9 

/dz 

(2.4) 

f (U - U) = d TiT 

/dz 

(2.5) 


where f is the Coriolis frequency and U g , V g are the geostrophic wind 
components. These equations attribute the ageostrophic component of the wind 
to momentum flux gradients which are produced by turbulence in the ABL. 

By integrating (2.4) and (2.5) with respect to z, using the boundary 

. O 

conditions -u V = u* and v V = 0, at z = 0, we arrive at the equations for 
the momentum fluxes at height z 

-uV (z) * / f [V (z) - V(z)]dz + u* 2 (2.6) 

o a 


j 


35 



HEIGHT (km) 


29 JUNE 1981 TINKER AFB 



Figure 2.5 Virtual potential 
temperature profile from rawinsonde 
data showing the height of the 
capping inversion. 


W (z) = j f[TT(z) - U (z)]dz (2.7) 

o 9 

Equations (2.6) and (2.7) were used to calculate the momentum fluxes at 
100 m intervals through the ABL. 

2.4. RESULTS 

2.4.1. Vertically Averaged Winds 

Using the observed wind profile, a vertical average over the whole ABL 
was computed for both components of wind. According to the hypothesis of Arya 
and Wyngaard (1975), the vertically averaged wind through the convective ABL 
can be regarded as independent of baroclinicity. To test this hypothesis, we 
also computed the vertically averaged wind components for the barotropic, 
convective ABL using the following relations for the similarity functions A-,- 
and B.: given by Arya (1977) based on the results of some recent, sophisticated 
numerical models for the convective ABL. 
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A.j - £n(-h/L) + £n(fh/u*) + 1.5 


( 2 . 8 ) 


B i = k(fh/u*) _1 = 1.8(fh/u*) exp (0.2h/L) 

(2.9) 

The similarity functions A^ and B n - are defined as 


A i = An ( h / z Q ) - k <U g >/u* 

(2.10) 

B. = -k <V g >/u* (sign f) 

(2.11) 

The values of A^ and B^ were first computed for the observed conditions of 
h/L, h/z Q , and fh/u* were determined from equations (2.10) and (2.11). The 
vertically averaged wind components <TJ>/u* and <V>/u* were computed from the 
following relations for them derived from the vel ocity -defect equations (2.4) 
and (2.5): 

<TT>/u* = <U g >/u* 

(2.12) 

<Y>/u* = <V g >/u* + u*/fh 

(2.13) 


The vertically averaged wind components and their counterparts for the 
barotropic case are given in Table 2.1. The agreement between the two sets of 
values is good, supporting the hypothesis of Arya and Wyngaard (1975) that the 
vertically averaged winds over the convective ABL can be assumed to be 
independent of barocl inicity. 

Unlike the vertically averaged winds, the vertical profile of wind is 
strongly affected by baroclinicity, as is evident from the numerical model 
results reported by Arya and Wyngaard (1975). Direct comparison of 
theobservational results on the vertical wind profile with the results of Arya 
and Wyngaard (1975) could not be made became of the difference in the values 
of the baroclinicity parameters = 77- (-—*■! and M = ~ (-r-^-1 , and also 

X U ^ d Z 0 jr ^ o 

the stability parameter h/L. 
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TABLE 2.1. Comparison of Vertically Averaged Winds 


<U>/u* 


<V>/u* 


Observed 

18.2 

-2.4 

From similarity functions 

18.0 

0.0 


2.4.2. Momentum Flux and Geostrophic Wind Profiles 

For the estimation of momentum fluxes from equations (2.6) and (2.7), the 
geostrophic wind components at the surface and their variation with height 
should be known. The surface geostrophic wind was estimated from surface 
pressure gradients. The actual wind above the boundary layer at the highest 
level (2 km) measured by the radars was assumed to be equal to the geostropic 
wind at that height; also the assumption of a linear variation of geostrophic 
wind with height was made. The computed momentum flux profiles are shown in 
Fig. 2.6. Also shown in the same figure are the momentum flux profiles com- 
puted from the higher order closure model of Wyngaard, et al. (.1974) for a 
barotropic convective ABL. There is considerable disagreement between the 
two, and the momentum fluxes computed from the radar-detected winds (hence- 
forth to be called "observed momentum fluxes") do not go to zero at the top of 
the boundary layer. This may well be due to the assumption of constant geo- 
strophic shear as well as the difficulty of estimating correctly the surface 
geostrophic wind components from surface pressure gradients. 

It can be seen from equation (2.7) that for v'w' to go to zero at the top 
of the boundary layer. 


<U> = <U g > (2.13) 

where the angle brackets denote the average over the whole boundary layer. 
Figure 2.6 shows clearly that this condition is not satisfied, indicating the 
possibility of inaccuracies in the computation of the geostrophic wind compo- 
nents. 

Let us assume that the computed geostropic wind components at the surface 
are accurate. It is not necessary to make any assumption regarding the 
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Figure 2.6 Momentum flux profiles 
(dashed lines) estimated from 
velocity defect equations using a 
geostrophic wind profile linear from 
the estimated surface geostrophic 
wind to the actual wind at 2 km. 
These flux profiles are compared to 
the numerical model results of 
Wyngaard, et al. (1974) for the 
barotropic, convective boundary 
layer. 

o 

NORMALIZED STRESS 

geostrophic wind components above the top of the boundary layer if we make use 
of equation (2.3) and the following equation derived from equation (2.6): 

<V> = <V g > + u* 2 /fh (2.14) 

for u V to be zero at the top of the boundary layer. But still the 
assumption of uniform geostrophic shear has to be made. As before, the 
momentum flux profiles could be computed from equations (2.6) and (2.7), and 
the results are shown in Figure 2.7. The agreement with the flux profiles 
reported by Wyngaard, et al. (1974) is better than that for Fig. 2.6. The 
observed momentum flux profiles show a shift to the right of the corresponding 
profiles for the barotropic case, as would be expected from the results of 
Arya and Wyngaard (1975) for the observed barocl ini city 
(3U g /3z = -6 m/sec/km, 3V g /3z = -4 m/sec/km). 

We can compute another set of momentum flux profiles without using the 
computed surface geostrophic wind components by assuming that the actual wind 
above the boundary layer at the highest level (2 km) detected by the radars 
equals the geostrophic wind at that height, together with the assumption of 
al inear variation of geostrophic wind with height. Figure 2.8 shows the 
computed momentum flux profiles based on these assumptions, and the agreement 
with the results of Wyngaard, et al. (1974) is much better, indicating that 
our assumption that the geostrophic wind equals the actual wind at the highest 
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Figure 2.7 Momentum flux profiles 
computed using a geostrophic wind 
profile linear from the estimated 
surface geostrophic wind with the 
constraint that the momentum fluxes 
go to zero at z * h. 
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Figure 2.8 Flux profiles assuming a 
linear profile of geostrophic wind, 
for which the geostrophic wind 
equals wind at 2 km, and the 
constraint that the fluxes go to 
zero at z = h. 


0 




level measured by the radars may be fairly accurate- The slight shift of the 
observed momentum flux profiles to the right of the numerical model results is 
again in the correct sense for the observed barocl inicity. 

2.4.3. Spectra of Horizontal Wind 

Space and time averaged spectra from radar and tower (see Fig. 1.9 of 
Part I) compare very well both in magnitude and shape, as did the average 
spectra (averages of spectra along 6 lines 500 m apart) from radar and lidar- 
estimated wind fields (see Fig. 1.8 of Part I). A consistent peak, although 
weak, near 4 km wavelength can be noticed in Figures 2.2, 2.3, and 2.4. This 
peak was also noticed in horizontal wind spectra at different heights through- 
out the period of data collection, suggesting the prevalence of some time- 
consistent three-dimensional structure that is symmetrical in the horizontal 
directions with a horizontal wavelength of approximately 4 km. It is proposed 
that this three-dimensional structure is the result of symmetric cells that 
have a wavelength 4 times the boundary layer height (h = 1140 m during the 
period of this study), as shown by Kuettner (1971) for the case of weak wind 
shear. 


2.5. CONCLUDING REMARKS 

The foregoing study of the structure of the convective ABL indicates that 
the vertically averaged winds over the boundary layer can be assumed to be 
independent of barocl inicity in accordance with the hypothesis of Arya and 
Wyngaard (1975). Our assumption that the geostrophic wind equals the actual 
wind measured above the top of the boundary layer along with a geostrophicwind 
profile that is linear throughout the boundary layer seems to yield momentum 
flux profiles which compare well with the numerical model prediction of 
Wyngaard, et al. (1974) and Arya and Wyngaard (1975). Likewise, this 
assumption yields surface geostrophic winds that are fairly consistent with 
those deduced from surface pressure gradients as can be seen in Fig. 2.1. 

Spectra of the longitudinal and transverse components of the wind 
computed from lidar, radar and tower data compare well both in shape and 
magnitude. A consistent peak near the 4 km wavelength found in all the 
computed spectra is proposed to be caused by horizontally symmetric cells with 
horizontal wavelength 4 times the boundary layer height as shown by Kuettner 
(1971). 
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This study has also demonstrated the feasibility of using the NASA 
airborne Doppler lidar system for observational studies of the structure of 
the clear air boundary layer. 



PART III: DOPPLER LIDAR OBSERVATIONS IN 

THUNDERSTORM ENVIRONMENTS 
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ABSTRACT 


On 30 June 1981, the wind fields around a variety of convective clouds, 
ranging from large thunderstorm complexes to isolated cumulus congestus, were 
observed in some detail using an airborne Doppler lidar operated by the 
National Aeronautics and Space Administration (NASA). Despite the presence of 
errors caused by problems in sampling some of the aircraft attitude and motion 
parameters used in the data analysis, reasonably clear pictures of the main 
features of the flow fields near several clouds have been obtained. 

Features of special interest seen in the lidar analyses are waves and 
vortices at the leading edge of a gust front marked by an arcus cloud forma- 
tion, and secondary surges of outflow behind another gust front which was 
mostly free of clouds and rain at flight level. Observations from one of the 
NSSL mesonet stations overflown by the lidar aircraft in the vicinity of the 
arcus cloud serve to confirm the approximate location and character of the 
gust front at the surface depicted by the lidar analyses in that area. The 
storm system, which contained the gust front vortices, moved eastward, inten- 
sified and within one hour produced a damaging gust front tornado in Norman, 
OK. 

Also seen by the lidar were clear air flows near the base and sides of an 
isolated cumulus congestus circumnavigated at two levels by the NASA air- 
craft. Prominent aspects of the observed flows near the isolated cumulus 
included 1 km-scale bands of enhanced convergence embedded in an overall pat- 
tern of convergence into the cloud region. These patterns were seen at both 
levels examined, 1100 m AGL (roughly at cloud base) and 2300 m AGL but must be 
regarded as suspect due to the presence of uncorrectable measurement errors 
associated with the navigation and aircraft motion sampling difficulties. 



3.1. INTRODUCTION 

The "Severe Storms Experiment" was designed specifically to study the 
clear air flows near convective clouds, especially those associated with 
severe thunderstorms, that might occur within range of one or both of the NSSL 
Doppler radars. Part III of this report concerns itself with the results of 
this experiment. 

With the advent of Doppler radar systems, it first became possible to 
observe the internal wind flow patterns of convective storms in some detail 
(Lhermitte, 1970, Ray, et al., 1975). However, in sampling air velocities at 
the edges of convective storms, the extremely small reflectivities in the 
"clear" air allow the much larger reflectivities of the adjacent storm, seen 
through the sidelobes, to contaminate Doppler measurements in the radar's 
resolution volume (Doviak and Zrnic', 1984), yielding radial velocity esti- 
mates which are biased toward storm-internal target motions, furthermore, 
radial velocity estimates at low radar elevation angle and close range may 
also be biased (toward zero velocity) by ground clutter. At longer ranges, 
the si delobe contamination and large radar resolution volume proportions can 
make it difficult to obtain adequately detailed pictures of the wind field 
near the edges of storms. In addition, storm features below cloud base are 
often shielded from the radar's view due to the earth's curvature when the 
storms are located beyond ranges of about 60 km (Doviak, et al., 1979, Doviak 
and Zrnic' , 1984). 

Doppler lidars not only feature fine resolution with very compact side- 
lobes, but also can be designed such that range and Doppler velocity ambigui- 
ties will not be encountered in normal use. However, the infrared frequencies 
employed in many systems suffer significant attenuation during passage through 
moist air due to absorption by water vapor and carbon dioxide. Thus, Doppler 
lidars are well-suited to the task of detecting clear air winds outside of 
storm clouds, but can only probe short distances into the clouds, while 
Doppler radars perform best in the hydrometeor -rich environment within the 
clouds. 

Although the wind field maps furnished by Doppler radars and radar arrays 
have been extremely helpful in revealing some aspects of the dynamics of 
severe storms, the limitations of such radars have prevented researchers from 
obtaining the definitive answers to many questions having to do with storm 
dynamics outside precipitation regions at cloud edge and cloud base. 
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For instance, better knowledge of the structure and interactions of 
thunderstorm inflow and outflow below cloud base and their implications with 
respect to production of damaging straight line winds, or of the amplification 
of vorticity during tornadogenesis, would be extremely useful to severe storm 
researchers and forecasters. While most severe storm investigators agree that 
the tilting and stretching of horizontal vorticity associated with vertical 
shear by differential rising motions and accelerations are of great importance 
in tornadic storms (Davies-Jones, 1982), there is still uncertainty regarding 
the impact on storm evolution of the many possible patterns of ambient verti- 
cal shear, buoyancy, and other parameters. Klemp and Wilhelmson (1978), using 
numerical simulation, have found storm evolution to be especially sensitive to 
the intensity and distribution of ambient vertical shear in the lower tropo- 
sphere. Storms simulated by Weisman and Klemp (1982) exhibited behavior which 
appeared to be governed by the environmental value of a bulk Richardson Number 
relating buoyancy to vertical shear. Although direct verification of such 
simulations is accomplished using detailed observations, the results and con- 
clusions drawn from them must be regarded as tentative. 

With respect to more general questions regarding cumulus convection dyna- 
mics, a better understanding of the distribution of entrainment over the sur- 
faces of convective clouds is eagerly sought by cumulus cloud modelers who 
continue to struggle with overly specific formulations which lack the desired 
general physical validity (Warner, 1970; Simpson et al 1982). 

Although there has been little doubt of the existence of entrainment in 
cumulus clouds since the pioneering work of Stommel (1947), there has been 
considerable controversy over the relative importance of lateral versus cloud- 
top entrainment mechanisms. A number of authors (Morton, 1957, Turner, 1962, 
Squires and Turner, 1962) have presented theoretical models of cumulus 
updrafts which require lateral entrainment during the cloud's growing phase, 
and some workers (e.g., McCarthy, 1974) have found confirmation of the 
inverse-cloud-diameter dependence of the entrainment rate predicted by most of 
these models. However, others (Sloss, 1967) have been unable to verify this 
inverse relationship between cloud size and entrainment rate. 

On the other hand. Squires (1958) and Fraser (1968) have proposed that a 
substantial proportion of total entrainment takes place through mixing at 
cloud top. Observational evidence in favor of their ideas has been obtained 
by Paluch (1979) and others. It seems likely that both lateral and cloud top 
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entrainment mechanisms work simultaneously on most convective clouds, though 
.he relative importance of each as a function of cloud size, shape and stage 
jf evolution remains an incompletely answered question. 

Mounted aboard a Convair CV-990 jet aircraft, NASA's airborne Doppler 
lidar system (ADLS) not only offers good mobility but also fore-aft beam 
scanning during data collection so that differing views of wind features can 
be obtained. Retrieval of the full horizontal wind vector at a given point 
involves solution of two simultaneous equations, one from each of two 
different but well-defined 1 ine-of -sight views of the wind field at that 
point. As will be seen later, precise knowledge of aircraft attitude and 
motion is necessary in order to establish accurate estimates of the wind 
vectors. 

It is the principal purpose of this study to see whether data collected 
by an airborne Doppler lidar system can provide useful insights into the wind 
flows near the edges and bases of a variety of convective clouds. 

3.2. METEOROLOGICAL SETTING FOR THE EXPERIMENT 

On 30 June, the second day of availability of the lidar, the morning 
meteorological situation in Central Oklahoma hinted at the possibility of 
strong thunderstorm activity late that day. A weak cold front had entered the 
northern part of the state the previous evening and had become quasi - 
stationary, lying roughly along an east-west line which meandered along the 
Kansas-Oklahoma border. McCaul (1985) gives a detailed description of the 
synoptic setting but we present here a few salient features. 

The 1200 GMT (0600 CST) sounding from Oklahoma City (0KC) indicated 
unusually deep moisture with little or no capping inversion above the boundary 
layer. Winds were light and veered from south at the surface to northwest 
above the 50 KPA level, indicating a deep layer of weak warm advection. The 
sounding had not changed significantly by the 1855 GMT (1255 CST) launch time 
of a special rawinsonde from Tinker Air Force Base (TIK) where the lidar 
aircraft was stationed. Fig. 3.1 depicts the temperature, moisture and wind 
profile of the special TIK sounding. Of special interest in Fig. 3.1 is the 
presence of convective instability in lower layers, the lifted index of this 
atmosphere was about -3, and the convective available potential energy (CAPE) 
per unit mass for boundary layer parcels was about 1830 m^s"^. Precipitable 
water amounted to 5.08 cm. The convective temperature was estimated to be 
29°C, which was just being reached at sounding time. 
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Figure 3.1 (a) Skew T-log p diagram; Tinker Air Force Base, 1855 GMT 

(1255 CST), 30 June 1981. (b) Hodograph of winds. Each circle represents 

5 kt of windspeed and the numbers next to the curve give the altitude in 
thousands of feet above mean sea level. 

The most striking feature of the wind profile is that at no point does 
the wind exceed 12 m s - * (25 kt; in magnitude. There is a tendency for the 
shear vector, which always points parallel to the tangent to the hodograph 
itself, to rotate to the right with increasing altitude, at least in the lower 
half of the troposphere. Because of this shear vector rotation, tilting by 
updrafts of the horizontal vortex tubes, which lie everywhere perpendicular to 
the shear vectors, could in principle lead to production of the cyclonically 
rotating updrafts often seen in “supercelT' storms (Browning, 1964; Barnes, 
1968). In this case, however, the winds and shears were so weak that the pro- 
duction of rotating updrafts was likely to be inhibited by the inability of 
the shear to relieve the updrafts of sufficient amounts of condensed water 
load. Excessive water loading in storm updrafts is believed to be a major 
factor constraining the buoyancy and hence, presumably, also the persistence 
of storm updrafts (Lilly, 1979). 

Around the time of the special TIK sounding, surface temperatures in many 
places in the state had risen above 30°C and deep convection was beginning 
west and southwest of Oklahoma City. At 2000 GMT (1400 CST) surface tempera- 
tures of 33°C were in southern Oklahoma at Ardmore, while a heavy thunder- 
shower was in progress at Ft. Sill (FSI) about 100 km to the southwest of 
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Norman, OK. Shortly before 2000 GMT it was apparent that the convection was 
becoming intense; thus, at about 2000 GMT the lidar aircraft departed TIK to 
collect data for the Severe Storms Experiment. 

For the reasons mentioned earlier, the storms did not develop into 
classic supercells, but rather exhibited the characteristics of "multicell" 
storms (Byers and Braham, 1949). However, according to a report in STORM DATA 
for June 1981, one of the storms produced damaging winds and a brief gust- 
front tornado in Norman, OK, around 2120 GMT (1520 CST), and thus satisfied 
the National Weather Service (NWS) criteria for being "severe". According to 
OKC radar summaries, storm tops were observed to reach 60 kft (18.3 km), which 
seems consistent with the thermodynamic profile of Fig. 3.1 in view of the 
tropopause height of 47 kft (14.3 km). 

3.3. DATA COLLECTION 

The CV990 collected 15 runs of lidar data during its 30 June 1981 flight, 
which lasted from 1959 GMT (1359 CST) to 2244 GMT (1644 CST). Only runs 2 
through 12 were used for this study because the other runs were taken before 
or between passes near convective activity. A map showing the flight track of 
the CV990 during the period when data runs 2 through 12 were being collected 
is given in Fig. 3.2. 

Run 2 commenced at 2018 GMT (1418 CST) and run 12 ended at 2100 GMT (1500 
CST). Thus, all the data used in this study were taken within a 42-minute 
period near the beginning of the flight. 

The CV990 flew southwest toward the vicinity of Chickasha, OK, (CHK) 
immediately after takeoff, taking a preliminary run of data during ascent. By 
the time the aircraft started to approach Chickasha, several very strong thun- 
derstorm cells had begun to develop in a loosely organized cluster just west 
and northwest of that City. Although the precipitation shafts were still 
distinct at this time, about 2015 GMT (1415 CST) outflow from the several 
cells was beginning to merge and move east-northeast behind a well-defined, 
nearly linear arcus cloud formation. 

At about 2018 GMT (1418 CST) the CV990 veered toward the northwest to 
make a data-gathering run (run 2) along a flight track roughly parallel to 
this advancing quasi -linear outflow feature. The aircraft was flying at 
approximately 720 m above ground level (AGL) gathering data during run 2 which 
lasted more than 7 minutes. During this period the lidar scanned the arcus 
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Figure 3.2 Flight track, NASA CV-990, 30 June 1981. Times (GMT) are shown 
every 20 min with aircraft altitude (AGL) printed just below the times. Radar 
reflectivities (dBZ) at 2037 GMT are shown in dashed contours. 

cloud formation or the edges of heavy precipitation continuously, with laser 
frame duty cycles ranging anywhere from 1.25 s to 1.5 s. 

Figure 3.3, one of a series of photographs (termed DF photos) taken by 
Dr. Dan Fitzjarrald on board the CV990 shows two things: (a) an isolated 

heavy thunderstorm rain core near Chickasha and (b) the arcus formation at the 
southern edge of another storm core to the northwest as viewed from the CV990 
at about 2019 GMT (1419 CST). 

Approximately 30 s after the time of Fig. 3.3, the CV990 came abreast of 
the arcus cloud formation surrounding the advancing edges of the northwestern 
storm. Judging from mission cloud photographs taken regularly from the 
aircraft at 1 s intervals and covering the lidar's scanning area, the base of 
the arcus cloud occasionally reached as low as an estimated 400 m above the 
ground. Although the mission cloud photos were useful in documenting the 
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Figure 3.3 Photo, looking southwest, ca. 201900 GMT (141900 CST), during 
run 2. Note rain core at center and arcus at right. Aircraft was moving 
toward the right. 

cloud forms scanned by the lidar, they were unfortunately damaged somewhat 
during processing and thus were unsuitable for reproduction here. Neverthe- 
less a DF photo taken just before 2020 GMT (1420 CST) and shown in Fig. 3.4 
reveals the structure of the arcus formation as broadly lobed at a horizontal 
wavelength of about 300 m, with irregularities in the altitude of the cloud 
base. Other mission cloud photos show that the arcus cloud did not display an 
especially straight front edge, but rather meandered somewhat at scales of 3-5 
km along its length. 

At about 2020 GMT (1420 CST) the CV990 flew just southwest of Tuttle, OK, 
and within 1 km of one of NSSL's mesonet stations, Tuttle-South (TTS), which 
was operating at that time in support of the lidar experiments. Data from 
that station indicate the passage of a moderately strong gust front at about 
2023 GMT (1423 CST), followed by heavy rain some 14 minutes later. These data 
are presented and discussed in Section 3.5.1. 

By 2025 GMT (1425 CST) the CV990 had reached the northernmost edge of the 
developing convective complex, and thus commenced a 90° turn to the left to 
scan the northern edge of the storm. Run 3 of lidar data began at about 
2026 GMT (1426 CST) and continued until just after 2030 GMT (1430 CST). 
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Figure 3.4 Photo, looking southwest, ca. 201930 GMT (141930. CST), during 
run 2. Note continuous, lobed arcus cloud in foreground, rain beyond it. 
Aircraft was moving toward the right. 


Typical frame cycles in run 3 ranged from 1.22 to 1.50 s. At about 2027 GMT 
(1427 CST) the aircraft passed near the front edge of a bulging gust front 
which, judging from the orientation of the edge of the outflow as seen in the 
photographs and in the lidar analyses, was pushing toward the north and 
northeast. The aircraft, flying at an altitude of 550 m AGL, obtained a 
number of lidar velocity measurements at the edge of and just behind this gust 
front. In contrast to the well-defined arcus cloud previously seen on the 
east side of the storms and studied in run 2, this gust front was marked only 
by scud having a ragged, shredded appearance, often arranged in irregularly- 
shaped vertical filaments (Fig. 3.5). 

The lidar aircraft began to pass out of the vicinity of the active 
convection just after 2030 GMT (1430 CST) as it continued toward the 
southwest. Consequently, at 2032 GMT (1432 CST) a left turn toward the 
southeast was made in order to fly along the back (or southwest) edge of the 
storm area. 

Run 4 of lidar data was taken starting at 2036 GMT (1436 CST) and ending 
at about 2044 GMT (1444 CST) at altitudes ranging near 800 m and with frame 
cycles ranging from 1.12 to 1.40 s. Little or no outflow cloudiness was seen 
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Figure 3.5 Photo, looking south, ca. 202700 GMT (142700 CST), during run 3. 
Note presence of detached scud instead of arcus clouds. Aircraft was moving 
toward the right. 

on the southwest side of the storms and only light rain was falling from an 
amorphous cloud base in that part of the complex. However, at 2039 GMT 
(1439 CST) the aircraft pulled away from this area of lingering light rain 
into an area of distinctly warmer air where cloud bases were noticeably 
better-defined, with bases at an estimated 1. 5-2.0 km AGL. Sensors aboard the 
aircraft confirmed not only a 1°C ambient temperature increase, but also a 
significant 4°C rise in dew point. The location of this apparent boundary 
between rain-cooled and undisturbed air was about 12 km southwest of 
Chickasha. At this point, the CV990 turned briefly toward the east in order 
to obtain data closer to the retreating storms, then turned at 2040 GMT 
(1440 CST) back to the southeast to pass by the rear of another rain cell. 

At about 2042 GMT (1442 CST) it was noticed that an isolated cumulus 
congestus was beginning to develop about 40 km to the southeast. 

Consequently, it was decided that several circumnavigations should be made of 
this new convective growth, which was centered about 10 km west of 
Elmore City, OK (ELC). 
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The CV990 approached the cumulus congestus from the north and began a 
square flight pattern at an altitude of about 1100 M AGL around the base of 
the cloud at 2046 GMT (1446 CST). Run 5 consisted of the data taken along the 
west side of the cloud as the aircraft headed south. Runs 6, 7 and 8 were 
taken on the other sides of the cloud after a succession of 90° left turns. 
Each run at the 1100 m level lasted an average of about 60 s. 

Figure 3.6, a DF photo taken during run 7 while looking west toward the 
cumulus, reveals that the cloud tilted toward the south with height, 
suggesting the updraft was being subjected to a vertical shear which was 
northerly, at least through mid-height of the cloud. 

At the conclusion of run 8 the CV990 was again at the northwest corner of 
the cloud, and the time was about 2052 GMT (1452 CST). The aircraft climbed 
as it headed south along the west side of the cloud in preparation for another 
circumnavigation at a higher level. An altitude of about 2300 m AGL had been 
attained when run 9 began at 2054 GMT (1454 CST) following a 90° left turn 
southwest of the cloud. Run 9 was taken as the aircraft headed east along the 
south side of the cloud, which could now be said to consist of a series of 



Figure 3.6 Photo, looking west, ca. 205030 GMT (145030 CST), during run 7. 
Note towering cumulus being undercut and sheared off to south. Aircraft was 
moving toward the right. 
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contiguous towers oriented along a northeast -southwest line. Runs 10 through 
12 were taken along the other edges of a square flight track around the cloud, 
each one lasting about 60 s average, as did runs 5-8. At the conclusion of 
run 12 the time was about 2100 GMT (1500 CST). 

The towering cumulus did not develop into a thunderstorm, instead the 
uppermost parts of the turrets were observed to shear off toward the east and 
southeast before reaching the 10 km level, and no significant anvil cirrus 
formed. 

3.4. DATA ANALYSIS 

NASA furnished a data tape for the 30 June experiment which was to 
contain the following information: 

A. Header information regarding dates and times of each run; 

B. Position, attitude, groundspeed and heading information from the CV990 
INS for each frame of each run, 

C. Aircraft altitude measured by down-looking radar during each frame, 

D. INS-derived wind vector at the aircraft for each frame of each run, 

E. Airspeed, ambient temperatures and dew points measured by the gust 
probe; 

F. Number of transmitted pulses per laser frame; 

G. COHO and VCO offset frequencies for each frame, 

H. Elevation angle of line-pf -sight and positions of scanner wedges for 
each frame, 

I. Times and line-of -sight azimuths for each lidar frame; 

J. Raw Doppler velocities, echo intensities and spectral widths for each 
resolution volume along each frame, 

K. Edited and smoothed Doppler velocities and range-normalized 
intensities for each resolution volume of each frame; 

L. Smoothing coefficients used at each resolution volume of each frame. 
These were computed according to local data quality based on signal - 
to-noise considerations (Lee, 1981); 

M. Vorticity, divergence, stretching and shearing deformation of the flow 
at each resolution volume of each frame; 

Upon inspection of the data from a number of runs, it was found that many 
of the data fields were either missing or erroneous, making it necessary to 
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ascertain whether or not the analysis of the data could be accomplished. The 
data that were missing, such as COHO and VCO frequencies, number of pulses per 
frame, pressure altitude, track angle and line of sight elevation, were either 
not critically important for the velocity analysis or were deducible from 
other given quantities. However, as discussed in Appendix A, the VCO fre- 
quency is an important parameter to record especially if the frequency is not 
tracking the Doppler shift associated with the measured air speed. The erro- 
neous data, including such other important quantities as the corrected and 
Smoothed lidar velocities and range-normalized lidar intensities, had to be 
recomputed anew from the raw data. Because the smoothing procedure (Lee, 

1981) u$ed by NASA was rather elaborate and would have been time-consuming to 
implement, a somewhat simpler approach to smoothing was ultimately taken, 
obviating the need to reconstruct the error-contaminated smoothing coefficient 
arrays furnished on tape. 

Thus, it became necessary to perform a number of "check 1 ' computations of 
certain of the raw data fields in order to verify those given and obtain those 
missing. In the process of doing these computations, other errors in the data 
became apparent. In the following sections these errors and omissions and the 
techniques employed for correcting them are briefly mentioned. A complete 
discussion of the verification and correction techniques relevant to CV990 
navigation data is given in Appendix C, while a similar discussion regarding 
verification and correction of the lidar moment data is given in Appendix D. 
3.4.1 Navigation Data Review and Editing 

Inertial navigation systems are subject to a variety of errors (Britting, 
1971), one of which involves oscillations at the 84.4-minute Schuler period. 
INS velocity errors having amplitude of about 4 m s"* and period near the 
Schuler period were discovered in the CV990 data and are described in 
Part I. Unfortunately, the approximately 30-minute long record of nadir 
photos which was available for the Severe Storms flight was too short to allow 
determination of whether or not Schuler periods were present in the INS errors 
which occurred in that flight. 

Based on the amount of apparent variability in the INS errors and the 
timing of changes in direction of the errors seen during the flight, the most 
likely cause of the INS errors in the present data is faulty measurement of 
accelerations during sharp turns. A review of Fig. C.2 in Appendix C.3 



reveals that aircraft position errors were about 1.6 km in magnitude during 
run 4; this is somewhat surprising since run 4 took place only about 40 min 
into the flight, and position errors are generally expected to be small during 
the early stages of a flight. Position errors of about 1.6 km were also noted 
as early as 20 min into the flight, which suggest that there was an inital 
position bias of 1.6 km ( = 1 nautical mile); this initial bias might have 
been caused by a discrepancy of 1 arc min of latitude in the starting coor- 
dinates entered into the INS at take-off. 

In reviewing the aircraft navigation data, the recorded values of air- 
craft latitude and longitude were found to contain errors which were the 
result of sampling the INS position output too infrequently. Thus, it was 
necessary to establish aircraft positions by forward-time integration of air- 
craft track and velocity data. Track angles had to be computed from drift 
angle and heading data, because track angle data were not furnished on the 
data tape. 

Obvious errors were noted in certain of the data arrays furnished on the 
tape, apparently caused by use of incorrect conversion factors during the 
writing of the tape (Appendix C.l). These errors had to be diagnosed and cor- 
rected. 

The times of the various lidar data frames were also subject to round-off 
errors, and had to be reconstructed by an interpolation process so that the 
actual distribution of lidar data points in space would be faithfully repre- 
sented in the analyses. These corrections are described in Appendix C.2. 

The corrected INS-derived aircraft positions were further tested against 
nadir-photography-derived positions for those portions of the flight where 
such photo data existed (Appendix C.3). Taking the smoothed sequences of 
photo-derived positions as true indicators of aircraft position, a variable 
southward position error of about 1.7 km average magnitude was noted (see Fig. 
C.2 of Appendix C.3); this position error could be described in terms of a 
constant bias in position plus superimposed time-dependent fluctuations in the 
position error. While a constant position bias would not corrupt wind field 
estimates, the time-dependent fluctuations suggested time-dependent errors in 
the ground velocity estimates which, in turn, can ruin wind estimates. Based 
on the behavior of the smoothed time-differentiated INS position errors, the 
ground velocity errors were found to be as large as 10 m s~* or more during 
turns, and perhaps 2-4 m s - * during the several km of straight flight fol- 
lowing turns. 



Because the 10 s interval between nadir photos was too large to allow 
high-resolution study of the INS errors, and since other time-varying errors 
of nearly equal magnitude (about 2 m s"-*-) were occurring along with the INS 
errors, only rough attempts at diagnosis and removal of these errors could be 
made. Such attempts were useful not so much for the marginal improvements in 
data quality which they produced, but rather to demonstrate quantitatively the 
magnitude of the harmful effects of the time-varying INS errors on the 
retrieved wind vector fields. 

Although some success was obtained in defining the magnitude of the INS 
errors via the techniques outlined in Appendix C, the discussions there make 
it clear that it is not yet feasible to make precise corrections. Hence, in 
all figures of wind analyses to be presented in Section 3.5, the INS errors 
have not been removed except when specially noted. The principal benefit to 
accrue from the study of the INS errors described in Appendix C has been to 
give some rough idea of the magnitude of the errors, their time variability 
and of the level of confidence which can be assigned to the details of the 
analyzed wind fields. 

3.4.2. Doppler Spectral Moment Data Review and Editing 

Review and editing were needed not only for navigation data but also for 
Doppler moment data. After insuring self-consistency in all aircraft position 
and attitude data, the next step in the analysis involved examining fields of 
Doppler moment data and removing any errors and inconsistencies affecting 
these data. Modifications to data were needed in some cases to prepare it for 
subsequent processing and analysis. 

The reflectivity estimates furnished by NASA showed evidence of being 
erroneous, so that these data had to be recomputed starting from the raw 
intensity data (Appendix D.l). The reflectivity estimation procedure employed 
in this study included correction for the inverse range-squared divergence 
experienced by backscattered photons as well as round-trip correction for 
absorption losses. The resulting normalized reflectivities were useful in 
identifying the positions of brightly reflective features such as cloud edges. 

The Doppler velocity data supplied by NASA were also contaminated by 
errors and a full discussion of these and methods used to correct them are 
presented in Appendix D.2. It appears that significant errors in radial 
velocity estimates were caused by delays in reporting time varying INS 
measured parameters. 
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An optimal correction for these delays would have required an exhaustive 
investigation into the system electronics and data transfer queueing, a task 
vhich was beyond the scope of this study. However, the delays which most sig- 
Tificantly affected the quality of the Doppler data were those involving the 
recording of drift angle (Appendix D.3) and ground velocity (Appendix D.4). 

The drift angle was observed to be a small but rapidly fluctuating quantity to 
which the radial velocities were quite sensitive, while the groundspeed was a 
large and relatively slowly varying quantity to which the radial velocities 
were somewhat less sensitive. 

For rapidly varying quantities such as drift angle, the effect of the 
delays was to produce seemingly random errors (i.e. frame bias; Appendix D.3) 
in the measured range-averaged radial velocities, such errors could be effec- 
tively eliminated by simple filtering techniques. Although filtering was nor- 
mally used to reduce the high-frequency errors caused by drift angle delays, 
tests of a delay correction algorithm were also routinely made; these tests 
met with limited success in some of the data runs. 

For the more slowly varying quantities, cross-correlation techniques 
could theoretically be used to determine delays of various parameters with 
respect to the other related parameters (McCaul, 1985). However, this tech- 
nique could not be applied to the present data in a straightforward manner due 
to simultaneous presence of several interacting errors, such as drift angle 
error and ground velocity error, each of which varied without well-defined 
periods. Even after attempts to remove ground velocity errors from the data, 
the cross-correlation technique failed because of the presence of lingering 
remnants of ground velocity error; techniques for elimination of the slowly 
varying ground velocity errors are not yet sophisticated enough to permit 
further analysis of those errors caused by delay of other slowly varying 
parameters. 

Spectral width values furnished by the system were found to be somewhat 
more subject to noise than the Doppler velocities and intensities. Neverthe- 
less, in many instances local maxima of spectral width were found to correlate 
well with positions of wind field boundaries or gradients. 

3.4.3. Velocity Errors Due to Beam Pointing Errors 

Velocity errors could be introduced through malfunctions of the wedge 
scanner which directs the beam onto a supposedly horizontal trajectory into 
the atmosphere. 



Based on Lee's findings of beam pointing biases, some uncertainty ini- 
tially surrounded the performance of the beam scanner. This uncertainty, 
caused by the failure of the scanner to compensate for aircraft roll, prompted 
a decision to delete all lidar data whenever the aircraft roll angle exceeded 
1° in magnitude. However, based on results to be discussed in Appendix E, it 
appears that these deletions may have been unnecessary. However, the amount 
of data deleted by this roll angle criterion was quite small. 

Doppler shifts recorded along radial s which lie out of the horizontal 
plane can introduce considerable error in presence of vertical shear, as well 
as have contamination from any vertical velocity that might be present. 

Because data are collected only on a surface (hopefully a horizontal one), 
insufficient data are available for proper separation of the horizontal and 
vertical components of air motion, and thus data taken "out-of -plane" cannot 
strictly be corrected to be compatible with data from the horizontal plane. 

Further problems related to beam elevation uncertainty arise because the 
lidar beam scanner is held fixed during the 0.4 s duration of each frame, 
allowing the sequence of pulses to sweep through a variety of elevation angles 
whenever the aircraft rolls in response to turbulence. 

Velocity errors can also result from beam pointing errors because they 
can cause incorrect components of aircraft ground velocity to be included in 
the radial velocity data. The most reliable way of assessing beam pointing is 
through study of the variation of range of ground strikes versus the known 
variation in form of some topographic feature by which the aircraft passed 
during data collection. Although no beam strikes on mountainsides occurred 
during the Oklahoma flights, several beam ground strikes were observed during 
one aircraft left turn; the variation of range of these ground strikes versus 
the range expected based on aircraft altitude and roll angle provided an 
opportunity to make a limited assessment of the accuracy of the beam poin- 
ting. The results of these comparisons, discussed more fully in Appendix E, 
indicate that beam pointing errors were limited to less than 1-2 deg in the 
vertical. Using a technique outlined by McCaul (1985, Appendix F) involving 
comparison of derived wind vectors from regions where data from three or more 
lidar lines-of -sight were collected, it was possible to conclude that there 
were no apparent significant errors in horizontal beam pointing. 

A conservative estimate of the windspeed error in these data caused by 
beam elevation errors is 0.5 m s - * except above the gust fronts of runs 2 and 



and along the cloud turrets of runs 9-12, where local errors could reach 2-3 
S'” 1 . Clearly, the problems of laser misalignment and failure of the scanner 
o provide the proper beam elevation compensation during aircraft rolls should 
e corrected before the next series of experiments. In addition, it would be 
lesirable if scanner positioning could be updated at least occasionally during 
;ach lidar frame so as to ensure horizontal beam trajectories for all pulses. 

5.4.4. Preparation of Meteorological Data Fields 

In order to synthesize meteorological data (i.e., reflectivity, wind, 
urbulence) fields from the edited moment data aligned along the laser beam 
tnd to present the results in a suitable form for display, it was convenient 
.o interpolate moment data to a Cartesian grid whose origin was placed at the 
ocation of NSSL's Cimarron Doppler radar because intercomparison of Doppler 
ind lidar wind fields were made (McCaul, 1985). Before the grid point values 
:ould be assigned, advection compensation had to be made. 

i) Data Advection 

The nonsimultaneity of the data collected by the lidar during each run, 
in some cases spanning about 500 s, required that attention be given to the 
problem of data advection. This was not so much a problem with respect to 
obtaining locally correct retrieval of wind or intensity patterns, which was 
governed by the maximum time delay (about 50 s) between fore and aft frame 
data used in analyzing fields at maximum ranges (about 10 km) from the air- 
craft, as it was a problem of obtaining a correct pattern over the larger 
scales. 

In order to avoid distortion of the data patterns found near the storm 
outflow fronts, it was decided that advection proportional to the storm's 
velocity vector would be uniformly applied to all data in each run. Although 
this decision undoubtedly fails to account for some of the complexity of the 
three-dimensionality of outflows observed in the data, it is believed that the 
resulting fields are less error-prone than they would have been had local 
advection been used. 

Tracking of the main storm cells on successive photographic images of 
NSSL's WSR-57 radar screen suggested a storm motion vector of 12.9 m s _1 from 
256°. Due to rapid development and evolution of cells in this multi cell 
complex as observed on the radar screen and in the field, this advection 



vector must not be considered an absolutely accurate measure of how the cells 
moved or how the windfield patterns propagated. 

Calculations of the mean wind from each data run used in the study of the 
main storm complex give values which differ substantially from the radar- 
derived motion vector. In fact, lidar run 4 showed an average wind vector of 
11.5 m s* 1 from 212°. This discrepancy may be taken as a rough indicator of 
the relative importance of propagation or differential motions in the wind 
fields studied here. In Appendix F, a model wind field, with a velocity 
(i.e., phase speed) of wind perturbations known but different than the mean 
wind in the layer, was sampled in exactly the same manner as the actual 
data. Wind vectors synthesized from radial velocity fields using the correct 
advection velocity and that equal to the mean wind in the layer show consider- 
able differences. 

ii) Data Interpolation and Wind Estimation 

Interpolation to a Cartesian grid was done using the circularly symmetric 
distance weighting function (Cressman 1959) and details are discussed by 
McCaul (1985). To give maximum resolution to the features of interest, 
slightly di fferent .radi i of influence were used for the interpolation velo- 
city, intensity and spectral width data. The radius of influence for interpo- 
lation of velocity data was chosen to be 600 m in order to ensure that data 
from three or four different frames would be interpolated, but that data from 
only two or three resolution volumes along the beam would be used. This is 
important because of the tendency of the radial velocity data to show more 
variability from frame to frame (due to the drift -angle-related errors 
described earlier) than from range bin to range bin. 

The radius of influence for the analysis of lidar echo intensity and 
spectral width was set at a smaller value, 400 m, in order to give maximum 
expression to the large gradients which occurred with these quantities in 
zones of strongest meteorological interest. 

Although the interpolation procedure for the scalar fields followed 
standard procedures (see McCaul 1985), several additional steps were needed in 
the process of interpolating data for wind estimates. Because wind is a 
vector quantity derived from the scalar radial velocity data, the inter- 
polation process had to be done four times: once for the forward frame velo- 

city data, once the forward line-of -sight data (i.e., angular direction of the 
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beam), once for the aft frame velocity data, and finally, once for the aft 
line-of -sight data within the neighborhood of each grid point. The full 
horizontal wind vector is synthesized only after the two differently directed 
sets of radial velocity data and two differently directed lines of sight data 
had been separately interpolated to the Cartesian grid. 

The geometric motivation behind the solution for the total vector is 
shown in Fig. 3.7. Solutions may be obtained in a variety of forms, the 
option used in this study retrieves the wind vector in terms of a magnitude 
"V" and an azimuth "AZ" according to 


AZ = ARCTAN 


f VRl*C0S(AZL0S2)-VR2*C0S(AZL0Sl) , 
IVR2*SIN(AZL0S1)-VR1*SIN(AZL0S2 ) j 


(3.1a) 


V = -VR2/C0S(AZ-AZL0S2) = -VR1/C0S(AZ-AZL0S1) (3.1b) 


where "AZL0S1" and "AZL0S2" are the fore and aft interpolated lines-of -sight 
and "VR1" and "VR2" are the fore and aft interpolated radial velocities valid 
at the grid point in question. In applying Eq. 3.1b, the solution for "V" is 
obtained using the form having the larger cosine value in the denominator. A 



Figure 3.7 Construction of wind 
vector from fore and aft radial 
velocities. Intersection of 
perpendiculars from the tips of the 
components defines the wind vector 
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sample plot of wind vectors thus obtained from the raw data is shown in Fig. 
3.8, along with the radial velocity data used in their calculation. 

3.5. METEOROLOGICAL INTERPRETATION 

In this section we present the results of the lidar data analysis in the 
vicinity of thunderstorm outflows and around an isolated cumulus congestus 
cloud. 


3.5.1. Interpretation of Thunderstorm Outflow Data 

In runs 2 and 3 of the lidar data, the horizontal structural features of 
several thunderstorm outflows were examined. Results from run 2 were obtained 
at an altitude of about 720 m AGL, while those from run 3 were obtained at 
about 510 m AGL. 

In both cases lidar range was variable, with a minimum of 3 km near 
dense, low clouds and a maximum of 8 km in other areas. Strong attenuation 
was noted when the lidar beam entered arcus or scud, as in run 2, but less 
attenuation was observed in run 3 when the outflow was free of low clouds. 

Run 2 was made on the northeast side of a complex of strong storm cells 
developing from Chickasha to El Reno, OK; the derived wind field is shown in 
Fig. 3.9 while the contoured fields of return intensity and spectral width are 
shown in Figs. 3.10 and 3.11. In Fig. 3.9, light southeast winds are evident 
in the area just east of the outflow front, and strong gusty southwesterly 
winds (sometimes 20 m s"*) are seen at and behind the gust front. Fig. 3.9 
shows a portion of the gust front which stretches from coordinates (X = -5, 

Y = -33) northwestward to coordinates (X = -28, Y = 0). Figs. 3.10 and 3.11 
show a "wall" of high intensity returns and large spectra 1 widths near the 
wind shift line seen in Fig. 3.9. As the evidence to be presented shortly 
will show, this "wall" of high intensity and spectral width was most likely 
due to the lidar signal impinging on the arcus cloud (Fig. 3.4) which occurred 
along much of the gust front in run 2. 

In Fig. 3.11, the large spectral widths associated with the arcus cloud 
appear on the southwest edge of the data swath and should not be confused with 
the line of large spectral widths lying along the aircraft track on the north- 
east edge of the data swath. This line is an artifact of width values which 
are reported in the first range bin of each frame, and has no meteorological 
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X-DISTANCE FROM CIM (km) 


Figure 3.8 Sample field of radial velocity data and synthesized wind 
vectors. Run 3. Shown are the successive aircraft positions and the 
interpolation circle. A scaled velocity vector is on the upper left of the 
figure. The radial velocity vectors are those aligned along the fore and aft 
directions indicated by the vector attached to each frame number. The 
synthesized wind vectors are those located on the cartesian grid with 250 m 
spacings. 
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Figure 3.9 Vector plot of lidar-deri ved winds of storm outflow observed 
during run 2. Note stronger winds along the front of the outflow at left edge 
of the data swath, and the irregular form of its leading edge. 
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Figure 3.10 Contour plot of the normalized reflectivity along the front of an 
advancing storm outflow. Note correspondence of high intensity areas with the 
leading edge of the outflow seen in Fig. 3.9. 
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Figure 3.11 Contour plot of spectral widths (turbulence) in the vicinity of a 
thunderstorm outflow. Note correspondence of high spectral widths with the 
leading edge of the outflow. 
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significance. Such spurious patterns in spectral width may be seen in the 
first range bin in every data run. These large widths, as well as those seen 
at cloud boundaries, might have been caused by the large echo powers (from 
targets close to the lidar or from highly reflective cloud particles) that 
exceed the dynamic range of the receiver (= 54 dB for the 10 bit analog-to- 
digital converter) causing the echo signal peaks to be clipped. When echo 
signals are clipped, odd harmonics are generated that have alternately posi- 
tive and negative frequencies about the mean Doppler shift of the echo and 
these harmonics can cause the spectral widths to be biased toward large values 
(Doviak and Zrnic', 1984, p. 158). 

Particularly noteworthy in Fig. 3.9 are the vortex patterns which appear 
at irregular 3-7 km intervals along the gust front. The best examples are 
indicated by circled "Vee's". These patterns could be caused either by 
(a) the tilting of vertical -shear-induced roll vortices in the near-storm 
boundary layer by the lifting action of the cold outflow at the edge of the 
gust front, or (b) to the convergence of the vorticity associated with the 
strong horizontal shears found along the gust front. 

Support for the former theory seems to be especially strong in the vicin- 
ity of the vortex located at (X = -9, Y = -24), where a surge of southeasterly 
wind is observed to extend from the aircraft's position back to the southwest 
into the vortex, producing a strong impression of a horizontal vortex tube 
seen in horizontal cross-section as it is tilted upward by the differential 
rising motion just ahead of the storm outflow. Analysis of this particular 
vortex pattern requires caution because spurious patterns can be produced by 
erroneous velocity data. Furthermore, no definite cloud pattern featuring a 
vortex in the arcus cloud could be seen in photographs (e.g., Fig. 3.4) taken 
at the time the lidar data were collected. The other vortices appear to be 
associated with notches or concavities seen both in the photographs (e.g.. 

Fig. 3.4) and in the reflectivity patterns and thus appear to be genuine 
despite the relatively weak horizontal winds in the warm air ahead of the gust 
front. 

Although lidar data were not collected at more than one level along this 
outflow (thus precluding a detailed three-dimensional diagnosis of the verti- 
cal shear distribution in the near-storm boundary layer), surface data from a 
NSSL mesonet station at Tuttle South (TTS) nearly under the aircraft flight 
track were available for comparison with the lidar-deri ved wind. Because the 



aircraft apparently was closest to TTS (within 1 km) at 201955 GMT (141955 
CST), the vertical shear at one point (TTS) could be estimated quite accu- 
rately by finding the difference between wind vectors obtained by the lidar 
and the surface station and dividing by aircraft altitude. Synthesis of the 
lidar data from bin 2 of frames 81 and 82 of the lidar data gave a wind of 

4.3 m s -1 from 135° at the TTS location at 201955 GMT, while the surface sta- 
tion gave a 60-s average wind of 1.9 m s"^ - from 134 J at 202000 GMT. Using the 
aircraft altitude estimate of 722 m from the CV99Q's down-looking radar and 
ignoring the tiny difference in wind azimuths, we obtain a vertical shear of 

3.3 x 10-V 1 . 

Thus, even though the winds were very light, there was vertical shear 
directed to the northwest in the immediate vicinity of the gust front that 
could have been tilted onto the vertical to produce the vorticity patterns 
seen in the run 2 data. To see whether the actual measured vorticities could 
be produced directly by the tilting mechanism, we examine the inviscid 
vorticity equation in height coordinates: 

3?/3t = • V U+f) - U+f) 6-k •M3^ H /3Z-k»VctXvP (3.2) 

where: ? = vertical vorticity 

f * planetary vorticity 
6 = horizontal divergence 
^ = total three-dimensional wind vector 
^ = horizontal wind vector 
W = vertical wind component 
a = specific volume 
P = pressure 

k = unit vector along Z-coordinate 

If we neglect the small terms involving advection of planetary vorticity and 
horizontal pressure-density solenoidal effects, Eq. 3.2 may be rewritten in 
terms of the rate of change of parcel vertical vorticity following the motion 
as follows: 


= - ( £+f ) 5 -k*$WX3if/3Z 


(3.3) 
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where the first term on the right represents convergence of pre-existing and 
planetary vorticity and the second represents tilting of pre-existing vertical 
shear. Assuming that vertical vorticity is initially negligible, and using 
observed convergence rates of 10 -2 s _1 and vertical shears of 10“ 2 S -1 along the 
gust front and an estimated gradient of rising motion of 10 m s -1 across 10 
km, or 10"3 S "1, we find a vorticity production rate of about 10"®s' 2 . Now 
assuming vorticity, vertical shear and gradients of rising motion are roughly 
independent quantities, it would take about lCr s to produce the observed vor- 
ticities of 10 _2 s"^, this time, which is approximately 3 hours, is too long to 
account for the observed vorticity development, which probably occurred in 
less than an hour. 

However, we note from Eq. 3.3 that once vorticity reaches about 10"^s"^, 
it no longer can be neglected as a factor in the total vorticity production 
rate. Considering tilting only, a vorticity of 10" 4 s"* could be generated in 
only 100 s (about 2 min). Once such a vorticity exists, it will amplify expo- 
nentially in the presence of the strong convergence observed. Amplification 
by a factor of 100 gives the required observed vorticity; this can be achieved 
in an additional time "T" given by the natural logarithm of 100 divided by the 
observed convergence of 10” 2 s - ^, or in only a time T * 460 s. Thus, we con- 
clude that in the vicinity of the gust front observed in this storm, total 
vorti cities of about 1CT 2 s'* could be produced in about 560 s, or less than 10 
min, due to tilting of initially horizontal vorticity and subsequent 
stretching by the strong frontal convergence field. 

In the preceding argument we have based our calculations on pre-existing 
vertical shears estimated from surface data and lidar data at an altitude 
above 700 m AGL. Simpson (1982) has argued that the evolution of some meso- 
vortices may be the result of tilting of the much larger vertical shears found 
very near the surface. Assuming such mechanisms to be acting on our data, it 
is quite likely that tilting alone could explain the bulk of the vorticity 
observed. 

If these vortices were indeed caused by the tilting of horizontal boun- 
dary layer vortex tubes by the lifting action of the outflow, then it is 
important to note that the relations between the boundary layer shear vector 
and the orientation of the lifting surface were almost ideal: the local shear 
vector was oriented from an azimuth of 135°, while the outflow was moving with 
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a vector almost perpendicular to the shear, providing maximum vorticity pro- 
duction all along the length of the front. 

Vortices qualitatively similar in appearance to those observed here have 
been observed by Carbone (1982, 1983) along an intense winter-time cold front 
in California. In Carbone's data, peak vorticities and divergences were simi- 
lar to those observed here, about lCT^s -1 , but overall wind speeds were much 
higher (e.g., storm motions of 28 m s"*), and the storm inflow was charac- 
terized by a strong jet parallel to the front. This jet was most intense at 
an altitude of 600 m AGL, and was hypothesized to be a storm-induced inflow 
feature. 

The vortices observed by Carbone apparently had typical scales of about 
13 km, somewhat larger than the 3-7 km observed here. Carbone attributed the 
vortices to Helmholtz instabilities in the regions of strong horizontal shear 
along the front. Barcilon and Drazin (1972) conclude, using standard linear 
stability analysis techniques, that shear zones which can be modeled as verti- 
cal vortex sheets are always unstable in the presence of superadiabatic lapse 
rates. Evidence for the presence of superadiabatic lapse rates near the out- 
flow shown in these data is described below. Thus, it is possible that 
Helmholtz instabilities may have been responsible for the development of the 
vortical circulations observed in the lidar data discussed here. Vortices 
resembling Carbone's and those observed here have also been obtained along 
outflow shear zones in numerical simulations of convective clouds (Rotunno, 
1984, private communication). 

It is especially important to note that the vortices observed here are 
strong enough to appear in the plots of total measured wind vectors, rather 
than merely in the plots of perturbation vectors. The lingering presence of 
systematic errors in the lidar measurements was seen in places in the pertur- 
bation wind fields and in contour plots of certain kinematical quantities 
inferred from the lidar winds. One important conclusion to be made from this 
study is that while the winds derived from the lidar measurements are only 
subject to relatively minor errors (usually about 1 m s"*), these errors are 
large enough to render suspicious many of the features seen in the perturba- 
tion wind fields. 

Unfortunately, it was not possible to confirm the presence of these vor- 
tical patterns by an independent observation from Doppler radar due to the 
entry of the gust front into the region of ground clutter by the time the CIM 
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Doppler radar was turned on. Details of the gust front which lay outside of 
the ground clutter appeared noisy and were apparently subject to sidelobe con- 
tamination from storm reflectivity in adjacent beam radials. However, the 
strength of the winds in the outflow, about 20 m s”* in spots, was observed to 
be similar to the CIM radial velocities reported in the core of the storm some 
10 km southwest. As mentioned earlier, such velocities would probably be good 
estimates of the wind due to the approximate alignment of the mean winds in 
the outflow layer toward CIM. 

Further confirmation of the strength of the outflow came when the front 
arrived at TTS. This station reported passage of the gust front at about 
202230 GMT (142230 CST), with winds veering to southwest and increasing to 
maximum values of 12.5 m s - * immediately after gust front passage. The 

surface winds at TTS at the time of gust front passage is shown in Fig. 3.12, 

along with other major meteorological parameters recorded. This figure, 
designed to depict the time series of TTS data in a space cross-section which 
might be seen across the outflow looking north-northwest, reveals not only a 
veering and surging of the wind associated with it, but also shows a sharp 9°C 
drop in temperature and a 1°C drop in a dew point behind the outflow. How- 
ever, this temperature decline was evidently not directly associated with the 
onset of heavy rain which was delayed by about 14 minutes with respect to the 
time of passage of the gust front. The onset of heavy rain, which reached 
intensities of 120 mm/hr briefly, was also evidently accompanied by a second 
surge of strong wind, with a peak gust of 18.3 m s~* from the west- 

southwest. This lag between the first wind surge and the onset of heavy rain 

provides strong evidence that the “wall" of high intensity returns seen by the 

lidar in the vicinity of the gust front was associated only with the arcus 

cloud and not with heavy precipitation. 

Additional information about the structure of this outflow may be gleaned 
from the data using appropriate assumptions about frontal motion. If the 

speed of propagation of the front is taken, following Goff (1975), to be about 

0.67 times the peak surface gust value at the leading edge, then it is pos- 
sible to deduce a frontal propagation speed of about 0.67*12.5 = 8.3 m s _1 . 
Using the speed and the observed delay (about 150 s) between 2020 GMT and the 
time of frontal passage at TTS, it appears that the surface front was roughly 
150*8.3 = 1250 m west-southwest of TTS at 2020 GMT. The lidar raw data taken 
from forward frame 82 at 201955 GMT indicate the frontal position at 722 m 
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Figure 3.12 Plot of meteorological observations at TTS (Tuttle South), 2000- 
2100 GMT, 30 June 1981. Note passage of the outflow at 2023 GMT (three min 
after the CV 990 over flight) and delayed onset of heavy rain. The vector 
labeled 20 m s“* scales the wind vectors; the bold arrows are 60 sec. averaged 
winds whereas the thin superimposed arrows are the peak winds in each 60 sec. 
period. Wind direction is referenced to North which is toward the page top. 
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altitude to be in resolution volume 12, at a distance of 192 + 11*320 = 

3712 m. Even if we account for as much as 1 km error in aircraft position 
with respect to the front, it is still reasonable to conclude that the front 
sloped back to the west -southwest at a slope of about 722/ (3200-1250) = 

0.37. While the front almost surely did not maintain a simple wedge shape at 
its leading edge, we can at least deduce that much of the lidar data collected 
in run 2 was collected just above the leading several km of a sloping storm 
outflow airmass. Unfortunately, the relative positions of the surface front 
and the lidar observed front cannot be given along the entire length of the 
front because the appropriate surface data were not available. 

Comparison of surface temperature data with temperature data taken by the 
aircraft reveals the likely presence of a slightly superadiabatic lapse rate 
of 10.6°C per km of altitude at TTS location at 2020 GMT, just two and one- 
half minutes before the arrival of the outflow. 

Additional outflow data were taken during run 3 along the northern edge 
of the storms. Plots of the derived winds are shown in Fig. 3.13, Note the 
presence of wavelike patterns in the wind field intensity and also the sinuous 
shape of the leading edge of the outflow as observed earlier along run 2. 

These features may be manifestations of the cellular character of the convec- 
tion or may be related to the vortical perturbations seen further southeast in 
run 2. Indeed there appears to be some continuity to the form of the outflows 
in the northeast corner (beginning) of run 3 and the northwest corner (ending) 
of run 2, both showing strong southerly or south-southwesterly flow at about 
20 ms"-'-. However, at least some of the perturbations observed during run 3, 
especially those behind the leading edge of the outflow, are probably arti- 
facts of lidar and INS system instrumental or coordination errors which were 
not completely or adequately handled by the filtering methods employed in 
these analyses. 

A feature of significance seen in Fig. 3.13 is the sharply defined change 
in velocity located just southwest of the gust front by a distance of about 
3 km and oriented northwest -southeast. This is evidently caused by an 
especially blatant case of the lidar system error mentioned earlier, because 
it is seen to align itself quite closely with the direction of aft-pointing 
lidar lines of sight in run 3. This particular lidar error may be associated 
with the "shock" experienced by the aircraft when it actually penetrated a 
portion of the gust front at coordinates (X = -36, Y = 4). INS measured wind 
and drift angle data from this run confirm the impact the gusty outflow had on 
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Figure 3.13 Vector plot of lidar-derived winds observed during run 3. 
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the aircraft: INS measured winds fluctuated by as much as 6 m*s"^ in magni- 

tude, and drift angles changed by more than 2° between successive lidar frames 
during the gust front penetration. Plots of aircraft-measured temperature and 
altitude are also provided (Fig. 3.13) to demonstrate the penetration of the 
outflow by the aircraft and to suggest that the temperature changes measured 
at the aircraft were likely to be due to more than just the slight altitude 
change it experienced. 

Clearly the results obtained by the lidar during the period while it was 
being buffeted strongly must be viewed with extreme skepticism. Thus, some of 
the larger wind shifts seen southwest of the gust front's leading edge are 
probably erroneous. Only the lidar data collected prior to penetration of 
thegust front seem fully credible; fortunately the orientation of the front 
was such that a relatively good view of the frontal wind structure was 
obtained while the aircraft was still in smooth air. 

Fig. 3.14 shows the field of reflectivity. The only significant features 
in the intensity field are noted near and just south of the outflow's leading 
edge. This small band of high intensity returns may be due to precipitation 
or to scud filaments which were observed along this portion of the storm 
(e.g., Fig. 3.5). No organized arcus cloud was observed during run 3. 

In Fig. 3.15 the spectral width is seen to have a well defined maximum at 
coordinates (X * -32, Y = 4), which correspond closely to the position of a 
portion of the leading edge of the outflow where it is strongly curved and 
apparently characterized by convergence and cyclonic vorticity. Also evident 
in Fig. 3.15 are the noisy widths, whose origins were explained earlier, in 
the nearest range bin to the aircraft. 

Based on the observations made in run 3 and confirmation from similar 
studies of the results from other runs, it is possible to make several general 
statements about the systematic errors found in the lidar data and their 
impact on the results of this study. First, wind shifts which show a prefer- 
ential orientation along either the fore or aft lidar lines-of -sight - in run 
3 (either 170° and 130°) must be considered suspect. Second, the presence and 
the intensity of the system errors seem to be quite sensitive to the relative 
motion of the aircraft and the ambient wind, with the worst errors noted when 
the aircraft was subject to a crosswind. Third, kinematical parameters, such 
as divergence, vorticity, and stretching and shearing deformation, are quite 
sensitive to the presence of these errors. 
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Figure 3.14 Contour plot of normalized reflectivity observed during run 3. 
Note relative lack of strong intensity areas near the gust front, due to the 
absence of the arcus cloud formation. 
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Figure 3.15 Contour plot of spectral widths (turbulence) observed during 
run 3. 



3.5.2. Interpretation of Data at the Edge of an Isolated Cumulus Congestus 

During runs 5 through 12, lidar observations were made of an isolated 
cumulus congestus near Elmore City (ELC), OK (Fig. 3.2), Observations were 
made at 1100 m AGL and at 2300 m AGL, the former being close to cloud base and 
the latter well above cloud base. Fields of velocity vectors, reflectivity, 
and spectral width were generated using an advection vector given by the mean 
of the four run-averaged wind vectors used in each circumnavigation. Thus, 
the advection vector for the lower level was not the same as that used for the 
upper level. Time differences between observations at each level were less 
than 350 s, but were as much as 750 s between levels, rendering the use of one 
common advection vector and objective analysis time inappropriate in view of 
the obvious convective growth which was taking place in the region of 
interest. 

The lidar range was about 6-8 km during the lower altitude passes and 
increased to more than 12 km during the higher altitude passes. This increase 
in range was, of course, probably attributable to the decrease in ambient 
water vapor with altitude. 

Figs. 3.16 and 3.17 show the lidar-deri ved wind vectors computed by com- 
bining forward and aft lidar data from within single runs only, one run at a 
time for altitudes of 1.1 and 2.3 km respectively. However, for small scale 
circumnavigation studies such as those involving runs 5-8 and 9-12, further 
information about the wind field, especially in the corner regions not well 
covered in the standard analysis (i.e., combining fore and aft data only in 
each run), can be obtained by combining lidar data from pairs of different but 
adjacent runs. In such analyses, for instance, wind vectors were computed 
using lidar data from forward frames in run 5 and aft frames in run 6, similar 
analyses were made for forward and aft data from runs 6 and 7, 7 and 8, and 8 
and 5 respectively (Fig. 3.18) . 

It should be noted that several forward and aft frames were deleted at 
the beginning of run 11 (at about X = 29, Y = -77) due to excessive roll by 
the aircraft. These deletions have produced a gap in the analysis shaped like 
an inverted "V" (Fig. 3.17). 

Note that in plots 3.16, 3.17 and 3.18 the superimposition of vectors 
computed from the various combinations of runs offered an opportunity to 
examine the lidar analyses for consistency and time continuity. In Fig. 3.16, 
agreement of vectors in areas of overlap between data collected in the four 
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Figure 3.16 Vector plot of lidar-derived winds observed during runs 5-8 (each 
run taken singly). The data were taken just below and around a towering 
cumulus at 1130 m AGL. 
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u 0 = 4.26 ms -1 
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WIND FIELD AROUND CUMULUS CONGESTUS 
30 JUNE 1981, 205640 GMT 
Z = 2290 m AGL 


X-DISTANCE FROM C1M (km) 


Figure 3.17 Vector plot of lidar-derived winds observed during runs 9-12 
(each run taken singly). The data were taken at mid-levels at an altitude of 
2290 m AGL. The dashed line is the outline of high reflectivity regions shown 
on Fig. 3.21. 
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Figure 3.18 Vector plot of lidar-derived winds observed during runs 5-8 
(using pairs of adjacent runs' data in the wind synthesis). 


RUN ? 



runs is variable, but consider the fact that the flow under examination is at 
the base of an active towering cumulus and that the time difference between 
observations used to construct the wind vector is as large as 90 s. Vector 
agreement ranges from near col linearity at coordinates (X = 27, Y = -88) to 
40° azimuth discrepancy at (X = 27, Y = -92). These differences are high- 
lighted in Fig. 3.19, which shows the perturbation vectors calculated relative 
to the assumed advection vector (see “UA“ and "VA" in the plot label) used for 
runs 5-8. The reasons for this and other disagreements are not known, but the 
actual time evolution of the features may be responsible. Another possibility 
is systematic instrument error in the lidar system, which as we have seen, 
hasthe capability to interfere significantly with the proper retrieval of the 
perturbation windfield. 

Accurate estimates of cloud-scale deformation, divergence, and vorticity 
can be made with an accuracy better than lO^s - * on a 250 m grid only if the 
velocities are obtained to within 0,25 m s~* accuracy. For these data, excur- 
sions, and presumably errors in drift angle measurements were on the order of 
0.2°, leading to uncertainties of about 0.5 m s~* in the retrieved wind data; 
for larger drift angle errors, large errors in retrieved winds are unavoid- 
able. Hence, many of the computed fields of divergence and vorticity obtained 
from the lidar data tend to be dominated by the presence of these systematic 
errors. 

Despite the differences in time and altitude, the two circumnavigations 
reveal the presence of several cloud-scale zones of convergence. Although 
there was considerable variation in the local pattern (Figs. 3.19 and 3.20), 
the general character of the environment near the towering cumulus was seen to 
be one of convergence at both levels. This would tend to confirm the observa- 
tions of Byers and Braham (1949), who found convergence in active thunder- 
storms often extended upward from above the low-level outflow up to beyond the 
mid-troposphere. This represents a substantially different type of cumulus 
cloud than that observed by Raymond and Wilkening (1982) in the New Mexico 
mountains, who found convergence below cloud base and divergence above it up 
to at least 6 km above mean sea level and inferred the presence of strong 
cloud-top mixing as the dominant entrainment mechanism at work. It should be 
noted that although our data show considerable evidence in support of the 
presence of lateral entrainment, the additional presence of cloud-top entrain- 
ment cannot be ruled out. 
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Figure 3.19 Vector plot of perturbation winds at 1130 m AGL 
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Figure 3.20 Vector plot of perturbation winds at 2290 AGL 
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This particular towering cumulus cloud failed to grow into an active 
thunderstorm; tops of the cloud were observed to tilt toward the east and 
southeast before reaching a level where an anvil could be produced or signifi- 
cant precipitation generated. Such behavior is inconsistent with the observed 
wind hodographs taken at OKC and TIK earlier in the day, which showed weak 
winds and weak shear at all levels. Entrainment, which quite likely prevailed 
at cloud edge based on the prevailing convergence observed at both levels 
examined, surely played a role in reducing the buoyancy of the updrafts. 

The approximate location of the edge of the cumulus may be deduced from 
an inspection of Fig. 3.21 which shows the pattern of reflectivity from 
thelidar data taken during runs 9-12. This pattern of reflectivity should be 
compared with that shown on Fig. 3.22 which is for data below cloud base and 
which does not show any high reflectivity boundaries. Although some data 
voids are present in Fig. 3.21, the edge of the cloud may be estimated to 
follow the dashed line which encloses all the anomalously high intensities. 

At this point it may be worthwhile comparing Fig. 3.21 with the photo of the 
cloud shown on Fig. 3.6, taken during run 7. Both figures show that the cloud 
is ragged, suggesting that the cloud is being eroded, probably by the entrain- 
ment of drier air above the mixed layer. This entrainment is consistent with 
the general pattern of convergence observed (Figs. 3.13 and 3.14). Also 
important is the fact that, from comparison of Figs. 3.16 and 3.17, the velo- 
cities at the higher altitude are slightly smaller than those at the lower 
altitude and suggest veering with height. Thus the shear vector averaged over 
the area was from the northeast, which may account for the northeast alignment 
of the cloud. 

Questions still remain as to why the larger thunderstorms were not sub- 
ject to the unfavorable mechanisms which prevented the observed cumulus cloud 
from growing into a thunderstorm. Although they were closer to the main fron- 
tal boundary, it seems hard to support the view that 40 km of distance from 
such a weak boundary could make such a large difference in the fate of the 
convective elements. Such increased distances from the front did not prevent 
Ft. Sill, in southwestern Oklahoma, from experiencing a heavy thundershower 
early on the afternoon of 30 June. Further observations with both lidar and 
radar, including more time continuity studies, along with temperature and 
pressure measurements and numerical simulation studies of convective clouds in 
specified environments will be needed to resolve the questions related to the 
success or failure of individual cumulus clouds to mature into thunderstorms, 
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Figure 3.21 Contour plot of normalized reflectivity at 2290 m AGL 
strong reflectivity gradients at cloud edges (dashed line). 








and to resolve the questions regarding the evolution of ordinary thunderstorn 
into severe thunderstorms. 

3.6. CONCLUSIONS AND RECOMMENDATIONS 

By the use of airborne Doppler lidar it has been shown that significant 

information not readily available from any other source can be obtained about 

the flow field in the storm environment. Details of the low-level outflow 
were mapped for a large multicell thunderstorm system which occurred in Okla- 
homa on 30 June 1981. In addition, ambient wind flow at cloud base and aloft 

was studied for an isolated cumulus congestus which occurred that same day. 

In the case of the outflow study, waves and vortices were noted along th 
leading edge of the flow, especially where the outflow air propagated perpen- 
dicular to the ambient low-level shear vector. The vortices appeared to be 
due to tilting of vertical -shear-related horizontal vortex sheets by the dif- 
ferential rising motions ahead of the outflow, to convergence of horizontal 
shear at the frontal boundary, and possibly to Helmholtz instability of the 
sheared interface between the warm air and the cool outflow. 

The flows around the cumulus congestus were generally characterized by 
cloud-scale convergence at cloud base and at mid-levels, with a number of con 
fluence lines outside the cloud and divergent flow patterns of a variety of 
scales at both levels. Flows at mid levels near the cloud edge appeared to bi 
a reflection of the expansion or contraction of active cloud turrets. 

An intercomparison with single Doppler radar was attempted for that por- 
tion of the lidar data most suitable for such a study, but noisiness and poor 
resolution in the radar data prevented any definite conclusions from being 
drawn other than that the instruments agreed to within 1 m*s _1 on the average 
velocities retrieved. 

Most of the meteorological data fields produced by the Doppler lidar mea 
surements seemed consistent and believable, but there remained evidence, even 
after intensive data editing, of some systematic errors in the data. These 
errors were apparently related in a complex way to poorly understood and 
inadequately measured fluctuations in the aircraft attitude during data col- 
lection. Although the basic wind vectors derived from the Doppler lidar data 
seemed to give a satisfactory and often revealing view of near-cloud flow pat- 
terns, systematic errors often prevented clear-cut views of the smaller per- 
turbations in the windfield. 



To take full advantage of the airborne Doppler Ildar's spatial resolving 
capability, Doppler velocities must be obtained with a 0.50 m s" 1 accuracy in 
order to estimate cumulus-scale (i.e., 500m) divergence and vorticity, which 
frequently have magnitudes of 10"^ s"*. For these data, even the quantization 
errors in the drift angle measurement on the order of 0.2 deg generated uncer- 
tainties of about 0.5 m s“* in the retrieved wind data; because larger drift 
angle errors also occurred, larger errors in retrieved winds, sometimes on the 
order of 3-4 m s~*, were observed. Hence many of the fields of divergence and 
vorticity obtained from the lidar data tend to be dominated by the presence of 
these systematic errors. 

One explanation for the reason the calculated kinematical parameters 
proved to be so sensitive to lidar errors has to do with the objective analy- 
sis scheme used in the interpolation of raw data to a regular grid. The 
Cressman weighting function used in the objective analyses of these lidar data 
was circularly symmetric. Such a weighting function is best suited for use in 
data fields which are rather evenly distributed. Although in general the 
lidar data network satisfied this requirement for evenness of distribution, 
there were notable instances where lidar data along certain noise-contaminated 
frames had to be deleted. These missing frames tended to occur, of course, in 
just those parts of the flight characterized by significant velocity fluctua- 
tion. The absence of data along these frame radial s caused the distribution 
of raw data points within the circular Cressman regions of influence around 
some grid points to be quite uneven. The uneven distributions of data within 
the circles of influence, along with the exaggerated error characteristics of 
the raw data obtained during turbulent flight, combined to produce especially 
biased estimates of gridded quantities near the "missing" radials. Applica- 
tion of finite differencing to pairs of gridded estimates which happened to be 
biased in opposite ways was one way in which the analysis technique used here 
might have amplified existing errors. 

In future work with this sort of data, an interpolation scheme which can 
deal effectively with asymmetric data distributions in the vicinity of gra- 
dients should be used. The interpolation scheme will not be urgently needed 
if steps are taken to reduce the laser's susceptibility to "moding" during 
turbulent flight, and if the roll angle compensation in the lidar optics is 
operating properly. 



Whenever the aircraft is experiencing large drift angle perturbations, 
systematic velocity errors (frame bias) may also enter the data even after 
applying the algorithm designed to remove frame biases if: (a) the duration 

of the lidar frame is long with respect to the period of drift angle perturba- 
tion, (b) the drift angle perturbation is of a form which allows significant 
variations in aircraft groundspeed to enter the actual lidar radial velo- 
cities, and (c) the values of drift angle and groundspeed assigned to the 
lidar frame do not correspond to the properly weighted values needed for an 
accurate removal of the aircraft groundspeed component. In run 3, a possible 
instance of this problem occurred when the aircraft penetrated the outflow 
which came from the storm under surveillance. Differences of 2° were noted 
between successive measurements or drift angle along with 6 m s - * changes in 
1 idar-derived radial velocity between successive like-directed lidar frames. 
These sharp changes could be associated with significant systematic error in 
the lidar-derived velocities, which were obtained at each range by averaging 
the returns from some 50 pulses fired during a 0.5 s period. This problem 
might be present to a lesser extent whenever the aircraft is flying through 
air containing milder turbulence. 

In order to eliminate this problem, it would be necessary to remove the 
aircraft's time-varying ground-relative component of Doppler shift in the echo 
of each transmitted pulse prior to the averaging which provides the estimate 
of radial velocities along a frame. If this proposed method of eliminating 
the problem is too burdensome to be compatible with real-time hard-wired data 
processing, then an alternative approach would be to record more frequently 
the drift angle, heading, and groundspeed. Further study of this potential 
source of error seems warranted to determine its importance for future work 
with this type of data. 

Another potentially significant source of error that could affect the 
computation of lidar-derived wind vectors is that related to possible dif- 
ferences in altitude of the fore and aft data used in solving for the wind 
vectors. Due to the beam pointing errors which the lidar system experienced 
as a result of the scanner's roll and pitch compensation problems, there was a 
1-2° uncertainty in the beam elevation angles. Obviously, fluctuations in 
lidar beam elevation angle will have their most injurious effect on retrieval 
of wind vectors at large ranges from the aircraft and in situations charac- 
terized by strong vertical shear. In fact, for a range of 10 km from the air- 
craft, the use of two different lidar frames having elevation angles differing 



by 2° implies use of radial wind data from levels separated by about 350 m. 
When vertical shear is strong (e.g., 10“^ s"*) this difference can produce 
errors comparable to the errors produced by the drift angle delay. 

A conservative estimate of the wind speed error caused by beam elevation 
errors is 0.5 m s"* except above the outflow on runs 2 and 3 and along the 
cloud turrets of runs 9-12, where local errors appeared to reach 2-3 m s - *. 
Clearly, the problems of laser misalignment and failure of the scanner to 
provide the proper beam elevation compensation during aircraft rolls should be 
corrected before the next series of experiments. In addition, it would 
bedesirable if scanner positioning could be updated at least occasionally 
during each lidar frame so as to ensure horizontal beam trajectories for all 
pulses. 

Another source of error contaminating the lidar data is that associated 
with uncertainty of aircraft groundspeed. If delays are present in the 
reporting of data during periods when the aircraft is neither accelerating nor 
decelerating, analysis shows that ground velocity uncertainty contributes only 
about 0.3 m s" 1 to radial velocity errors in the lidar wind data. However, if 
the aircraft is accelerating or decelerating, then ground velocity errors can 
be as large as the product of the acceleration and the delay. This can be a 
significant source of error during reconstruction of the wind fields. 

Several other problems were encountered during the data anslysis with 
respect to improper sampling or recording of various aircraft motion and 
attitude parameters. Improved accuracy and detail in the recording of lati- 
tudes, longitudes, velocities, angles and times will be helpful in future work 
with the data from the NASA's airborne Dopier lidar. In view of the magni- 
tudes of some of the velocity errors found in some parts of the flight, 
special care should be taken to insure that a detailed record of the INS per- 
formance is maintained. It is especially important that all portions of data 
collection flights be documented with nadir photography, so that the aircraft 
positions computed by the INS can be compared with the actual positions 
obtained from detailed maps. The timing of the photographs should be perhaps 
slightly faster on short flight legs (such as those found in cloud circumnavi- 
gations) in order to get a sufficiently large number of photos per run. 

One possible alternative to the use of nadir photography in determining 
whether or not the INS is furnishing erroneous data involves a scheme in which 
ground strike data are obtained after every fore/aft frame pair and the mea- 



sured Doppler velocity of the ground compared with its expected value of 
zero. While this idea shows considerable promise for providing a direct way 
of measuring the INS errors, the nadir photography should not be dispensed 
with until the merits of any newly proposed technique have been established 
experimentally. 

With respect to the determination of the delays which tend to occur in 
the data recording process, one possible way of studying this problem would be 
to fly a specially designed test track at both the beginning and the end of 
each day's data collection effort. The track would consist of level flight, 
first at a constant but slow ground speed of say 120 m s~* for a period of 
about 30 s; this would be followed by a 30 s period of steady acceleration to 
a ground speed of 150 m s'"*, then followed by 30 s of flight at that high 
velocity; finally another 30 s period of steady deceleration back to 120 m s~^ 
would ensue. After another 30 s of straight flight at 120 m s"*, a final 60-s 
leg featuring shallow weaving and banking maneuvers to the left and the right 
would be made, with banking in each direction lasting 10 s or so. Nadir 
photographs would be taken every 3-5 s to provide adequate coverage of each 
leg of the test run, and the lidar would operate in its normal data collection 
mode. Inspection of the wind field analyses from such test runs would be 
undoubtedly helpful in showing whether any errors or inconsistencies are 
affecting the data. 

In summary, the lidar thus far has demonstrated considerable promise in 
portraying the details of the clear air flows near an assortment of convective 
clouds. Not only can the wind fields be obtained in the clear air near cloud 
edges, but the probable locations of those cloud edges can be determined from 
inspection of the echo intensities. Doppler spectral widths often are useful 
in locating wind field shifts and perturbations. Revisions to the lidar sys- 
tem which are currently underway will probably correct most of the problems 
noted above, and result in even more detailed and revealing insights of the 
flows near convective clouds. 
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APPENDIX A 


DESCRIPTION OF INSTRUMENTATION 


\.l Doppler Lidar 

A description of the lidar system has been given in Bilbro, et al . 

(1984), nevertheless, certain features of the system merit discussion due to 
the novelty of the instrument's use in this type of experiment and the ways in 
tfhich system performance can affect data quality. Fig. A.l depicts in block 
jiagram form the system's main components and their interrelationships. 

The NASA lidar employs an 8 W continuous-wave C02 laser which emits 
infrared radiation at a wavelength of 10.6 un. This laser, which acts as the 
master oscillator in a master-oscillator power amplifier system, sends its 
output to a beam splitter, which transmits most of the infrared beam along a 
pathway toward the atmosphere, and diverts the small remaining portion for use 
as a reference frequency during detection of the Doppler shifted echoes. 

That portion of the beam which is bound for the atmosphere is first sent 
to an electro-optic modulator where it is converted from a vertically 
polarized continuous wave into a train of horizontally polarized 2.0 us pulses 
at a rate of fp of 100 Hz. The electro-optic modulator consists of a mercury 
cadmium telluride crystal flanked by polarizing filters of vertical and hori- 
zontal orientation at its input and output sides respectively. The actual 
conversion o? the CW radiation into a pulse train with new polarization occurs 
because the crystal, under the influence of a transverse electric field 
(applied for the specified pulse duration at the specified pulse repetition 
rate), becomes birefringent and causes the vertically polarized radiation from 
the input filter to acquire horizontal polarization during passage through the 
crystal, thus permitting unobstructed passage through the output filter only 
during those time intervals designated for pulse transmission. The pulses 
produced by the electro-optic modulator then pass through a Faraday optical 
isolator which prevents scattered radiation from re-entering the master laser 
and causing any frequency shifts ("pulling") in that Laser's output. 

Next the pulses enter a beam expander and a series of amplifiers where 
they experience some 36 dB of gain, giving them typical energies of about 
15 mJ per pulse. The pulses then pass through a Brewster window aligned to 




Figure A.l Block diagram of NASA's Airborne COg Doppler Lidar System 
(ADLS) , BW is the width of the beam. Symbols are defined in the text 














pass only horizontally polarized radiation. Next the pulses pass through a 
quarter-wave plate made of transparent material having an index of refraction 
which is a function of the polarization orientation. The quarter-wave plate 
is oriented so that the mutually perpendicular axes for which indices of 
refraction are minimal and maximal are both 45° from the horizontally oriented 
axis of polarization of the outbound laser radiation. Thus, after traversing 
the quarter-wave plate, one component of the initially horizontally polarized 
radiation lags its normal component by 90° of phase, and the resulting net 
radiation transmitted by the plate is circularly polarized. 

After passing through the quarter-wave plate the pulsed beam is expanded 
and collimated by an all-aluminum folded parabolic off-axis telescope and 
directed through a scanner which deflects the beam out an exit window on the 
left side of the aircraft and onto a specified path through the atmosphere. 
Because of the series of expansions the laser beam undergoes while on its way 
through the system optics, the width of the beam increases from about 8.5 mm 
at the master oscillator to about 20 cm at the point of entry into the atmos- 
phere, as measured between the (1/e^) points of a Gaussian function which 
approximates the power distribution across the beam. 

The scanner consists of a pair of separately rotatable germanium wedges, 
the angular orientations of which determine the beam's trajectory within a 
cone having a 19, 6° half -angle and axis normal to the left side of the air- 
craft. For this experiment the system was programmed so that the scanner 
wedges would, for any value of aircraft pitch and for all values of aircraft 
roll smaller than 9°, direct the beam in a true horizontal plane at the maxi- 
mum possible forward and aft deviations from the centerline of the cone, which 
itself always points normal to the left side of the aircraft. For roll angles 
greater than 9°, the wedges were programmed simply to direct the beams out the 
aircraft at the angle (relative to the aircraft) that would be chosen if the 
aircraft were assumed to be in roll-free flight. Thus, in the absence of sig- 
nificant aircraft roll, the pulses were to be horizontally directed approxi- 
mately 20° fore or aft of a line extending from the left side of the aircraft 
and having an orientation normal to the aircraft heading vector (Fig. 1.1), 

Lee (1982) found evidence from some lidar data taken in California in 
1981 that the scanner was directing the lidar beam slightly out of the hori- 
zontal plane during quasi -roll -free flight. He measured elevation biases of 
about -0.4° for aft shots and +1.0 or more for fore shots. Because such beam 
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elevation biases can lead to significant differences in the altitudes of the 
various resolution volumes of data used in the objective analyses, attempts 
were made in this study to confirm whether or not the scanner was operating 
according to specification. The technique used to assess scanner performance 
and the results and implications of that assessment are discussed in Appendix 
E. 

Typical operation of the data sampling cycle in this experiment involved 
the transmission of about 50 pulses, requiring about 0.5 s in the forward 
direction, followed by a time gap of at least 0.6 s for the reorientation of 
the wedges, followed by the transmission of about 50 pulses in the aft direc- 
tion, and again reorientation of the wedges. This pattern was repeated more 
or less continuously for the several minutes of each of the many data runs. 

For the purposes of subsequent discussion, each distinct sequence of 50 simi- 
larly directed transmitted pulses is denoted by the term "frame". The wedges 
remain locked in position during the interval in which the pulses of each 
frame are transmitted. In this experiment, the period between frames varied 
from 1.1 to 1.5 s. Although generally each frame was followed by a counter- 
directed frame, exceptions to this rule were occasionally noted. 

The sense of rotation of the polarization vector of the backscattered 
radiation, viewed along the direction of propagation, is reversed by 
spherically symmetric scatterers (Newell et al., 1955). Thus, the radiation 
which was originally right-circularly polarized prior to being scattered 
becomes left -circularly polarized after being scattered (see Fig. A.l). This 
returned radiation is also characterized by a Doppler shift in frequency given 
by f d = -2V p /x, where \ is the wavelength of the transmitted radiation and V r 
is the target's radial velocity (positive when motion is away from the lidar) 
relative to the aircraft. The received signal (quasi -continuous in time if 
aerosols are quasi -uniformly distributed in space) when sampled at any range- 
time delay (i.e., time interval following the transmission of the pulse) gives 
a sample amplitude proportional to the sum of the incoherent echoes from all 
scatterers within the resolution volume associated with that particular sample 
(Doviak and Zrnic', Chapter 4, 1984). 

After returning to the telescope, the Doppler-shifted "echoes" travel 
back through the quarter-wave plate where they have their polarization con- 
verted from circular to planar (Fig. A.l). During this final change of 
polarization, the lagging component of the circularly polarized backscattered 



wave experiences a further retardation of phase by 90°, causing the received 
signals to acquire vertical polarization instead of the horizontal polariza- 
tion they originally possessed prior to traversing the quarter-wave plate on 
the outbound portion of their excursion into the atmosphere. The Brewster 
window which the returns encounter next thus reflects them toward a beam 
splitter which combines them with the original master oscillator signal and 
directs both to the surface of a cryogenically cooled mercury cadmium 
telluride detector. 

The detector takes the resulting time-varying interference patterns and 
converts them into an output signal having a modulation frequency equal to 
that of the Doppler shift f d of the echo signals. Fig. A. 2a shows a pair of 
transmitted pulses and (b) the interference pattern, at the detector input, 
caused by mixing the Doppler shifted echo signals (assumed to be returned from 
a uniform distribution of targets all having the same 10 m s“* velocity) with 
the master oscillator signal. The repetition time (PRT) of the pulses and the 
corresponding spatial separation are indicated on the figure as well as the 
approximate pulse shape, but the schematic exaggerates the period f t _ * of the 
radiation. The Doppler shift f d (r) associated with the echoes from the 
various resolution volumes is due to the aerosol target speed added to the 
aircraft airspeed along the lidar line-of -sight (L.O.S.). An aircraft 
airspeed of 150 m s" * causes about a ±10 MHz shift for the fore- or aft- 
directed beams. 

Although the optical detector cannot distinguish positive from negative 
Doppler shifts (that is the reason for absolute signs around f d in Fig. A.l), 
ambiguities in the sign of f d (r) are avoided because of the dominance of the 
Doppler shift due to aircraft airspeed along the lidar L.O.S., a shift whose 
sign can always be determined from a knowledge of the direction of the lidar 
scan. Now f d (r) may be expressed in terms of the Doppler shift due to the 
aircraft airspeed along the L.O.S., namely f d (0), and some superimposed range- 
dependent deviations Af d (r) from f d (0) due to the L.O.S. speed of aerosols 
relative to the L.O.S. speed of air at the aircraft. 

Following extraction of the Doppler-shifted signal by the detector, an 
amplified version of the detector's output is mixed with a signal having a 
frequency given by the sum of f co = 60 MHz and j f d ( 0 ) | , the latter being near 
10 MHz, to produce a difference frequency at an intermediate frequency fjp 
near 60 MHz but shifted slightly from it by the amount ±Af d (r) for the aft or 
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(a) TRANSMITTED PULSES 



(b) DOPPLER-SHIFTED ECHO, OPTICAL DETECTOR INPUT 
TYP RES. VOL 



Figure A. 2 Transmitted pulses (a) and the interference pattern (b) at the 
input to the optical detector assuming airspeed is zero and aerosols are 
uniformly distributed and all move at a radial velocity of 10 m s . (c) The 

inphase (or quadrature phase) signal input to the poly pulse pair processor's 
(PPPP) analog to digital (A/D) converter. The bold vertical lines indicate 
the spacings of signal samples digitized by the A/D Converter. 


100 


fore beams respectively. The 60 MHz signal is supplied by a coherent oscil- 
lator (COHO), while the signal having frequency |f cJ (0)| is supplied by a 
voltage-controlled oscillator (VCO) tuned to the most recent value of the air- 
craft airspeed component along the lidar line-of -sight. The signal 1^(0)! is 
added to f in the mixing process in order to keep the mixer's output near 60 
MHZ, and to obviate the need for estimating directly the inconveniently large 
Doppler shift f d (r), The aircraft's airspeed is measured by a probe and moni- 
tored by the Central Timing and Control System (CTCS) which provides commands 
to the VCO. 

The signal, now at fjp, is next sent to a poly-pulse-pair processor 
(PPPP) where it is mixed with inphase (0) and quadrature phase (t r /2) compo- 
nents of the COHO voltage in order to generate the complex video (I and Q) 
version of the range-dependent Doppler shift Af d ( r ) , the sign of which is 
determined by the relative phase of I and Q. Figure A. 2c shows the I or Q 
signal at the input of the PPPP or, if airspeed is zero, at the output of the 
optical detector. In the PPPP, the complex video, usually having a frequency 
of not more than several MHz, is converted from an analog signal to a digital 
sample at a 15 MHz rate and then subjected to the autocovariance and spectral 
processing which produces the estimates of the signal's first three moments: 
intensity (zeroth moment M Q ), mean frequency (first moment M^) and spectral 
width (the square root of the second moment M 2 ). 

The poly-pulse-pair processor represents a nat.ural extension of the 
familiar pulse-pair processor (PPP) employed in most Doppler weather radar 
systems. In the PPP, the mean frequency shift of the returns from a resolu- 
tion volume is estimated from the argument of the first lag of the auto- 
covariance function (acvf) of the digitized complex video signal. However, 
because of the smallness of the phase shifts that occur during the timespans 
of the pulse durations typically used by Doppler radar, the autocovariance 
estimation technique for such systems usually involves formation of lagged 
products of echo signal samples associated with sequentially transmitted 
pulses. In contrast, because of the larger Doppler phase shifts associated 
with the use of the lidar's higher transmitted frequencies, the calculations 
of autocovariance made by the lidar's poly-pulse pair processor is accom- 
plished using echo samples associated with individual transmitted pulse. 
Furthermore, because the lidar has a wavelength much shorter than the radar, 
echos from successive transmitted pulses are uncorrelated for the PRTs used. 
However, for each resolution volume the PPPP averages 32 lagged products at 
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each lag of the acvf examined in order to reduce the statistical fluctuations 
of the estimated acvf. Although 32 contiguous samples (spaced x $ = 66.7 ns 
apart) for each lag product (Fig. A. 2c) are averaged to reduce the statistical 
uncertainty in the covariance estimates, the equivalent number of independent 
samples effective in reducing variance is significantly less than 32 for large 
signal to noise ratios (SNR) because the contiguous lag product samples, when 
spaced within the time span of the transmitted pulse, are highly correlated 
(Doviak and Zrnic', 1984). However, for low SNR samples of signal plus noise 
are not highly correlated and hence a large number of nearly independent 
estimates are obtained to reduce variance in the estimate of the acvf. 

The PPP uses only the zeroth and first lag of the acvf to estimate the 
mean frequency and width of the Doppler spectrum, whereas the PPPP 
estimates the first 8 lag (i.e. lx <n x <8t ) products of the acvf in order to 

O v 5 

improve the accuracy of thee spectral moment estimates. The acvf is then 
Fourier transformed after setting to zero those values of the acvf that are 
not estimated at longer lags (i.e., 1024x s >nx s >9x $ ). This zero weighted 
extension of the lag domain increases the density of lines in the spectral 
domain so that the peak of the Doppler spectrum is better resolved. In NASA's 
PPPP the 8 lag acvf is unweighted before the Fourier transformation is per- 
formed but it has been shown (Lee and Lee, 1980) that an optimum weighting is 
obtained when each lag product is multiplied by the magnitude of the acvf at 
its respective lag before the Fourier transformation. The first moment esti- 

A 

mate Mj (the caret symbol is used to signify an estimate) is then the Doppler 
shift of the line which has the largest power. The intensity M q of this line 
is also recorded in order to estimate echo power given the estimate Mg and 
noise power, as well as the assumption of a Gaussian shaped Doppler spec- 
trum. However, no attempt was made to incorporate the second moment esti- 

A 

mate Mg in the computation of echo power. The truncation of the acvf 
at 8x g generates a spectral window that is usually much broader than the width 
of the Doppler spectrum associated with turbulent velocities inside the reso- 
lution volume. In this case, the intensity of the peak spectral line is pro- 
portional to the signal plus noise power. Although the zeroth lag product is 
not computed, errors in power estimates are usually small (i.e., <ldB) when 
the Doppler spectral width is small compared to the width of the spectral 
window. The second moment M 2 is estimated by least-squares fitting the 
measured autocovariance magnitude to the assumed Gaussian function. The noise 
power may be estimated by examining the behavior of the peak line intensity as 



a function of range where targets are not detected and may be removed from 
each intensity estimate via an algorithm described in Appendix D.l. The actual 
estimates of intensity, velocity and width assigned to a given resolution 
volume are the average values in the appropriate “range bin" (range interval 
corresponding to the resolution volume, numbered and ordered with respect to 
distance from the aircraft) obtained from processing echo samples from each of 
the 50 transmitted pulses that comprise each frame of moment data. 

After calculating the spectral moments in all resolution volumes of a 
given data frame, the PPPP sends the results to an on-board minicomputer which 
computes the target's Doppler velocity relative to ground using the aircraft 
ground velocity and airspeed estimates derived from the Inertial Navigation 
System (INS). Fig. A. 3 depicts the relation between the various components of 
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ai rcraft -measured velocities and the PPPP estimate M^. This figure is drawn 
for the hypothetical case in which the Ildar is measuring wind at zero range 
(i.e., at the aircraft itself). In the case depicted, the true airspeed vec- 
tor equals the sum of the ground velocity and wind vectors, but the drawing is 
equally valid if the wind vector is measured at nonzero range, however, in 
this latter case the sum of the ground velocity and wind vectors does not 
equal the true airspeed vector. Sensors measure the airspeed along the 
heading direction but the true airspeed vector is the vector sum of the mea- 
sured airspeed and the speed of air perpendicular to the heading ("side- 
slip"), Although sideslip is not measured by the aircraft sensors, the accu- 
racy of the lidar-deri ved radial wind data is not compromised, since any 
sideslip velocity components are automatically included in the lidar's radial 
velocity estimates. However, the INS estimates of wind at the aircraft will 
be in error if the sideslip is nonzero. 

For accurate retrieval of the Doppler velocities, the VCO (Fig. A.l) must 
accurately track the measured airspeed such that f (J (0) always equals: 

I f d ( 0) | = 2 (REF V r )A (A.l) 

for the fore and aft radial components of measured airspeed which establish 
reference Doppler velocities "REF V r ", However, it is not crticial that the 
VCO frequency satisfy Eq. (A.l) exactly, since the L.O.S. airspeed component 
REF V p removed from the raw Doppler signal during generation of f j p is 
restored intact during subsequent calculations (see Eq. (A. 2) below) in which 
account is taken of the L.O.S. components of aircraft airspeed and ground 
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Figure A. 3 The various velocity components used to derive the horizontal wind 
vector from lidar measurements. Thus schematic assumes the measurement of 
wind is made at r=0 (i.e., at the aircraft) where the sum of aircraft's ground 
velocity and wind equals the true airspeed vector. 
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velocity. As mentioned earlier the REF V p is used in order to place the 
observed Doppler shifts into the bandwidth of the PPPP, which has a Nyquist 
velocity V a of 39.8 m s -1 . Because the fore and aft data are not acquired 
simultaneously the fore and aft REF V r are generally not equal to each other 
as they are in Fig. A. 3. 

The quantity aV^ is the velocity to be estimated by the PPPP 
(i.e., A V r = AMj/2. Using Fig. A. 3 it is easy to deduce the relations 
between the fore and aft radial components of the wind and those measuredquan- 

A 

tities AV , INS-estimated ground velocity, and REF V_. The equations deduced 
r * 

from this figure are. 


r,a ■ ‘Va ' < REF V r ' 

(A. 2a) 

r.f ' "r.f + < REF V r ' V> 

(A. 2b) 


where V p a and V r ^ are the L.O.S. components of the wind vector along the aft 
and fore directions respectively, and V g>a and V g>f are the L.O.S. components 
of the aircraft ground velocity in those respective directions. These equa- 
tions are equally valid for measurements of wind far from the aircraft where 
the sum of ground velocity and wind does not equal the true airspeed vector. 

The digital sampling rate, 15 MHz, used in the PPPP is such that the 
highest Doppler frequency unambiguous to the processor, 7.5 MHz, corresponds 
to a Doppler velocity of 39.8 m s"* relative to the measured aircraft airspeed 
REF V p along the lidar line-of -sight (see Fig. A. 3). Signals having Doppler 
frequencies larger than 7.5 MHz relative to f d (Q) are filtered at the input 
side of the PPPP's digital converter. Nevertheless, ground-relative radial 
wind velocities of magnitude greater than 39.8 m s"* can be measured by the 
system if V g differs from REF V p such that AV p remains within the Nyquist 
interval (-V a ,+V a ). However, if the wind vector produces a L.O.S. velo- 
city aV larger than 39.8 m s"* an accurate estimate of aV will not be made 
r r 

because the signal with the Doppler frequency implied by such a large devia- 
tion will be filtered at the input to the PPPP and then only noise will remain 
to be analyzed by the processor. This has implications for the strategies 
which will be needed if the ADLS is to be used successfully to study such 
intense phenomena as tornadoes and tropical cyclones, where the variations of 
radial velocity along certain 1 ines -of -sight could exceed 39.8 m s _1 . On the 
other hand, it also indicates the importance of tracking REF V r w -j^h t ^e VCO 



in order that the Doppler-shifted echoes do not fall outside the Nyquist 
interval . 

A diagram showing the typical arrangement of resolution volumes in a dat 
frame and of lines of sight associated with successive frames is given in 
Fig. A. 4. Networks of fore- and aft-directed data are generated by the pulsei 
lidar beam as it scans the area to the left side of the aircraft. Both net- 
works are presumed to lie in the same horizontal plane in the data analyses 
described in the main body of this report. Because many of the lidar lines o 
sight deviated slightly from true horizontal due to aircraft roll and surveyoi 
error, some of the fore and aft resolution volumes which appear to overlap 
when viewed in plan actually differ in altitude by amounts which increase witl 
range. This point receives further discussion in Section 3.4.3. A tabular 
listing of all important parameters governing the operation and performance o 
the Doppler lidar system is given in Table A.l. 

A. 2 DOPPLER RADAR 

The National Severe Storms Laboratory (NSSL) operates two 10 cm Doppler 
radars, one located in Norman (NRO) and the other at Page Airport (CIM, for- 
merly Cimarron airport) 40 km northwest of NRO (Fig. 3.2). Parameters from 
these Doppler radars are compared with parameters of the lidar system in Tabl< 
A.l. The obvious advantage of the lidar system is that it has excellent angu 
lar resolution and low power requirements. However, its range is limited and 
because of strong atmospheric attenuation, increases in transmitted power do 
not increase range in direct proportion to the increase of transmitted power. 

The Doppler radar operation principles and methods have been described 
elsewhere in detail; see for example Doviak and Zrnic' (1984), or Doviak 
et al (1979). 

A. 3 OTHER INSTRUMENTS ON BOARD THE AIRCRAFT 

In addition to the lidar itself, a number of additional instruments on 
board the CV-990 took measurements which were important to the success of the 
experiments. For the sake of completeness, those instruments which furnished 
data mentioned (excluding cameras), are briefly described here. Where avail- 
able, information regarding the resolution and probable accuracy of the 
various system components is included, for further details see Bilbro (1982) 
and Ames Research Center (1981). 




Figure A. 4 Typical arrangement of data samples in space. Raw data values 
(i.e.. Mg, M^, M2) for each resolution volume are the result of Poly Pulse 
Pair Processing echos from 50 transmitted pulses. 
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TABLE A.l 


Doppler Radar and Lidar Parameters 


Parameter 

NRO 

CIM 

Lidar 

Beamwidth (3 dB width) 

O 

• 

00 
♦ — 1 
o 

0.85° 

20 cm 

Resolution Volume Range Length 

150 or 450 m 

150 or 450 m 

320 m 

Peak Transmitted Power 

750 kW 

600 kW 

6 kW 

# of Transmitted Pulses for a 
Velocity Estimate 

64 

64 

50 

Pulse Width 

1 or 3 ys 

1 or 3 ys 

2 yS 

Unambiguous Range 

115 or 345 km 

115 or 345 km 

1364 km 

Nyquist Velocity 

34 . 3 or 
11.4 m s" 1 

35.4 or 
11.8 m s' 1 

39.8 m s“ 

Wavelength 

10.5 cm 

10.9 cm 

10.6 ym 

Maximum Range of Detected 
Targets (clear air) 

=80 km 

=80 km 

=10 km 

Maximum Range Sampled 

all 

all 

28.8 km 


a. INERTIAL NAVIGATION SYSTEM (INS) 

Two Litton LTN-51 units were used. These gyro-stabilized, four-gimbal , 
all-attitude units employ two degree-of -freedom gyros with gas bearings. The 
data furnished to the data acquisition system are: 

position (latitude, longitude, resolution 0.1 min), 

true heading (resolution 0.1 deg, accuracy 0.4 deg), 

track angle (resolution 0.1 deg, accuracy 0.5 deg); 

drift angle (resolution 0.1 deg, accuracy 0.5 deg), 

pitch angle (resolution 0.044 deg, accuracy 0.5 deg), 

roll angle (resolution 0.044 deg, accuracy 0.5 deg); 

groundspeed (resolution 0.51 m $ ■, accuracy 0.51 m s _1 ), 



wind speed (resolution 0.51 m s”V)i 
wind azimuth (resolution 1.0 deg), 


b. TOTAL AIR TEMPERATURE PROBE 

The Rosemount 102 AH2AG probe furnishes total air temperature and static 
air temperature data (resolution 0.1 deg) to an on-board computer which 
calculates airspeed (resolution 0.51 m s"*) in the heading direction, the 
airspeed data are furnished to the INS computer be used in calculating the INS 
windspeed and azimuth at the aircraft (see Fig. A. 3). 

C. DEW-FROST POINT HYGROMETER 

General Eastern 1011 thermoelectric hygrometer (response 1-3 deg per 
second). 

d. RADAR ALTIMETER 

Stewart -Warner APN-159, resolution 0.61 mi accuracy 1 percent. 

A. 4 STATIONARY AUTOMATED MESONET INSTRUMENTATION 

Although the lidar experiments were not within the time period of the 
major portion of the data collection effort by NSSL for 1981, four automated 
mesonet stations were operating. Below is a brief summary of the instrumenta- 
tion used to collect meteorological data near the surface (see Doviak (1981) 
for further details): 

a. WIND SPEED 

Specially calibrated NWS F-420-C rotating cup DC generator-type ane- 
mometer (resolution 0.4 m s" 1 ) 

b. WIND DIRECTION 

NWS F-420-C splayed-tail wind vane with modified direction transmitter 
(resolution 1.4 deg). 

c. DRY-BULB/WET-BULB TEMPERATURE 

Linearized Yellow Spring Model 44202 Thermistors, with self-wetting wet 
bulb, housed in standard NWS Stevenson screen (resolution 0.2 K). 
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d. STATION PRESSURE 


Texas Electronics Aneroid/LVOT Unit (resolution 0.5 mb), 
e. RAINFALL 

Belfort Model 5-780 weighing bucket raingage with weight-sensing load 
cell (resolution 0.6 mm). 

Dew point temperatures were computed from the dry-bulb and wet -bulb 
temperature data by integration of the Clausius-Clapeyron equation, assuming 
latent heat of vaporization to be constant. 



APPENDIX B 


ESTIMATION OF THE AVERAGE SURFACE HEAT FLUX OVER AN 
INHOMOGENEOUS TERRAIN FROM THE VERTICAL VELOCITY VARIANCE 

B.l INTRODUCTION 

It is straightforward to compute the surface fluxes of heat and momentum 
for the horizontally homogeneous atmospheric boundary layer (ABL) from the 
surface-layer profiles of wind and potential temperature by the profile method 
which uses the integrated form of the flux-profile relationships (Paulson, 
1970); it is also possible to compute both the friction velocity and surface 
heat flux from the surface-layer wind profiles alone, using the method sug- 
gested by Klug (1967). For the convective ABL, the surface heat flux can also 
be computed from just the surface-layer potential temperature profiles using 
the free convection formula, and this approach seems to give a heat flux esti- 
mate that is better than Klug's method and is as good as the profile method 
(Sundara-Rajan and Macklin, 1976). 

But it is not that easy to get a reliable estimate of the surface fluxes 

of heat and momentum if we don't have surface-layer profile data. This paper 

describes a method of determining the surface heat flux from the observations 

— 2 

of the vertical wind velocity variance (w 1 ) at the lower mid-levels of the 
convective ABL. This method is based on the assumption that the vertical wind 
velocity variance remains constant between the levels 0.2h and 0.6h where h is 
the height of the convective ABL - an assumption adequately supported by 
observations of the convective ABL over horizontally homogeneous sites (Kaimal 
et al., 1976). The heat flux estimates obtained with this method were sur- 
prisingly good, considering the fact that this method was used over a terrain 
which was far from being homogeneous. 

B.2 THE METHOD USED AND RESULTS 

Convective ABL observations over horizontally homogeneous sites (Kaimal 
et al., 1976) indicate, between levels 0.2h and 0.6h, 

« 0.39 (B.l) 

w* 
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1/3 

where w* is the convective scaling velocity; w* = ( 6hQ Q ) ' where Q 0 is the 
kinematic heat flux at the surface, and e = (g/0 ) is the buoyancy pa ram- 

'jyQ 

eter. We assumed the above relationship for w ,<: and computed the values for 
Q 0 from the observed values of the vertical velocity variance at mid-levels of 
the ABL, and the height (h) of the ABU 

The data we used were from the NSSL-KTVY tower located about 10 km north- 
northeast of Oklahoma City on a rolling, generally undeveloped terrain that is 
bounded on the south by developed urban terrain. A detailed description of 
the terrain and tower is given by Sanders and Weber (1970) and Lee and Stokes 
(1978). The highest level of measurement was at 444 m, and this restricted us 
to use the data collected during the early morning hours when the height of 
the ABL was still below 444 m. The height of the ABL was determined from the 
vertical profiles of virtual potential temperature, and only those runs for 
which the ABL height could be accurately determined were included in the 
analysis. 

The vertical velocity measurements were made at six levels (26, 45, 89, 
177, 266, 444 m) with R.M. Young Model 27103 Propeller anemometers at a 
sampling rate of 10 seconds. The somewhat slow response of the propeller 
anemometer can lead to an underestimation of vertical velocity variance in the 
surface-layer (Kaimal , 1975), but the lowest level of measurement used for the 
analysis reported here was 89 m. Also, since the free convection similarity 
implied in relation (1) holds good for highly convective conditions, only data 
for which h/L > -10 were used. The Obukhov length L(= -u* /'k.&Q ) for each run 
was computed from the values of Q 0 and the friction velocity (u*) determined 
from the mean wind speed values at 26 m and 45 m using the flux-profile 
relationships reported by Businger, et al. (1971) for the surface-layer. 
Although the u* estimates from the mean wind speed values at heights above the 
surface-layer should be considered relatively crude, they are sufficient for 
the computation of the values of the stability parameter h/L which are used 
here only for the selection of runs. The values of z, z/h, h/L, and Q 0 are 
given in Table B.l. 

We wish we had independent measurements of surface heat flux to test the 
accuracy of the Q 0 estimates using this method. We did not have any, but we 
estimated the surface heat flux indirectly from the observed boundary layer 
heating rate. The heating rate equation is given by (Tennekes, 1973): 
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(B.2) 


TABLE B.l, Surface kinematic heat flux estimates from the vertical 
velocity variance and the boundary layer heating rate. 


Qp (m s'k) 


DATE 

Time 

(LST) 

(m) 

z/h 

-h/L 

Variance 

method 

Heating 

rate 

method 

7-7-77 

8:15 

89 

0.34 

10.4 

0.104 

0.093 

7-10-77 

8:15 

89 

0.36 

17.63 

0.080 

0.089 

7-10-77 

8:45 

89 

0.34 

33.5 

0.084 

0.122 

7-14-77 

8:45 

89 

0.34 

24.5 

0.073 

0.082 

7-15-77 

8:45 

177 

0.40 

54.8 

0.093 

0.102 

7-18-77 

8:15 

89 

0.34 

11.5 

0.132 

0.119 

7-18-77 

8:45 

177 

0.40 

10.2 

0.152 

0.183 

7-19-77 

8:15 

89 

0.36 

49.1 

0,104 

0.112 

7-19-77 

8:45 

89 

0.29 

62.3 

0.080 

0.084 

7-25-77 

9:15 

177 

0.40 

23.4 

0,184 

0.163 

7-28-77 

9:15 

177 

0.50 

47.3 

0.121 

0.147 

7-30-77 

8:45 

89 

0.44 

61,3 

0.149 

0.100 

8-26-77 

8:45 

89 

0.20 

44.9 

0.260 

0.305 

9-3-77 

8:15 

89 

0.58 

12.23 

0.118 

0.148 

9-3-77 

8.45 

89 

0.39 

21.0 

0.118 

0.120 

5-11-79 

7:15 

89 

0.37 

80.3 

0.119 

0.097 

5-17-79 

8:15 

177 

0.50 

16.2 

0.154 

0.149 

5-19-79 

8:15 

177 

0.42 

10.6 

0.149 

0.141 

6-3-79 

10:15 

177 

0.43 

577.5 

0.332 

0.274 

6-4-79 

7:15 

89 

0.60 

12.6 

0.073 

0.077 

6-4-79 

7:45 

89 

0.49 

54.8 

0.094 

0.087 

6-4-79 

8:15 

89 

0.34 

11.2 

0.113 

0.126 

6-9-79 

10:15 

177 

0.44 

10.4 

0.275 

0.296 

6-10-79 

7:45 

89 

0.27 

36.5 

0.106 

0.123 

6-10-79 

8:15 

89 

0.27 

10.2 

0.127 

0.144 
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It was found that for all the runs included here, the strength of inversion 
was negligible, and hence the downward kinematic heat flux at the inversion 
base could be neglected. 

The Q 0 values computed using the variance method and the boundary layer 
heating rate method are compared in Fig. B.l. The agreement between the two 
sets of values is very good, the standard error being 0.024 m s"*K which is 
almost the same order as that for the profile method (Sundara-Rajan and 
Macklin, 1976). 

B.3 CONCLUDING REMARKS 

It seems that this method of computing the surface kinematic heat flux 
from the vertical velocity variance gives a fairly reliable estimate of the 
average surface heat flux values over an inhomogeneous terrain. Unfortu- 
nately, independent and direct estimates of surface heat flux were not avail- 
able to test rigorously the validity of this method. What was done was a com- 
parison with the surface heat flux estimates from the boundary layer heating 
rates, and the results are encouraging. Further documentation from other 
inhomogeneous sites is required to confirm the results of this study. 



Figure B.l Comparison of surface 
heat flux estimates from the 
variance method and the heating rate 
method. 
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APPENDIX C 

INSPECTION AND VERIFICATION OF CV990 INS DATA 


In this Appendix the process of verifying the aircraft position and atti- 
tude data, and correcting any errors found in those data is described. The 
aircraft position was defined in terms of latitude and longitude recorded at 
the time of each frame in each run. However, inspection of the raw data indi- 
cated that these latitude and longitude measurements were either not being 
updated frequently enough or were subject to some sort of round-off error, as 
the values tended to remain constant for periods as long as 6-8 s at a time. 
Although it was clear that many of the specific positions suggested by these 
latitude and longitude values were in error, it was also evident that if the 
latitude and longitude measurements were averaged with a "running mean" filter 
extending about 10 data points either side of each given datum in the respec- 
tive series of latitudes and longitudes, the resulting set of aircraft posi- 
tions seemed to give a good representation of the aircraft flight path. 

This was verified by comparing the smoothed INS-derived positions with 
those obtained from nadir photographs taken at about 10 s intervals during 
several of the lidar data runs. 

Although the INS-derived attitude and velocity data contained significant 
error, they showed sufficient internal self-consistency that they could be 
used to generate a first guess of aircraft positions, which would later be 
compared with photo-derived positions to yield estimates of the larger-scale 
INS biases in the lidar velocities. Thus it was deemed satisfactory for the 
overall analysis to obtain the aircraft latitudes and longitudes by integra- 
tion forward in time using the starting positions of the aircraft at the 
beginning of each run as reference points and the INS-measured groundspeeds 
and directions as the velocity. This procedure would also permit a verifica- 
tion to be made of the interrelationships (Fig. C.l) among the airspeeds, 
groundspeeds, heading angles and drift angles reported on tape. It was 
especially important to verify these data since they would have to be used to 
compute the track angle which was not furnished on the data tape. In addi- 
tion, verification of drift angle values was crucially important due to the 
sensitivity of the lidar-derived wind vectors to drift angle fluctuations. 
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Figure C.l Aircraft flight 
parameters relevant to data 
analysis. 


C.l ERRORS IN CONVERSION FACTORS 

Prior to discussing this verification process, one obvious error in the 
INS derived wind azimuth data (INS derived wind azimuth describes the azimuth 
of the wind vector estimated by the aircraft INS and is recorded at the time 
of each lidar frame) had to be corrected before beginning the verifications. 
Examination of the INS derived wind azimuth values showed some angles larger 
than 360 °, which was easily traced to the use of an incorrect conversion 
factor during the creation of the data tape. Multiplication by a suitable 
ratio of factors remedied this problem and gave believable INS wind direction 
values. 

In the first step of checking the aircraft attitude and position data, it 
became apparent that forward integration of the aircraft velocity vector, as 
defined by the groundspeed and track angle (= true heading + drift angle), 
yielded discrepancies in aircraft position which increased steadily with 
time. Since the INS wind azimuths were previously found to be in error by a 
factor given by the ratio of two conversion factors, it was suspected that a 
similar problem might be contaminating the drift angle values. Several test 
solutions of the vector triangle relating aircraft ground-relative motion, 
air-relative motion and INS-measured winds quickly confirmed this 
hypothesis. Although the INS wind azimuths had been found to be in error. 
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confidence in the INS wind speeds was warranted because of their agreement 
with the magnitudes suggested from the soundings taken on the day of the 
experiment. This impression was confirmed during subsequent analysis of the 
radar and lidar wind data (see section 1.6). In fact, it was ultimately 
ascertained that the errors in INS wind azimuth and drift angle were comple- 
mentary: the conversion factors for each quantity had been inadvertently 

switched during the writing of the tape data. It was thus a simple matter to 
recompute the correct values of the aircraft drift angles by multiplying them 
by the ratio of the correct conversion factor to the incorrect one, as had 
been done for the INS wind directions. 

Late in the analysis stage a similar conversion factor problem was 
detected in the recorded values of roll angle, which were all too small by a 
factor of 2.27. Evidence of the roll angle errors came from a sequence of 13 
ground returns observed in run 4 during a left turn in which lidar data were 
collected; confirmation of the errors was obtained by comparison of the roll 
angle values recorded on the lidar data tape with those recorded independently 
by the aircraft INS system. Because of the lateness of the discovery of the 
roll angle error and because the roll angle errors had no impact on the calcu- 
lations of aircraft position, no correction was applied to erroneous roll 
angle values. 

C.2 CORRECTION TO THE RECORDED INS DERIVED POSITION DATA 

After modifying values of aircraft drift angle to reflect use of the 
proper conversion factor, it was possible to compute consistent track angles 
and to carry out the time integration of the ground-relative velocity vector 
in order to establish a consistent set of latitudes and longitudes. The 
actual integration was performed in a local x-y coordinate system having an 
origin at the starting point of each run. The calculations were made using 
double precision arithmetic, as the use of single precision arithmetic was 
found to cause position biases large enough to interfere with accurate assess- 
ment of the INS position errors. The revised computed aircraft x and y posi- 
tion coordinates were then transformed back to latitudes and longitudes by 
means of: 

LAT(I) = LAT(I-l) + A*(X( I ) - X(I-l) ) (C.la) 

LON(I) = LON(I-l) + A*(X( I) - X(I-l) ) / (COS (LAT { I ) ) (C.lb) 

where "A" is a factor which converts distance to degrees. 
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Latitude-longitude values obtained in this manner followed quite closely 
the track suggested by the smoothed version of the data given on the tape, and 
since these computed positions were obtained in a deterministic manner, they 
were deemed appropriate for use in all subsequent calculations. 

A sample of the uncorrected tape data from the first 62 frames of run 2 
dealing with aircraft position and attitude is given in Table C.l. The posi- 
tion problems discussed above are evident in the latitudes ("LAT") and longi- 
tudes ("LON") as well as the calculated x and y coordinates of the aircraft 
relative to the start of the run ("XCV990" and "YCV990" ). Note in Table C.l 
the tendency of all of these quantities to remain constant for periods of time 
lasting up to 4-5 s. A typical example is seen in the sequence of XCV990 
values from frames 36 through 39. An even more disturbing sequence is seen in 
frames 59-61, where both the XCV990 and YCV990 values, as well as LAT and LON, 
remain constant, suggesting that the aircraft remained stationary for several 
seconds! It should be noted that although the drift angle data (DRANGL) 
appear to be reasonable, they still contain the conversion factor error 
described above, and are too small by more than a factor of two. 

Another portion of the output of tape data is given in Table C.2; this 
table again deals with the first 62 frames from lidar run 2. Of particular 
interest are the values of INS wind azimuth (WINDAZ), some of which (frames 
1-3 for example) are seen to exceed 360°. 

The effect of the corrections made during the data review stage of the 
analysis can be gauged by examining Tables C.3 and C.4, which are revised 
versions of Tables C.l and C.2, with some format changes in the case of Table 
C.4 to accommodate some new calculated parameters. Note the increased smooth- 
ness of rate of change of the LAT, LON, XCV990 (renamed "C1X990" in Tables 
C.3-C.4) and YCV990 (renamed "C1Y990") arrays in Table C.3, and the conversion 
of the erroneous WINDAZ values in Table C.2 to credible values of "AZINS" in 
Table C.4. Note also that in Table C.4 the values labelled "VINS" represent 
INS-derived wind magnitudes in meters per second obtained by converting the 
"WINDSP" values in Table C.2 from knots. In Table C.3 values of drift angle 
can also be seen to have increased, some hint of problems to be discussed 
later is seen in the fact that the drift angles are not resolved to a preci- 
sion beyond 0.1 deg. Note also that a number of other parameters have been 
revised as well; the reasons for these changes and the method of making the 
changes will be discussed below. 



FRAME TIME 

LAT LON 

ZPHiiS 

ZR ADR 

XCV990 

Y CV 99 0 

TRUHEAO 

ORAN OLE 

1 

201810. 

35. 113 -97.737 

0. 

2424 . 

0. 

0 . 

310.999 

0 .703 

2 

201811. 

35.113 -97.737 

0. 

2404. 

0. 

0. 

311. 262 

0.57 1 

3 

201812. 

35. 1 15 -97.740 

0. 

24 10 • 

-200. 

180. 

31 1 .262 

0.571 

4 

201815. 

35.117 -97. 740 

0. 

2470 . 

-280. 

360. 

312.185 

0.439 

5 

201617. 

35.118 -97.745 

0. 

2394 . 

-740. 

560 . 

313.503 

0.044 

6 

201616. 

35.118 -97.747 

0. 

2376. 

-900. 

560 . 

313. 767 

0.044 

7 

201619. 

35.122 -97.748 

0. 

2390 . 

-1 0-4 0 • 

920. 

313.767 

0.132 

8 

20182 1 . 

35.122 -97.740 

0. 

24 24 . 

-1040. 

920. 

314.690 

0.396 

9 

20 1622. 

35.123 -97.740 

0. 

2402. 

-1 040. 

1 1 00 . 

314 . 690 

0 .3 0 3 

10 

20 1824. 

35.123 -97.752 

0. 

2386. 

-1340. 

1100. 

315.569 

0.3 0 6 

1 1 

201825. 

35.127 -97.755 

0. 

2414. 

-1640. 

1480 . 

316.272 

0.3 06 

12 

201827. 

35.120 -97.757 

0. 

2398. 

-1800. 

1660. 

316.799 

0.439 

13 

201828. 

35.132 -97.750 

0. 

2400. 

-1960. 

204 0. 

316.799 

0 .439 

14 

201830 . 

35.132 -97.758 

0. 

2424. 

-i960. 

204 0. 

317.371 

0.396 

15 

201831. 

35. 132 -97.750 

0. 

2456. 

-1960. 

2040. 

317.678 

0.483 

16 

20 1833. 

35.133 -97.762 

0. 

2394. 

-2260. 

22 20. 

317.986 

0.5 7 1 

17 

20 1834. 

35.135 -97.762 

0. 

2 372 . 

-2260. 

2400. 

318.301 

0.571 

18 

201836. 

35. 137 -97.767 

0. 

2338. 

-2720. 

2500. 

318.689 

0.527 

19 

20 1 837 . 

35.138 -97.768 

0. 

2348. 

-2060. 

2780. 

318.609 

0.5 27 

20 

201839. 

35.130 -97.760 

0. 

2334 . 

-2860. 

2780. 

318.689 

0.43 9 

21 

201640. 

35. 130 -97.768 

0. 

2354. 

-2860. 

2780. 

319. 568 

0.220 

22 

20 1842. 

35.143 -97.772 

0. 

2398. 

-3160. 

33 20. 

319. 568 

0 .2 2 0 

23 

201843. 

35.145 -97.772 

0. 

2402. 

-3 160. 

35 20 . 

319. 56 e 

0.220 

24 

20 1845. 

35.145 -97.775 

0. 

2374. 

-3460. 

3520. 

319.568 

0.264 

25 

201846. 

35. 148 -97.770 

0. 

2352. 

-3780. 

3800. 

319.963 

0.3 52 

26 

201847. 

35.140 -97.778 

0. 

2362. 

-3700. 

3080. 

320. 402 

0.0 

27 

20184 9 . 

35.140 -97.780 

0. 

2322. 

-3920. 

3880. 

320. 402 

0.0 

28 

201050. 

35.148 -97.700 

0. 

2338. 

-3 92 0. 

3880 • 

320 . 403 

0.0 

29 

201851 • 

35.152 -97.782 

0 • 

2352. 

-4 080. 

4260. 

320.402 

0.132 

30 

201852. 

35.152 -97.702 

0. 

23 70. 

-4080. 

4 2 60 • 

320.491 

0.132 

31 

20 1 654. 

35.157 -97.703 

0. 

2356. 

-4220. 

48 00. 

320.579 

0.176 

32 

201855. 

35.157 -97.703 

0. 

2372 . 

-4220. 

4 8 00. 

320.579 

0. 1 76 

33 

201856. 

35.150 -97.785 

0. 

2368. 

-4380. 

5000 • 

320.579 

0.0 80 

34 

201858. 

35. 158 -97.785 

0. 

2392. 

-4360. 

50 00 . 

320. 579 

0 .0 88 

35 

20 1859. 

35.158 -97.705 

0. 

24 52. 

-4380. 

50 0 0 . 

320.796 

0.060 

36 

20 1900. 

35.162 -97.792 

0 . 

2426. 

-4980. 

53 6 0 . 

320.798 

0.264 

37 

20190 1 . 

35.162 -97.792 

0. 

2444. 

-4980. 

5360. 

321 .194 

0 .264 

36 

201903. 

35.163 -97.792 

0. 

2390 . 

-4 980. 

5540 . 

321 . 194 

0.308 

39 

201904. 

35.163 -97.792 

0. 

2382. 

-4980. 

554 0. 

321 .370 

0 .220 

4 0 

201905. 

35.167 -97.795 

0. 

2372. 

-5300. 

5920. 

321.501 

0. 1 32 

4 1 

201906 . 

35.160 -97.795 

0. 

23 70 . 

-5300. 

6100. 

321.501 

0. 1 32 

42 

201908. 

35.168 -97.795 

0. 

2378. 

-5300. 

6100. 

321 .501 

0.044 

43 

201909. 

35.170 -97.798 

0. 

2364. 

-5600. 

62 80. 

321 . 50 1 

0.176 

44 

201910. 

35.172 -97.798 

0. 

238 4. 

-5600. 

64 80. 

321.677 

0.132 

45 

201912. 

35.172 -97.002 

0 • 

2376. 

-5900. 

64 80. 

321 . 677 

0 . 1 32 

46 

201913. 

35.172 -97.802 

0. 

23 78. 

-5900. 

64Q0. 

321 .677 

0.176 

47 

20191 4. 

35.172 -97.802 

0. 

2404. 

-5900. 

6480 . 

322.600 

0.176 

4 8 

20191 5. 

35.178 -97.805 

0. 

24 18. 

-6200. 

7220 . 

322.600 

0.1 76 

49 

201917. 

35.170 -97.805 

0. 

2 492. 

-6200. 

7220. 

323.39 1 

0.176 

50 

201918. 

35.180 -97.807 

0. 

2462 • 

-6360. 

74 00 . 

323.39 1 

-0.264 

51 

20 1919. 

35.182 -97.800 

0. 

2454. 

-6500. 

75 80. 

323.171 

0.132 

52 

201921 . 

35.183 -97.612 

0 • 

2462 . 

-6800. 

7760. 

323. 1 71 

0. 1 32 

53 

201922. 

35.183 -97.812 

0. 

24 12. 

-6800. 

7760. 

323.075 

-0.088 

54 

201923. 

35. 185 -97.613 

0. 

2404. 

-6960. 

79 60 . 

323.567 

0.0 

55 

201924. 

35.185 -97.013 

0. 

2422. 

-6960. 

7960. 

323.567 

0 .0 

56 

20 1926. 

35.185 -97.613 

0. 

2436. 

-6960. 

7960 • 

323. 075 

0.088 

57 

201927. 

35.190 -97.813 

0. 

2400. 

-6960. 

8.5 00. 

324. 182 

0.088 

53 

201928. 

35. 192 -97.017 

0. 

2446. 

-7260. 

8700. 

324.490 

-0.088 

59 

201930. 

35. 193 -97.010 

0. 

2430 . 

-7420. 

8880. 

324.490 

-0.008 

60 

201931 . 

35. 193 -97.810 

0. 

24 22 . 

-7420. 

06 80 . 

324 .490 

-0.068 

6 1 

201932. 

35.193 -97.810 

0. 

24 36. 

-7420. 

8880 . 

324 .402 

0.132 

62 

201933. 

35.197 -97.823 

0 . 

24 1 6 • 

-7800. 

9240. 

324 .402 

0.1 32 


Table C.l Sample of raw position data NASA CV 990. 
See text for definition of terms. 
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WINOSP 

WINOAZ - TTOTC ., 

T STAT Oil WPT _ 

TSRFC 

TSEN 5 

T 1 MP OS 

kr£qlo 

UF K 50 T 

- 11.00 

3 75 .8 0 . 

,33 * 90 _ 

.23 .9 Q — 10.90 

31 .90 

0.10 

0.0 

0 . 0 

0 . 0 

ji.oo 

375.50 

33.90 

23.90 19.10 

31.20 

0 . 1 0 

0.0 

0.0 

0.0 

i i .oo 

375 .50 

33.90 

23.90 10.60 

32.00 

0 . 10 

0 . 0 

0.0 

0.0 

.. 10*00 

318 . 60 - 

. 33 . 60 .. 

23 . 60 . 20, 0 0 . 

29. 70 

0.10 

0.0 

0 .0 

0.0 

- 10.00 

31 9 . 6 0 - 

. 33 . 70 , 

. 23.0 0 _ 18 . 00 . 

.29 • 90 . 

— 0.10 

0.0 

0.0 

0. 0 

10.00 

355 .00 

33.60 

23.80 19.40 

30.4 0 

0.10 

0. 0 

0.0 

0.0 

1 0 . 00 

355.00 

33. 70 

23.70 20.30 

30.40 

0.10 

0.0 

0.0 

0.0 

10 . 00 . 

,3 55 . 0 0 , 

. 33.70 

23 . 70 _ 19 . AO - 

.3 0.10 

- 0.10 

0 . 0 

0 . 0 

0 . 0 

_ i 3. 00 

- 357 . 30 . 

. 33. 70 

23 . 70 ^ 19. 50 

30.70 

- 0.10 

0 . 0 

0 .0 

0. 0 

1 J . 00 

373.20 

33.60 

23.50 , 19.80 

31.80 

0 . 10 

0. 0 

0. 0 

0. 0 

14.00 

361 #80 

33 .70 

23.50 19.30 

31.70 

0 . 10 

0.0 

0.0 

0. 0 

... i J .00 

. 366 .4 0 

33 .60 

23 . 50 — 1 9 . 20 - 

3 1 .0 0 . 

— 0,10 

- 0.0 

0-0 

0 . 0 

.-. 14 . 00 - 

_ J 57 .3 0 

33.00 

. 23 . 50 — 20 . 50 .. 

3 1 .30 

0.10 

0. 0 

0.0 

0.0 

14.00 

357.30 

33.80 

23 , 70 . 19.30 

30. 5 0 

0.10 

0.0 

0.0 

0. 0 

14.00 

368.60 

33.90 

23.70 1 0.00 

29.50 

0.10 

0.0 

o.o . 

0.0 

14 .00 

. 368.60 

33.90 

23 , 70 . 1 9.60 

3 0.40 

0.10 

0.0 

o.o 

0 . 0 

- 15 . 00 . 

- 366.40 

33 . 90 . 

- 23 . 70 - 19.00 

29.90 

. 0.10 

0 • 0 

0.0 

0.0 

14.00 

368 .60 

34.00 

23 * 60 , 19.90 

3 0.70 

0.10 

0 . 0 

0.0 

0. 0 

14.00 

368.60 

34. 00 

23.60 19.50 

3 0.40 

0.10 

0.0 

0 . 0 

0 . 0 

14 . 00 . 

- 36 6 .60 

34 . 10 . 

- 23 . 60 — 1 9.80 

. 31.10 

0 . 10 

0.0 

o.o 

0. 0 

-1 1 .00 

„ 3 39 • 1 0 

34.10 

23 . 70 - 20.90 

. 30.50 

0.10 

0.0 

o.o 

0.0 

11.00 

339. 1 0 

34.20 

23.60 18.80 

29.90 

0,10 

0.0 

0.0 

0.0 

11.00 

366.40 

34.20 

23.60 19.50 

30.50 

0.10 

0.0 

0.0 

0. 0 

.. 10 . 00 

355 . 00 - 

. 34 .00 

23 . 70 — 19.00 

30.50 

-. 0.10 

0.0 

0.0 

0. 0 

- 11.00 

316.30 

34.00 

23 . 60 - 20.60 

29.40 

, 0.10 

0. 0 

0.0 

0.0 

i l . 00 

316 .30 

33.90 

23.50 19.70 

3 0.40 

0.10 

0.0 

0 .0 

0.0 

la . oo 

316.30 

33.90 

23.60 20.90 

30.40 

0.1 0 

0.0 

0.0 

o.o 

- 12,00 

325.40 

33.90 

23.60 19.80 

31.00 

0. 10 

0.0 

0 . 0 

0. 0 

. 13.00 

343-6 0 . 

, 33.90 

23.70 19.00 

29.90 

0 . 10 

0.0 

0 .0 

0. 0 

l J .00 

343 .60 

33.80 

23.70 20.10 

29.00 

0 . 10 

0.0 

0.0 

0.0 

14.00 

343 .60 

33.80 

23.50 19.50 

30.70 

0 . 10 

0. 0 

0 .0 

0.0 

- 14.00 

343 . 60 - 

. 33.60 

23 . 60 - 20 .10 

. 30.00 

- 0.10 

0.0 

0 .0 

0. 0 

- 14.00 

332 . 20 . 

. 33.8 0 

23 * 60 - 21 . 20 - 

30.60 

. 0.10 

0.0 

0.0 

0.0 

14.00 

34 1 .3 0 

33.90 

23.60 19.90 

29.80 

0.10 

0.0 

0.0 

0.0 

14.00 

350.40 

33.90 

23.60 20.00 

29.20 

0 . 10 

0.0 

0.0 

0. 0 

.- 14.00 

350 . 40 . 

- 33.90 

23 . 70 , 21 . 50 , 

. 29 . 20 , 

- 0.10 

0.0 

0.0 

0.0 

-. 14.00 

350.40 

34.00 

23 . 60 . 19.60 

29.90 

0,10 

0.0 

0.0 

0. 0 

14.00 

352.70 

34 . 00 

23.60 19.70 

2 9.80 

0.10 

0.0 

0.0 

0.0 

14.00 

352.70 

34 . 00 

23.70 18.50 

3 0.00 

0.10 

0.0 

0.0 

0.0 

13.00 

334 . 50 . 

34 .00 

23 . 70 . 10.80 

3 0.00 

0.10 

0.0 

0.0 

0.0 

13 . 00 

334 .50 

34,00 

23 • 70 _ 20.30 

30.80 

0.10 

0. 0 

0.0 

0. 0 

13.00 

334 .5 0 

33.90 

23.50 19.50 

29.20 

0.10 

0.0 

0.0 

0. 0 

13 . 00 

334.50 

34 . 1 0 

23.50 19.40 

29 . 10 

0.10 

0.0 

0.0 

0. 0 

_ 12.00 

. 345 . 90 - 

. 34 . 00 , 

- 23 . 30-2 0 . 80 . 

29.90 

- 0.10 

. .. 0.0 

0.0 

0.0 

- 12 . 00 . 

- 345.90 

34.10 

23 . 60 - 18.80 

3 0.90 

0 . 10 

0.0 

0.0 

0 , 0 

12.00 

345.90 

34.10 

23,50 19.80 

30.20 

0.10 

0. 0 

0.0 

0.0 

' 13. 00 

304.90 

34.10 

23.70 21.10 

31.00 

0 . 10 

0.0 

0 .0 

0.0 

— 13. 00 

304 . 90 . 

- 34.00 

23 . 60 — 1 9*50 

31 .20 

0.10 

0.0 

0.0 

0.0 

__ 1 3 • 00 

_ 3 04 .90 

34.30 

23 * 50__2 0 . 20 

29.40 

0.10 

0.0 

0 .0 

0.0 

13.00 

3 04 .90 

34.30 

23.70 19.90 

29.90 

0.10 

0 , 0 

0.0 

0.0 

12.00 

336.80 

34.40 

23.50 19.00 

29.00 

0. 10 

0.0 

0.0 

0.0 

.. 10.00 

310.30 

_ 34 . 60 

23.4 0 1 9.80 

. 30.80 

0 . 10 

0.0 

0.0 

0.0 

- 10.00 

316.30 

34.60 

23.60 19.30 

31.10 

~ 0 . 10 

0.0 

0.0 

0.0 

H . 00 

336.60 

34.70 

23.50 19.80 

31 .30 

0 . 1 0 

0 . 0 

0.0 

0.0 

8 . 00 

336.00 

34.60 

23.50 20.40 

30.90 

0.10 

0.0 

0.0 

0.0 

0. 00 

336.00 

34.50 

23.50 20.50 

30.40 

0.10 

0.0 

0.0 

0.0 

- 8 . 00 

“ 336.00 

34 .40 

23 . 30 ^ 20.40 

2 9.10 

0.10 

0 . 0 

0.0 

0.0 

0.00 

269 .00 

34.30 

" 23 * 10 " 20.00 

29.4 0 

0.10 

0.0 

0.0 

0. 0 

8.00 

289.00 

34.40 

23.20 20 V . 60 

30.20 

0. 10 

0.0 

0 . 0 * 

0.0 

8. 00 

209.00 

34 . 1-0 

23 . JO 19.50 

32.40 

0. 10 

o.o 

0.0 

0. 0 

* . 0 0 

345.90 

34 . 30 

2 J . 20 _ 20.20 

34. 10 

0.10 

0.0 

0.0 

0. 0 

«) . 0 0 

3 1 y .60 

34 . 30 

23.30 20.50 

34 .00 

0.10 

0.0 

0 .0 

0.0 


Table C 2. Sample of raw meteorological data, NASA C V 990. 
See text for definition of terms. 
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OltSSflAL PAGE IS 

OF POOR QUAUTY 


FRAME TIME 
— 1 201810. 
2 201011 . 

3 201812. 

4 201015. 
_ 5_ 201817. 

6 201010 . 

7 201 01 0. 

8 201021 . 
-9 201822. 

10 201823. 

11 201025. 

12 201827. 
- 1 3—2 01828. 

14 201030. 

15 201031. 

16 201833. 
- 17-201034. 

18 201836. 

19 201837. 

20 201839, 
-21-2 0 1840. 

22 201842. 

23 201843. 

24 201845. 
-25—20 1846* 

26 201847. 

27 201849. 

28 201850. 
^29-20-1851. 

30 201852. 

31 201854. 

32 201855. 
>23-2018 56, 

34 201857. 

35 201859. 

36 20 1500. 
-37-20.1.901 , 

38 20 1902. 

39 201904. 

40 201905. 
-4 1.-20 1906. 

42 201908. 

43 201909. 

44 201910. 
-4 5.-2019.1 1. 

46 201913. 

47 201914. 

48 201915. 
-49-201917. 

50 201918. 

51 201519. 

52 201920. 
-53—20 1922. 

54 201923. 

55 201924. 

56 201926. 

57 2 01927. 

58 201928. 

59 201929. 

60 201931. 
-6 1 20 19 32. 

6 2 2019 33. 


TS 
-0 
1 
2 
5 
7 
0 
9 
1 1 
12 

14 

15 
17 

-18 
20 
21 
2 3 
-24 
26 
27 
29 
5—30 

0 32 

5 33 

1 35 
4 — 3 6 

7 37 

0 3 9 

2 40 

6 4 1 

8 4 2 

1 44 

4 4 5 

7- 4 6 

9 48 

2 49 

6 50 

.7 — -SI 
9 53 

2 54 

5 55 

7 56 

1 58 

3 59 

6 60 

9 62 

2 63 

5 64 

t>5 
67 


EC CL AT CLON ZHADfi 
.0 -35. 1 13 -- 97. 737- 739. 
.40 35.115 -97.738 733. 

.80 15.116 -97.740 735. 

.60 35. 1 10 -97.743 753. 

. 00 3 5.-1 19--97. 74 5—730. 
.40 35. 121 -97.746 725. 

.80 35.122 -97.748 728. 

. 20 3 5. 123 -97. 749 739. 

.60 3 5. 124—97.751 —732. 
.00 35.126 -97.752 727. 

. 40 35. 127 -97. 754 736, 

.00 35.129 -97.756 731. 

, 50-35. 1 30 — 97,757 — 732. 
.00 35. 1 31 -97.759 739. 

. 50 35. 1 33 -97. 760 749. 

. 00 3 5. 1 34 -97. 76 2 7 JO. 

. 50 3 5. 1 36—97.-764 — 723. 
.00 35, 137 -97.765 713. 

.50 35,139 -97.767 716. 

• 00 35. 140 -97.768 7 1 1. 

. 50 35. 14 2—97-7 7 0— 718. 
.00 35.143 -97.771 731. 

. 50 3 5. 145 -97. 773 732. 

.14 35.146 -97.774 724. 

. 4 3-3 5^1.48—97.-7-76 7 17.. 

.71 35.149 -97.777 720. 

.00 35.150 -97.778 708. 

. 29 3 5. 152 -97.700 7 13. 

.57-^-5.453— 37,-78 t 7.17, 

. 86 35.154 -97.782 724. 

.14 35.156 -97.784 
. 43 3 5. 157 -97.785 
. 7 I -3 5. 45 6-^.9 7. J7 86- 
.00 35. 160 -97.788 
. 29 35. 161 -97.789 
.57 35.162 -97.790 
. 86-35.4 63 — =57 .-7 9-1. 


C 1.X 99.0 C 1 Y 0 5 0 THUH EAD 


C 

-14 6. 
-293. 
-58 4. 
- --72 9. 

-87 3. 
-101 7. 
-1159. 
-1299. 
-143 9. 
-1578. 
-1734. 
-- 188 0. 
-2025. 
-2169. 
-2312. 
-2454. 
-2595. 
-27 36. 
-287 0. 


718. 

723, 

722.- 

729. 

74 7. 

739. 

7 45.- 

728. 

726. 

72 3. 

722,. 

725. 


310. 999 

311. 26 2 

311. 262 

540. 3 1 2. 185 
676. 313,503 
3 14, 3 1 3. 76 7 
""" 3 13. 7b7 

314,690 
314. 65 0 
1381. 315. 569 
1528. 316. 272 
16 98. 3 1 6. 799 
18 58. 3 1 6. 79 5 
2019. 317.371 
2182. 317.673 
23 47. 31 7. 98 C 
2513. 318. 331 
26 80. 318, 639 
2848. 318.689 
30 15. 318. 638 
3019. —31 83.-319. 568 
-316 0. 33 51. 319. 568 

-3300. 35 19. 31 9. 568 

-3454. 3703. .31 9. 568 

-3574. .-..38 4 fl, 319. 563 
-3693. 3993. 320.403 

41 38. 32 0. 4 03 

4283. 320.403 
4428._320. 403 
4574. 320. 491 

4721. 320.5/9 
4867. ““ 

5013 


0. 
135. 
269. 


9 53. 
1094. 
12 37. 


.00 35.165 -97.793 
. 29 3 5. 166 -97.794 
.57 35.167 -97.795 
.86 .-3 S.J 69-^97,796. 

. 14 35. 170 -97.798 

.43 35. 171 -97.799 721. 

.71 35.173 -97.800 727. 

. 0 Q_3 5, 1.7 4— =.97 ,8 (LI 724, 

. 29 3 5. 175 -97.003 725. 

.57 35.177 -97.804 733. 

.86 35.178 -97.805 737. 

. 1 4 -3 5 .>179—- 9 7 ,'80 6_7 6 0 . 
3 68.43 35.181 -97.808 757. 

6 69.71 35. 182 -97. 009 748. 

9 71.00 3 5. 103 -97. 810 750. 

2 — 72 -2 9— 3-5.-1 85-^9 7,8 33 735, 

5 73. 57 35. 106 -97.313 733. 

7 74. 86 35. 103 -97.814 738. 

1 76.14 35.189 -97.015 743. 

3 -77. 43-35. 190 .-97.016 - 756, 

6 78.71 35.192 -97.818 746. 

9 80.00 3 5. 193 -97.819 74 1. 

2 01.29 35.194 - 97.820 738 . 

5 — 62 , 5 7-3 5.-1 96 -- 97,821 — 7 43 . 

7 HJ. 86 35. 197 -97. 822 736, 


-301 3 _ 
-3933. 
-=4052,. 
-4172. 
-4291. 
-4409. 

-=4528- 
-464 8. 
-4767. 
-4085. 
>.-5003, 
-5107. 
-522 4. 
-534 0. 
.-54 57, 
-5575. 
-5692. 
-5809. 
>-59 2 5 , 
-6041. 
-6157. 
-6270. 
>-6383, 
-6497. 
-661 1. 


320. 579 
32 C. 57S 
5160. 320.579 
5306. 320. 798 
54 54. 320, 798 
-5.6.0.2.-32.1,1.9 4. 
5734. 32 1. 194 
5883. 321. 370 
6032. 321. 501 
-6_1 8.0,-321,501. 
6329. 321. 501 
6477. 32 1.501 
6626. 321,677 
_6775._32 1, 677, 
6925. 321. 677, 
7075. 322.600 
7226. 322.600 
.7378. 323. 391 
7530. 323.391 
7682. 323. 171 
7835. 323. 171 


-6724. 

6 83-7 -79 8 8 ._3 2 3 ,8 7 5. 

-695 0. 81 41. 32 3. 567 

'8294* 323. 567 
8448. ““ 


-7063. 
-717 5. 
-7287.> 
-739 3. 
-7510. 
-7622. 
-- 773 3.- 
-7844. 


. 323.875 
.8603, .3 24.182 
8758. 324.495 
89 14. 324. 4*90 
9069. 324.490 
-9225,-324. 402 
9332. 324. 402 


Dfi ANGLE 
1,6 00 
1. 300 
1.300 
1.0 00 

n. loo 

o. no 

0.300 
9.9 90 
0.700 
0.700 
0.7 00 
1.000 
1.000 
0.9 00 
1.100 

1.3 0 r 

1 . 3 09 
1.200 
1 .209 
1.0 00 
0.500 
0.509-- 
9. 500 
0.6 00 
0.800 
0.0 
0.0 
0.0 

. 0.300- 
0.3 09 
0.4Q0 
0.490 
0 . 200 - 
0.200 
0.200 
0,600 
- 0.6 00 _ 
0.700 
0.500 
0.300 
- 0 . J00 

o. mo 

0.409 

0.300 

0 .300 

0.400 
0.4 00 
0.400 
> 0.400- 
-0.600 
0,300 
0.300 

^Q.2aau 
0.0 
0.0 
0,200 
G. 200 
- 0.200 
- 0.200 
- 0.200 
— C .300 
0.300 


Table C 3. Sample of corrected position data NASA CV-990. 
See text for definition of terms. 
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AZ990 

fi A N 9 9 0 

EPS990 

IFCfie" A2CCS 

VINS 

AZINS 

VR INS 

TTCT 

SAT* 

TCEK 

TSFC~ 

170.2 

40755. 

C. 65 l 

1 

240 

80 

5.7 

1 65.0 

-5.6 

33 .9 

23 .9 

18 .9 

2 l .9 

170.2 

40599 . 

0.646 

- 1 

2 C 2 

10 

5.7 

165.0 

-5.6 

22.5 

22 .5 

15.1 

21 .2 

170.5 

40442 • 

0.840 

1 

24 i 

10 

5.7 

1 65 .0 

-5 .6 

32 .9 

22 .5 

16.6 

2 2 .C 

170.9 

4 0 1 2 C . 

C .6e2 


242 

oc 

5. 1 

140.0 

— 4.4 

33 .6 

23 .6 

20 .0 

29 .7 

171.0 

39973 . 

C .862 


2 C 4 

3 C 

5 . 1 

1 4 C . 0 

-4.4 

22. 7 

2 3*8 

16.8 

2 5.5 

17 1.2 

29816. 

C. 663 

1 

243 

60 

5.1 

l 56.0 

-5.0 

33 .6 

23 .8 

19 . 4 

20 .4 

171 .4 

29659 . 

0.865 


2 04 

60 

5 • L 

156.0 

-5.0 

32.7 

23.7 

2C .2 

20. 4 

171.5 

39499. 

C.888 


244 

50 

5 . 1 

156 .0 

-5.0 

32.7 

22 .7 

15.4 

2 C • 1 

17 1.7 

25236. 

c. eee 

— 1 

2C5 

5 C 

6*7 

1 57.0 

-6.5 

33.7 

23 .7 

15 .5 

30.7 

171 .9 

39177. 

0 .80 1 


245 

4Q 

6.7 

164.0 

-6.6 

32.6 

23.5 

15.8 

3 1.8 

172.1 

2 5 C 14. 

C.5C5 

— i 

2 C 7 

1 C 

7.2 

1 55 .0 

-7.0 

33.7 

23 .5 

19 .2 

21.7 

172 .3 

30825. 

0 .502 

1 

2 4 6 

60 

6.7 

16 1.0 

-6.6 

32.6 

23.5 

15.2 

2 1.0 

172.4 

38648. 

C • 9 1 1 

— 1 

2C7 

70 

7.2 

157 .0 

-6.5 

22 .6 

22 • £ 

2 C • £ 

2 1.2 

1 72 .6 

28470 . 

0.515 

1 

24 7 

2 C 

7 . 2 

i £7. C 

-6.5 

33.6 

23 .7 

19.3 

30 .5 

1 72 • 0 

38292 . 

0.945 

— 1 
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Table C.4. Sample of corrected meteorological data, NASA CV-990 
See text for definition of terms. 



In examining the times of each of the frames as reported on tape, it was 
found that although the average spacing of frames in time was just slightly 
over 1 s, the resolution of the time measurements was limited to whole 
seconds. This problem is seen in Table C.l in the column marked "TIME", 
although no digits of time beyond the decimal were printed in this particular 
table, such digits of the raw time data were examined and always found to be 
zero. This "time quantization" was capable of introducing a highly patterned 
error into the aircraft positions during the time integrations, for instance, 
aircraft positions would be in error by amounts ranging up to 120 m every 
fifth frame if frames were taken every 1.2 s and the aircraft were flying at 
150 m s' 1 , because the reported times would be in error by as much as 0.8 s. 

As a consequence of these time errors, it was possible for velocity errors as 
large as 1 m s - * to be generated by the objective analysis scheme if a frame 
of data were missing in an area having significant velocity gradients; such a 
situation was actually observed to occur occasionally in the data collected 
near gust fronts (see Section 3.5.1). Thus, it was deemed necessary to con- 
struct a special algorithm which could produce more accurate estimates of the 
time at which each data frame was collected. 

The design of this time correction algorithm was nontrivial in that the 
lidar duty cycle between frames was not constant for more than about 20 or 30 
seconds at a time on average, and often much less. This made it necessary for 
the algorithm to scan backward and forward locally through the data in each 
run until best straight line fits could be made locally to the recorded time 
values. The ultimate aim of the algorithm was to establish corrected times 
spaced as evenly (locally) as possible which, when truncated by dropping the 
fractional parts of each second, would yield the original recorded times. A 
sample of the raw recorded time data and the corrected time data output by the 
algorithm is shown in Table C.5, which consists of a "before corrections" and 
an "after corrections" section, again for the first 62 frames of run 2. The 
quantity "TSEC" is the relative time measured in seconds from the beginning of 
the run, and is calculated in the "before" section by merely comparing the 
values of the "TIME" array with the initial value of "TIME". In the "after" 
section, however, TSEC is the basic quantity; it has been computed in the time 
correction algorithm and then used to create a corrected version of the TIME 
array. Note the significant improvement in time resolution with much smoother 
rates of change of the time parameters. The small discrepancies between the 
digits to the right of the decimal in the final TSEC and TIME arrays is a 
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Table C.5. Samplp of raw and corrected time data. See text 
for further details. 
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consequence of single-precision arithmetic round-off of the TIME values. This 
round-off error has no effect on any other computations, since only the TSEC 
arrays were used in estimating relative times in the analyses. 

APPENDIX C.3 INS GENERATED ERRORS IN AIRCRAFT POSITION AND GROUND VELOCITY 

Low-frequency frame biases (i.e., changes in average radial velocity 
error from frame to frame) enter the lidar data as a result of time-varying 
errors in the output from the aircraft INS. These errors cause incorrect com- 
ponents of aircraft ground velocity to be removed from the raw radial velocity 
estimates during the operation of referencing the Doppler velocity to ground- 
relative coordinates (Appendix A.l). Components of the ground velocity errors 
along the lidar lines-of -sight then contaminate the PPPP estimated radial air 
velocities and, after synthesis of the wind field, ultimately show up as any 
of a variety of vector error patterns superimposed on the true wind field. 

The problem may be either minor or serious depending on both the amplitude of 
the INS errors and on their time rate of change; the INS error amplitude is 
found to influence the amplitude of the error introduced into the calculated 
wind field, while the time rate of change of the INS error is found to 
influence the spatial scale of the error patterns produced. 

Nadir photographs taken at 10 s intervals during portions of the Severe 
Storms flight could be used to establish reasonably accurate aircraft posi- 
tions during selected time intervals, allowing at least rough estimation of 
the magnitude and direction of the INS position and velocity errors for por- 
tions of the flight. Although at first glance it might seem that the nadir 
photographs would resolve conclusively any controversy about the aircraft's 
exact position, it was found that sufficient uncertainty existed in the photo- 
graphically dervied positions to render their interpretation difficult. These 
uncertainties (caused by the relative absence of landmarks easily identified 
in the photos), the relatively coarse scale of the U.S. Geological Survey 7,5 
Minute Topographic Quadrangle maps which were compared to the photos, the 
relatively long time period between photos, and the sensitivity of camera 
orientation to aircraft pitch and roll perturbations, produced random uncer- 
tainties of about 30 m (standard deviation) in the photographically derived 
positions. These uncertainties imply lateral aircraft velocity uncertainties 
of about 30 m each 10 s, or about 3 m s"* from photo to photo, far too large 
to be of help in studying any short-time-scale INS errors in aircraft velo- 
city. 
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Nevertheless, the nadir photo-derived positions were useful in studying 
the INS velocity errors over larger spatial scales (i.e., over many frames). 
Such studies were made for those runs where the photos were available, runs 2- 
4, and the results indicated the INS-derived ground velocity errors were often 
as large as 2 m s" 1 when averaged over each run. 

The technique employed to estimate the INS-derived aircraft position and 
ground velocity errors was made complicated by the differing sampling rates of 
the photos and the INS reports on aircraft positions, the variability of the 
position sampling rate, and the inherent non-straightness of many of the 
flight track segments. There are two ways to define the INS-derived position 
error in these data: one involves interpolating the INS positions at the 
times of lidar frames to the much wider spaced times of the nadir photos, then 
subtracting the x and y coordinates of the two sets of positions to obtain 
position error estimates at the nadir photo times; the other involves inter- 
polating the nadir photo positions to the much closer spaced times of the 
lidar frames, followed again by subtraction of coordinates to obtain position 
error estimates at the lidar frame times. 

Attempts to use either technique to determine the best estimates of posi- 
tion error required smoothing of the photo derived positions. After some 
experimentation with various smoothing techniques, a method based on 
detrending followed by Fourier filtering was finally developed which seemed to 
give satisfactory results. In this method, the position errors at the times 
of the lidar frames are detrended such that the remaining time series of 
"residual" position error are forced to be periodic within the data sequence; 
these residual values are then Fourier analyzed, then subsequently Fourier 
synthesized using only harmonics 0, 1 and 2. Addition of the originally 
diagnosed trend back to these smoothed residual position error data produced 
the smoothed sequence of INS position errors (Fig. C.2). These smoothed posi- 
tion errors were then finite-di fferenced with respect to the time intervals 
between lidar frames in order to produce estimates of the aircraft's ground 
velocity error occurring in each frame. Additional Fourier filtering of the 
velocity error results was needed because of new noise introduced during the 
differencing process as a result of the non-constancy of the time intervals 
between lidar frames. 

Although a number of simplifying assumptions were made in these INS error 
analyses, the results of the technique seem to show good fidelity to the 
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Figure C.2 INS derived position error vectors for run 4. Also shown are the 
aircraft's ground velocity errors throughout the run. 
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original photo-derived position error data. The vector error in INS position 
is seen in Fig. C.2 to average about 1.6 km south-southeast of the "true" 
(photo-derived) position, with variations of about 200-300 m around this 
average. However, it is just the time rates of change of these position error 
variations which are of interest, because they define the INS-derived ground 
velocity errors. As noted at the bottom of the figure, the position error 
vector shifted, from its initial coordinates, toward the east-northeast (81.4 
deg) by some 171 m during the course of the 255.1 s-long data run, implying an 
average INS velocity error of 0.67 m s - ^ from a direction of 261.4 deg during 
the run. The fact that the aircraft made a nearly 30 deg left turn during the 
run is shown by the ground velocity vectors which have been added to the dia- 
gram; the effects of this turn on the quality of the INS data will be shown to 
be quite dramatic. 

In order to analyze the ground velocity errors in more detail, the x- and 
y-components of the position and velocity errors were computed and studied 
separately. In Figs. C.3 and C.4 are depicted the x- and y-position errors at 
the times of lidar frames. 

The Fourier-filtering was performed separately over the three main parts 
of the time series of position data for run 4 (Fig. 3.2): (a) the initial 

straight flight leg prior to the left turn; (b) the left turn (as identified 
from the time series of aicraft roll angles), and (c) the final straight 
flight leg following the left turn. Following reconstruction of the Fourier- 
filtered position time series, additional smoothing using a simple three-point 
filter was applied to remove cusps which occurred at the junction points 
between the three parts of the run 4 data. 

The final filtered sequences of x- and y-position errors used to compute 
the ground velocity errors in run 4 are given in Figs. C.5, C.6. Note the 
removal of the high-frequency noise and the retention of the basic pattern 
that had been evident in the unfiltered data. In Figs. C.7 and C.8 are given 
the ground velocity estimates which were deduced from the time variations of 
the data in Figs. C.5, C.6. Fourier-filtering techniques similar to those 
applied to the position data have been used. Note that it is apparent both 
from the velocity curves and also from the position curves that the ground 
velocity errors may exceed 2 m s“* during straight flight, 10 m s“* during 
turns, and may vary with time. The time variability of these errors has par- 
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FRAME NUMBER 



Figure C.5 Fourier-filtered version 
of x position error data in 
Fig. C.3. Only harmonics 0, 1 and 2 
plus the observed trend were used to 
create these position data. 


Figure C.6 The same as Fig. C.5 
except for the y position errors. 
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Figure C.7 The x component of 
ground velocity errors estimated 
using data in Fig. C.5. After 
differentiation with respect to the 
uneven frame time intervals, further 
Fourier filtering was used. Mean 
error = 0.70 m s . 
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Figure C.8 The same as Fig. C.7 
except this is y component of ground 
velocity errors estimated using data 
in Fig. C.6. 
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ticularly important consequences for the credibility of the 1 idar-measured 
winds as is discussed in Section C.4. 

Note also the appearance of large ground velocity errors at the beginning 
of the run; this pattern could be the result of the method of detrending the 
raw position errors, which did not employ any least-squares straight-line 
fitting. The detrending technique used in this study, while insuring pure 
periodicity of the residual data within the data domain, is vulnerable to the 
effects of large random data errors at the end point* However, close examina- 
tion of the raw data associated with the beginning of run 4 failed to reveal 
any noteworthy irregularities in the photo position estimates that might 
account for the large ground velocity estimates obtained there, rather it 
appeared that these errors came directly from a large-scale error trend which 
was clearly apparent in the raw data. In summary, some skepticism is war- 
ranted for the values of ground velocity error obtained near the endpoints of 
data sequences, although for the data shown here evidence also exists which 
tends to support the values calculated. 

The estimated ground velocity errors are well outside the manufacturer's 
specifications of a 3-sigma velocity error of 0.5 m s" 1 , and are capable of 
introducing significant bias into the lidar radial velocity measurements with 
respect to the ground. Substantial further investigation, including addi- 
tional field measurements, will be needed in order to obtain more definite 
information regarding the magnitudes and variability of the ground velocity 
errors on time scales that affect the wind estimates. The effect of these 
time-dependent errors on the wind vectors retrieved by the lidar system are 
demonstrated in Section C.4. 

C.4 EFFECTS OF NAVIGATION ERRORS ON DOPPLER LIDAR DERIVED WIND FIELDS 

Important to the wind field estimate accuracy is the velocity error asso- 
ciated with the time rate of change of the position error. Samples of esti- 
mates of the two components of the velocity errror were presented in Figs. C.7 
and C.8. The most conspicuous feature of these two plots was the large but 
relatively brief y-velocity error which occured between frames 132 and 153, 
during a left turn made by the aircraft. This velocity error attained a 
magnitude of more than 10 m s’ 1 for a portion of the turn, a magnitude which 
was roughly confirmed by inspection of the behavior of the y-position error in 
Fig. C.6. Most of the velocity error was confined to the y-direction during 
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the turn, and since the aircraft track was nearly eastward at the time, it 
appeared that most of the ground velocity was due to inaccurate measurement of 
lateral accelerations. It is clear that INS errors are especially likely to 
occur in the vicinity of aircraft turns. This suggests strongly that biases 
might be present in the lidar-measured winds from the towering cumulus circum- 
navigations (runs 5-12), which featured numerous sharply banked turns. Unfor- 
tunately nadir photo documentation of the aircraft's movement was not avail- 
able during the cumulus circumnavigations due to a temporary malfunction of 
the camera; therefore the performance of the INS during these contorted 
maneuvers could not be assessed. 

The exact manner in which these large errors are introduced into the INS 
output is not yet known, although it is possible that internal filtering 
processes employed in the INS computation sequence might be capable of pro- 
ducing large transient errors during turns and accelerations. Further inves- 
tigation will be needed in order to isolate the true cause or causes of the 
errors. 

In order to establish what sort of wind vector errors might be introduced 
into the lidar analyses as a result of neglect of or incorrect compensation 
for the INS errors, artificial data or "measurements" consisting only of the 
components of estimated ground velocity errors along the lidar lines-of -sight, 
were synthesized to generate wind vectors. Hence the raw velocity data used 
in these analyses consisted of "data" which were strictly independent of range 
for each lidar frame. The results, a sample of which is shown in Fig. C.9 for 
run 4 data, reveal the presence of a number of vortical and divergent "flows" 
generated by patterns in the time-varying ground velocity error interacting 
with the overlapping of fore and aft lines -of -sight at various ranges. Some 
of these error vectors have magnitudes of 3-4 m s* 1 , quite significant in view 
of the fact that the actual wind field itself was only about 10 m s"* in 
strength for much of the data examined. 

It is most important to realize that the magnitudes of the wind field 
errors caused by these time-varying errors in radial velocity are not limited 
to the actual magnitudes of the ground velocity errors which caused the radial 
velocity errors. Rather, the components of the ground velocity errors along 
the lidar 1 ines-of-sight are "distributed" throughout the data collection area 
as the aircraft continues along its flight track, only to be combined with 
other erroneous radial velocity components resulting from later or earlier 
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Figure C.9 Wind field corresponding to INS velocity errors from Figs. C.7 and 
C.8. Note the presence of vortical and divergent patterns resulting from the 
time-varying INS errors. 


errors. If the error vectors should happen to lie outbound along a fore line- 
of -sight at one moment, to be followed a moment later by an error of the same 
magnitude lying inbound along an aft line-of-sight, the error in the wind 
field vector estimate obtained from combining these two radial velocity errors 
would be almost three times as large in magnitude as either of the original 
ground velocity errors. 

Although the error in wind vectors of Fig. C.9 was obtained using esti- 
mates of ground velocity error which varied relatively smoothly with time, it 
is clear that any time-varying errors introduced into the radial velocity data 
have the potential to generate wind error. Generally, the vector error pat- 
terns will have a spatial scale which corresponds with the time scale of the 
radial velocity errors, and will have magnitudes which correspond roughly to 
the magnitudes of the radial velocity errors. Lingering errors caused by 
incorrect estimation of drift angle biases could, for instance, produce small- 
scale, weak vortical and divergent patterns in the wind field analyses; at the 
same time, ground velocity errors which are larger in magnitude but more 
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slowly varying in time could produce a superimposed pattern of larger scale 
vorti cities and divergences. 

Thus it would appear that any wind perturbations having magnitude of 
3-4 m s -1 or less obtained from the Doppler data analysis must be viewed with 
caution and skepticism, as perturbations of that magnitude could have been 
caused solely by any of the various forms of bias errors in the radial velocity 
measurements. Fortunately, many of the gust front observations made during the 
Severe Storms Experiment contained such strong wind field features that even 
large errors of the type discussed here could not obscure the real flow 
pattern. 



APPENDIX D 

INSPECTION AND VERIFICATION OF DOPPLER SPECTRAL MOMENT DATA 


Inspection of listings of the data fields from the tape furnished by NASA 
showed that the NASA-corrected velocities and intensities contained occasional 
inconsistencies that were suggestive of data processing errors. The raw data 
fields, however, seemed to be unaffected by these errors. 

Table D.l shows a sample of the raw tape contents from the first frame of 
run 2, featuring the corrected data supplied by NASA as well as the raw 
data. The columns marked "XBEAM" and "YBEAM" give the X and Y coordinates of 
each resolution volume center relative to the start of the run. The raw lidar 
velocities are given under "RAWVEL", while the corrected and smoothed lidar 
velocities furnished by NASA are given under "CORVEL". All velocities are 
given in m*s“*. Raw lidar intensities are given in units of dB above noise. 
Spectral widths are given in logarithmic coded form under "WIDTH". Conversion 
of width values to m s"* format is described below. 

Of particular note in Table D.l is the susceptibility of CORVEL to large 
errors; see resolution volumes ("BINS") 12 to 14 and 16 for examples. In 
these same resolution volumes the parameter RAWVEL is seen to be much better 
behaved; of the first 18 resolution volumes only the first contained noisy 
data, and this was expected based on the fact that the transmission of the 
signal was still in progress during the time the first resolution volume was 
being sampled. 

Values of CORAMP seem fairly reasonable, except that the steady decrease 
seen from resolution volumes 2 to 14 suggests that a range-squared correction 
may not have been incorporated. The coded values of width appeared to be 
quite consistent and credible; note the near constancy of width through 
resolution volume 13, followed by a sudden degeneration to large noisy values 
at larger ranges. Note also the fact that the RAWVEL array retains consis- 
tency out to larger ranges than width; this is a consequence of the way in 
which spectral width is estimated in the poly-pulse-pair processor (PPPP) 

(Lee, 1980) and is not sufficient reason to consider the velocities lying in 
the region of large width to be suspect. Furthermore it was ascertained that 
the values of XBEAM and YBEAM (Table D.l) were not consistent with ranges to 
the resolution volumes (bins) deduced from the measured time positions of the 
sampling gates in the PPPP (McCaul , 1985). 
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Table D.l, Sample of raw lidar moment data from frame 
1 of run 2. See text for further details. 



Inspection of data such as that shown in Table 0.1 made it clear that all 
data preparation for meteorological interpretation would have to be accom- 
plished starting with the raw data and that resolution volume coordinates 
would have to be recomputed. Furthermore, raw data also made it clear that 
much raw data at longer ranges from the aircraft were noisy and should be 
ignored when creating disk files of raw data to be used in later analyses. 

The range beyond which data degenerated to noise was found to depend on air- 
craft altitude, and hence probably was controlled by levels of ambient mois- 
ture. For runs 2 through 8, only the first 28 range bins out of 45 available 
were retained for subsequent analysis; for the higher altitude runs 9 through 
12, the first 40 range bins were retained. 

In order to prepare the data for meteorological interpretation, it was 
necessary to range-normal ize the raw intensity data to produce estimates of 
reflectivity, convert the width data to m*s"* format, and remove any "frame 
bias" from the measured Doppler (radial) velocity data. This frame bias con- 
sisted of velocity errors which were approximately constant within each frame 
of data and were associated with incomplete or erroneous documentation of 
small scale variations in aircraft attitude and motion.. The removal of frame 
bias from the velocities in each frame is done only after point editing, and 
is discussed in Section D.3. 

D.l PREPARATION OF REFLECTIVITY FIELDS 

Range normalization of the raw intensity furnished by NASA was accom- 
plished using a two-step algorithm and is fully discussed by McCaul (1985). 
Therefore, only an outline is given here. The first step was to determine for 
each frame the background noise level of the intensity field. This noise 
level was found by averaging the intensity from the most distant 20 resolution 
volumes, with any intensity features rising more than 3 dB above the minimum 
of all resolution volumes excluded from the averaging. In cases where the 
minimum intensity in any frame occurred at a range closer than the outermost 
20 resolution volumes, the bin number of that minimum value was used as the 
inner range bin limit in the search through maximum range for values to be 
used in computing the noise level for that frame. The exclusion of any high 
intensity features from the noise level calculations was made necessary by the 
occasional occurrence of highly reflective cloud edges in the most distant 20 
resolution volumes. The noise level was determined for each frame because 
small random changes of 1-2 dB were noted in the background noise values 
during the course of the typical data run. 
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The second step in the range normalization was to subtract the noise 
level from all recorded values of intensity and to multiply the result by the 
sum of the squares of the resolution volume range and of a "Rayleigh Range" 
factor related to laser beam divergence (Murty and Bilbro, 1978). An exponen- 
tial range-dependent attenuation function is included in the Murty and Bilbro 
formulation. 

It was appropriate to estimate the attenuation constant k by linearly 
fitting the range-dependent decibel variations in intensity. Although the 
actual intensity attenuates exponentially with range if k is constant and 
independent of range, the intensity, expressed in decibel units, will decrease 
linearly with range. 

Only relative intensities were computed because the information regarding 
emitted laser power was not available for each frame. An example of estimated 
relative reflectivity fields computed from raw intensity data using the above 
described method is shown on Figure D.2 whereas Fig. D.l shows the raw inten- 
sity data. The raw intensity profile from frame 12 of run 9 (Fig. D.l) is 
seen to settle down to a noise-level value beyond a range of 5 km from the 
aircraft, after showing a large jump at a cloud edge which was located at a 
range of about 4 km. The cloud-strike reflectivities (Fig. D.2) are seen to 
stand out clearly from the reflectivities in the cloud free air at near 
ranges. The cloud attenuates the lidar pulse so that no echoes are received 
beyond 4.8 km range although cloud may be present there. 

The results of all the revisions documented above are visible in Table 
D.2, which represents the corrected version of Table D.l, with the decoded 
versions of Doppler spectral width (McCaul 1985) listed under the "WIDTH” 
column. Note that not only are spectral widths converted to m*s"* format and 
reflectivities estimated, but also the positions of the lidar resolution 
volumes with respect to the aircraft have been modified so as to be consistent 
with resolution volume positions deduced from timing measurements. 

D.2. THRESHOLDING AND EDITING OF DOPPLER MOMENT DATA 

The next stage in preparing the data for analysis involved "point- 
editing", or removal of noisy data point-by-point based on detailed inspection 
of the data. Point editing was performed on velocity, intensity, and width 
fields separately. The first step was to apply a "thresholding" algorithm 
which could set to "missing" the large numbers of obviously useless values 
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Table D.2. Sample of corrected lidar wind, intensity and 
width data for data in Table D.l. See text for further 
details. 


j 


140 



tNIENSltr na 


RRNUE SIN MS 
<? 1:5 18 


24 


21 


30 




Figure D.l Typical raw intensity profile 
for a lidar frame containing a cloud 
strike (frame 12 of run 9). 


Figure D.2 The profile of relative 
reflectivity estimated from data in 
Fig. D.l. Note that reflectivity 
is uniform in range until the lidar 
pulse strikes the cloud at 4 km. 
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encountered at ranges beyond the lidar's functional limits. The editing of 
width and intensity data required no special threshold routine because all 
valid (i.e., non-spurious) reflectivity estimates were found to fall between 
the fixed thresholds of 20 and 80 dB above noise and all valid Doppler spec- 
tral widths were found to fall between 0 and 6.0 m*s'^. The velocity thres- 
hold limits had to be established individually for each data run; for a list 
of the values used, consult Table D.3. The spurious velocities which were 
removed via the thresholding method were easily identified in listings of raw 
data because of their relatively large magnitudes; examples are visible in 
Table D.l under the RAWVEL (Raw Velocity) column, beginning at resolution 
volume (BIN) 19. 

Further editing was often needed to remove spurious point values. To 
simplify this sometimes tedious process, a graphical routine was designed 
which displayed each of the fields to be edited in a way which would facili- 
tate finding spurious data quickly. Samples of the graphical output are shown 
in Figs. D.3-D.6, for data from run 10 of the Severe Storms experiment. 

In Fig. D.3 the raw radial velocities in each fore and aft frame are 
depicted in approximately their true relative positions in space; the 
"vectors" representing the radial velocities are drawn at an angle with 
respect to the line-of -sight of each lidar frame for clarity and ease of pat- 
tern observation. 

Fig. D.4 shows both the lidar intensity and spectral width for the data 
of run 10. The circles which represent the locations of the data points in 
space have diameters proportional to relative reflectivity, while the line- 
weight with which each circle is rendered is roughly proportional to spectral 
width. The band of large circles seen at mid-range in Fig. 0.4 thus repre- 
sents a highly reflective target, in this case the edge of a cumulus congestus 
cloud. The lightweight line with which the large circles are drawn indicates 
small measured spectral widths, and presumably, good definition of mean 
Doppler velocity. 

Figures 0.5 and D.6 depict the observed distributions of observed radial 
velocities for each fore (Fig. D.5) and aft (Fig. D.6) lidar frame. The velo- 
city scale used in these figures runs from -20 m*s -1 at the bottom of the plot 
to +20 m*s"l at the top. In both figures valid meteorological data appear as 
swaths of loosely distributed points near the centers of the plots, with 
spurious returns arranged in swaths at the edges of the plots. In Fig. D.6 
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Summary of Characteristics of Lidar Data Runs 
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VMIN = Minimum radial velocity used in thresholding 
VMAX = Maximum radial velocity used in thresholding 
(JADV = Zonal motion component used in data advection 
VADV = Meridional motion component used in advection 

TSTART = Time (GMT) data collection began in the run 

TSTOP * Time (GMT) data collection ended in the run 

OATIME = Time (GMT) of objective analysis of the run 





Figure D.3 Graphical output for 
editing data: raw lidar radial 

winds from run 10. Wind data 
from fore frames are shown 
above; aft data shown below. 
Along-track distance scale is 
expanded for clarity. 





> o£pOO OO»OQPOOO 

oooo -cxjdooo QMO fleacxrocccQOOcoo o o o o 
o oooooo 0©0©000 0^0000 o o 

OOO CD-OGDO (JD(33)©<2DC1DQD<I)<30<3DC3DC2) O O O 
O O 0(2)0 <JD OGDCDCIDGD® ®(2>®<5>(2><2)0 O O 

o o oo oo ©<0<n>©O(ii><n>(Q)®®(i3)<n) o o 

O OODGDO- OQOCDOOGDODODOXjOGDQDOOCDO O 

OOOO oooooooooooooo 

OOOOOOO 00,0000000000^000000000^00000 ^ 


FRAME NUMBER 



0 

10 

DB 

0 

10 

-20 

OB 

O 

20 

•30 

OB 

0 

30 

-yo 

OB 

O 

40 

-50 

OB 

o 

50 

-60 

OS 

0 

60 

70 

OB 

O 

70 

-80 

08 

o 

0,0- 

•l, 0 

M/S 

o 

a. 0- 

•2.0 

M/S 

0 

2,0- 

■3.0 

M/S 

o 

3.0- 

y.o 

M/S 

o 

4.0- 

5.0 

M/5 


Figure D.4 Graphical output (run 10) for editing data: relative reflectivity 
is proportional to circle diameter, and spectrum width values are scaled in 
proportion to line weight. Small dashes indicate low-intensity, large-width 
returns. 
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Figure D.5 Graphical output for editing velocity data: velocity distribution 
is for the same fore frame data presented in Fig. D.3. Data points near top 
or bottom margins are due to moding in the lidar transmitter. 
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Figure D.6 Same as Fig. D.5 except for the aft frame data 
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the meteorologically useful data is the swath of positive velocities whereas 
the swatch of negative velocities is associated with resolution volumes at 
unusably long rangaes (Fig. 0.3) as is the data swatch near +17 m s"* in Fig. 
D.5. Velocity estimates at long ranges, where echo intensity is less than 
noise, should be uniformly distributed across the Nyquist interval. These 
coherent velocity swaths in echo free regions might be caused by leakage of 
coherent signals into the PPPP. 

Thus, with the help of graphical displays, it was possible to examine the 
data from each resolution volume of each lidar frame and to determine quickly 
which data points were consistent with neighboring data points. Such a visual 
examination was made for every data point in every run, and if an inconsistent 
data value was found, it was included in a list of "bad" points contained 
directly in the text of a point edit subroutine whose function was to replace 
such data points with a "missing" data parameter value, usually -999.99. On 
some occasions it was necessary to delete entire frames from the data set due 
to noise contamination. 

D.3. DIAGNOSIS AND REMOVAL OF RADIAL VELOCITY BIASES 

The next step in preparing the data for meteorological interpretation 
involved deletion of velocity biases which contaminated radial velocity data 
in each frame. Correction of the measured radial velocities for this "frame 
bias" was a difficult and troublesome task. Lee (1982), McCaul (1985), and 
other investigators using the NASA data have described various approaches 
taken. 

Lidar frame bias manifests itself in the tendency for the variability of 
the measured velocities to be much larger from frame to frame at a given range- 
than from resolution volume to resolution volume within a given frame. The 
underlying physical cause of the bias is not completely understood, but 
apparently is related to inadequate measurement of changes in aircraft atti- 
tude and relative airspeed, especially those associated with quasi -periodic 
fugoid motions and cyclic feedback effects in the aircraft flight control 
system, and to larger-than-expected delays in recording certain INS-deri ved 
data such as drift angle (Lee, 1982). 

Lee (1982) found evidence in a detailed study of lidar data taken in 1981 
in California that the frame biases were most likely due to delays in 
reporting INS-measured drift angle values to the data recording system. Lee 



(1982) examined various possible delays in drift angle and found a remarkable 
correlation between the excursions of the range-averaged radial velocities and 
the excursions of the lag-2 differences in measured drift angles. An attempt 
was made in this study to duplicate Lee's findings by means similar to those 
he employed. However, for the runs examined in the 30 June data, no clear cut 
relation between velocity excursions and drift differences emerged. 

In general, the large variability of the true wind speed in the convec- 
tive situation found on this day prevented the discovery of the specific 
physical mechanisms responsible for the biases. However, examination of the 
behavior of the velocity in all the data runs suggested that the magnitude of 
the problem was sensitive to the local turbulence of the flow near the air- 
craft and to the relative orientation of the airspeed and ambient wind 
vectors. 

Based on a drift angle study conducted by McCaul (1985), and the absence 
of any clear-cut evidence as to the physical mechanism responsible for the 
frame bias, it was decided to handle the problem by means of filtering only. 

Two important additional causes of low-frequency velocity errors (i.e., 
frame bias) were found in the present data: ground velocity bias associated 

with data recording delays, and ground velocity errors produced by the INS 
(appendix C.4). Both of these error sources were found to be capable of 
causing radial velocity errors of 2 m s"* or more, and are thus comparable in 
importance to the high-frequency drift angle delay errors already discussed. 
Unfortunately, neither the errors due to delays in recording ground velocity 
(discussed in the next paragraph) nor errors generated by the INS could be 
eliminated by using simple filtering techniques because their spectral content' 
was too similar to that of the true wind field. These low-frequency frame 
biases therefore had to be estimated using independent means. 

A portion of the radial velocity errors (frame biases) were found to 
result from delays of recording INS-derived ground velocity during periods of 
aircraft acceleration or deceleration. The transmission of the ground 
velocity errors "EGS" through the data collection chain for an assumed calm 
wind field and accelerating aircraft is depicted in Fig. 0.7. From the figure 
it. is clear that an erroneous wind vector equal in magnitude and direction to 
the ground velocity error vector EGS is "measured" by the lidar system as long 
as the aircraft acceleration is maintained constant. Addition of a constant 
nonzero wind field to this hypothetical situation will not alter the error. 



TRUE GROUND SPEED 



Figure D.7 Geometrical explanation 
of the effects of ground velocity 
errors (due to delays in recording 
data) on the radial velocity 
computation. 


However, depending on the relative magnitudes of the ambient wind vector and 
the error EGS, and on the orientation of the ambient wind vector with respect 
to the aircraft track, an apparent rotation of the lidar-derived wind vectors 
from their true orientation will be observed. Thus for an ambient wind 
blowing from right to left across the aircraft track and a situation where the 
aircraft is accelerating, the data synthesis will produce wind vectors whose 
tips are rotated rearward of their true orientation with respect to the air- 
craft; for a decelerating aircraft subject to the same wind, the synthesized 
wind vectors will be rotated forward of their true aircraft-relative orienta- 
tion. Some complication of the pattern will occur at the endpoints of the 
aircraft acceleration, where the fore and aft lidar beams are unequally 
affected by delays in reporting the ground velocity. 

Although the delays in reporting the aircraft's ground velocity could not 
be estimated in detail, it is believed that they contribute significantly to 
the error budget of the lidar system for the data reported here. 
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APPENDIX E 

ASSESSMENT OF SCANNER PERFORMANCE 


A series of 13 consecutive lidar beam ground strikes was observed during 
a left turn executed by the aircraft at about 2039 GMT (1439 CST) during run 4 
data collection. Through study of this sequence of ground strikes, it was 
possible to obtain some insight into how well the scanner system was opera- 
ting. Figure E.l shows plots of the sequences of corrected (see Appendix C.l) 
roll angles (curve in top plot) and range bins in which ground strikes, iden- 
tified by peaks in lidar intensity, were observed (boxes in bottom plot). It 
is apparent from a comparison of the variations of roll angle and observed 
ground strike ranges that the measured values of roll angle corresponding to 
the various lidar frames were being recorded approximately one frame, or about 
1.33 s, too late. When roll angle and ground strike range data are appro- 
priately paired (for instance, roll angle from frame 139 with ground strike 
range from frame 138), a rough quantitative assessment of two aspects of beam 
pointing performance can be made. 

First it can be seen that, in accord with design specifications, no 
ground returns were observed for roll angles more negative than -9°; the roll 
angle of -8.4° reported in frame 134 might have been expected to produce a 
ground return in range bin 20 of frame 133 (accounting for the 1-frame delay 
in reporting of roll) if the scanner were not compensating properly for roll 
angles between +9° and -9°, but because no ground return was found in frame 
133 there was no evidence that the beam scanner was failing to compensate for 
roll angles. This of course does not conclusively rule out possible beam 
pointing errors in near-zero roll situations, where the scanner operates in a 
different manner in its attempts to maintain horizontal beam trajectories. 
Although a range of 20 bins is considerable, it is not too far for the recep- 
tion of signals from highly reflective targets. 

Additionally, it is possible to predict from a knowledge of the roll and 
scan angles and aircraft altitude the range bin in which a ground strike 
should occur, assuming flat terrain and operation of the scanner according to 
specifications (i.e., return the beam to point along a line perpendicular to 
the aircraft heading and in the plane of the wings) at roll angles greater in 
magnitude than 9°. By comparing the predicted ground return range with the 
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Figure E.l Ground strike ranges and aircraft roll angles as a function of 
frame number. Top curve depicts aircraft roll angle; bottom depicts predicted 
(line) and observed (boxes) range of ground returns. 





observed value, limits can be placed on the magnitude of the beam pointing 
biases reported by Lee (1982), assuming the biases are always present in the 
scanner performance regardless of the roll angle. 

Such a comparison is made in the lower plot of Fig. E.l, where the solid 
curve superimposed on the observed values of ground strike range represents 
the values of predicted ground strike range, taking the 1-frame delay in roll 
angle properly into account. An averaged aircraft altitude of 847 m was used 
in all calculations of ground strike range in order to remove the influence of 
topographical variations, observed to be about 10-20 m in magnitude, from the 
computed range values. In general, the agreement between the predicted and 
observed ground strike ranges is excellent, with all but three frames showing 
exact agreement. Of the three frames not showing exact agreement, two are in 
"error" by only one range bin, and the other lies at the edge of the data 
sequence, where the roll angle magnitude was decreasing rapidly and the ground 
strike range was growing even more rapidly in response. 

For the ground strike ranges shown in Fig, E.l, typically 2-3 km, and the 
aircraft altitudes associated with those ground strikes., 850 m, the appearance 
of a ground strike in a given range bin defines the beam elevation angle ('mea- 
sured relative to the horizon) only to within about a degree. Although the 
range to ground targets is quite sensitive to errors in beam elevation angle, 
the nearness of some for the ground strikes to the aircraft, the large range 
extent of the range bins and the relative sparseness of sampling in time pre- 
vent a more definitive bracketing of the elevation errors. 

In conclusion, it appears that if the elevation errors of the scanner at 
small roll angles are similar to the errors seen at larger roll angles, then 
those errors are apparently limited to about 1. 5-2.0° for both fore and aft 
shots. The biases found by Lee (1982) fall well within these limits. 

Although the elevation error for both fore and aft shots could be zero, it is 
likely that some nonzero elevation error, probably of the same order as that 
observed by Lee, was present in the lidar data. Such elevation angle errors, 
of the order of 1°, could introduce significant errors in wind estimates in 
zones of strong vertical wind shear. 



APPENDIX F 


TESTS USING MODEL WIND FIELD 

It was also found desirable to test the wind retrieval algorithms on a 
simulated monochromatic wind field designed to have properties which would be 
easy to identify on plotter output. To perform this test, starting position 
information from runs 9 through 12 were used to initialize the aircraft 
location at the times given for the beginning of those runs. Then, using 
actually measured aircraft airspeed data from each of the runs and the model 
wind field specified by: 

U(X,Y,T) = 10 + 2*C0S[(2tt/ 5000)*(X + Y - 5T)] (F.l) 

V(X,Y,T) = 10 + 3*C0S[(2ir/5000)*(X + Y - 5T)] (F.2) 

the aircraft was forced to fly through the simulated atmosphere along a track 
obtained by solving the vector triangle of airspeed/heading, ground 
velocity/track and the model windspeed/di recti on for each frame in each of the 
four runs using INS-derived headings and initial positions for each run. The 
solution involved computation of the magnitude and direction of the ground 
velocity vector, because the other two quantities were assumed to be known. 

As the aircraft progressed through the model wind field, sample lidar 
frames were acquired at times specified to be identical to those given in the 
actual real data runs, and "observations" of the appropriate radial component 
of the true wind in the various resolution volumes were taken based on the 
behavior of the model wind field as a function of space and time. The 
simulated observed radial velocities and line of sight angles were then 
subjected to the analysis algorithms described in section 3.4.4. 

For the analyses it was necessary to estimate the phase speed at which 
the analytical wind field features are being "advected". Because the "true" 
phase speed was specified in the form of the model wind field itself (see the 
coefficients of time T in the sinusoids of Eqs. F.l and F.2, i.e., C x = C y = 
5.0 m s" 1 , the phase lines propagated to the northeast at 3.53 m s~* and it 
was possible to study the impact of an incorrect guess of the advection speed 
on the analyzed wind vectors. 
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Figure F.l Wind field synthesized from simulated lidar observations of a 
modeled wind field. Data are adverted using an advection velocity equal to 
the true one (i.e., U a = V a = 2.5 m s -1 ). 
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Figure F.2 Perturbation wind field of Fig. F.l (mean of wind shown on Fig. 
F.l subtracted). 
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Figure F.3 Velocity vector differences between the true perturbation wind 
field and that synthesized from the simulated observation. 
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Figure F.6 Same as Fig. F.3 except U a = V a 10 m s' 1 . 
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The results of the analysis for an assumed advection which matched the 

true "advection" (i.e., U a = V g = 2.5 m s -1 ) are given in Fig. F.l. It is 

clear that the true wind field pattern is quite well depicted, although there 
are some errors near the edges of the analysis area which can be easily seen 
in the perturbation wind field pattern shown in Fig. F.2. Note that the 
errors are largest along the edges of the field as can be seen in Fig. F.3 
which shows the differences in the true perturbation wind field and that syn- 
thesized wind from simulated observations. These errors are due to the inade- 
quacies of the Cressman interpolation scheme when data are not uniformly 
distributed over the interpolation volume. Note also the smooth blending of 
the retrieved vectors in areas where two runs of data overlap in space. 

When the estimated advection vector is set equal to the mean wind 

(U a = V a = 10 m s" 1 ) in the layer, the results, shown in Figs. F.4-F.6 are 

significantly different. The basic pattern of the wind vectors is still quite 
accurately represented, but some inconsistencies arise in areas of overlap 
between two data runs. Note the displacement of data areas by comparing Figs. 
F.l and F.4. In Fig. F.5, displacements are noted between the convergence 
lines retrieved from the synthesis of data from the different runs. The dif- 
ference fields shown in Fig. F.6 are considerably larger than in the case when 
the correct advection velocity is used. These inconsistencies are clearly the 
result of the improper shift of the data points in space caused by the advec- 
tion correction algorithm having to use an erroneous advection velocity. 
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